R语言for循环中简单向量存储问题排查

我很少使用for循环，在实验处理随机化任务中遇到一个简单问题：实验有4种处理，持续12周，需将4种处理随机化3次，生成样本的每周处理分配表。我需要通过for循环为45个样本批量生成处理设计数据集，但循环内的存储逻辑无法正确执行，单独运行循环内代码正常。

原代码：

treatments = as.numeric(c(1:4))
sample_reps = 45

store = c()
for (i in 1:sample_reps){
  weeks1_4 = sample(treatments)
  weeks5_8 = sample(treatments)
  weeks9_12 = sample(treatments)
  
  set = c(weeks1_4,weeks5_8,weeks9_12)
  store[i] = rbind(store,set)
}

运行时出现警告：

In store[i] <- rbind(store, set) :
  number of items to replace is not a multiple of replacement length

期望输出示例：

[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
  1    3    4    2    4    2    3    1    2     3     1    4
  4    2    1    3    2    1    3    4    3     4     2    1
...etc

• 初始化错误：store = c() 创建的是一维向量，但我们需要存储的是45行12列的二维数据结构（矩阵或数据框），向量无法容纳多行多列的内容。
• 赋值逻辑错误：store[i] = rbind(store, set) 试图把rbind生成的多行矩阵赋值给向量store的第i个元素，这不符合R的赋值规则——向量每个位置只能存储单个元素，而rbind的结果长度远大于1，因此出现警告。

方案1：预分配矩阵（推荐，效率更高） #

提前创建指定维度的矩阵，避免循环中频繁修改数据结构：

treatments = 1:4  # 1:4本身就是数值型，无需as.numeric转换
sample_reps = 45

# 预分配45行12列的空矩阵
store = matrix(nrow = sample_reps, ncol = 12)

for (i in 1:sample_reps){
  weeks1_4 = sample(treatments)
  weeks5_8 = sample(treatments)
  weeks9_12 = sample(treatments)
  
  set = c(weeks1_4, weeks5_8, weeks9_12)
  store[i, ] = set  # 直接为矩阵的第i行赋值
}

# 查看前6行结果
head(store)

方案2：循环合并（简单直观，小数据量可用） #

初始化为空矩阵，每次循环将新生成的处理序列作为行合并进去：

treatments = 1:4
sample_reps = 45

# 初始化为0行12列的空矩阵
store = matrix(nrow = 0, ncol = 12)

for (i in 1:sample_reps){
  weeks1_4 = sample(treatments)
  weeks5_8 = sample(treatments)
  weeks9_12 = sample(treatments)
  
  set = c(weeks1_4, weeks5_8, weeks9_12)
  store = rbind(store, set)  # 将当前序列合并为新行
}

# 查看前6行结果
head(store)

内容的提问来源于stack exchange，提问作者lmbradley