如何用R生成包含0的a与b所有值对组合的完整数据集？

需求：获取变量a和b所有可能值对组合的数据集，其中a的取值范围为0到30（共31个值），b的取值范围也为0到30（共31个值）。

问题：使用如下R代码生成的数据集d1仅包含900条数据（30×30），缺失了含0的组合，无法得到预期的961条（31×31）数据，需获取正确的R代码。

原代码：

a <- numeric()
b <- numeric()

for (i in 0:30) {
  for (j in 0:30) {
    a[i] <- i
    b[j] <- j
    d1 <- expand.grid(y = a, x = b)
  }
}

错误原因 #

• R的向量索引从1开始，当i=0时，a[i] <- i试图给向量的第0位赋值，这是无效操作，导致0这个值从未被存入a、b向量，最终a、b仅包含1-30的30个值，生成的组合自然只有30×30=900条。
• 嵌套循环完全多余，expand.grid本身就是用来生成所有笛卡尔积组合的工具，无需手动循环赋值。

正确代码 #

方案1：直接使用expand.grid（推荐） #

这是最简洁高效的方式，直接传入0-30的序列即可生成所有组合：

d1 <- expand.grid(a = 0:30, b = 0:30)
# 验证数据行数：nrow(d1) 会返回961，符合预期

方案2：修正循环逻辑（仅作原理演示） #

如果一定要用循环实现，需修正索引问题，利用i+1作为向量索引（对应0-30到1-31的位置）：

# 预先初始化长度为31的空向量
a <- numeric(31)
b <- numeric(31)

# 填充a的所有值
for (i in 0:30) {
  a[i + 1] <- i
}

# 填充b的所有值
for (j in 0:30) {
  b[j + 1] <- j
}

# 生成所有组合
d1 <- expand.grid(a = a, b = b)

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