如何将qbinom循环计算结果存入R语言向量？

问题场景 #

需要计算一组不同p参数对应的二项分布中位数（通过qbinom(0.5, 5000, p)实现），并将结果存入向量。原for循环代码无法正确存储结果，输出全0，但直接打印计算结果是正确的。

原错误代码：

list.of.p.parameter <- c(0.32, 0.32, 0.30, 0.24, 0.29, 0.28, 0.25, 0.29, 0.32, 0.27)                             

for (i in list.of.p.parameter){
  my_vector <- vector(mode="numeric", length=length(list.of.p.parameter))
  median.binom <- qbinom(0.5,5000,i)
  my_vector[i] <- median.binom
  print(my_vector)
}

输出全0，但直接打印median.binom能得到正确结果：

for (i in list.of.p.parameter){
  median.binom <- qbinom(0.5,5000,i)
    print(median.binom)
}

结果：

[1] 1600
[1] 1600
[1] 1500
[1] 1200
[1] 1450
[1] 1400
[1] 1250
[1] 1450
[1] 1600
[1] 1350

错误原因 #

1. 向量初始化位置错误：my_vector的初始化放在循环内部，每次循环都会重置向量为全0，之前的计算结果被完全覆盖。
2. 索引使用错误：用p参数的数值（如0.32）作为向量索引，R会将小数索引自动取整为0，而向量索引从1开始，0是无效索引，导致无法赋值，向量始终保持全0。

修正方案 #

方案1：正确使用for循环 #

将向量初始化放在循环外，用整数索引遍历参数列表：

list.of.p.parameter <- c(0.32, 0.32, 0.30, 0.24, 0.29, 0.28, 0.25, 0.29, 0.32, 0.27)

# 提前初始化向量，长度与参数列表一致
my_vector <- vector(mode="numeric", length=length(list.of.p.parameter))

# 用整数索引遍历，seq_along生成1到10的索引
for (idx in seq_along(list.of.p.parameter)) {
  current_p <- list.of.p.parameter[idx]
  my_vector[idx] <- qbinom(0.5, 5000, current_p)
}

print(my_vector)

输出：

[1] 1600 1600 1500 1200 1450 1400 1250 1450 1600 1350

方案2：向量化操作（更高效） #

R的qbinom支持向量化输入，无需for循环，直接传入整个p参数向量即可：

list.of.p.parameter <- c(0.32, 0.32, 0.30, 0.24, 0.29, 0.28, 0.25, 0.29, 0.32, 0.27)
my_vector <- qbinom(0.5, 5000, list.of.p.parameter)

print(my_vector)

此方法代码更简洁，计算速度更快，尤其适合长参数列表。

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