TA希望我创建一个新的随机变量Y_n=sum(X_i)，其中X_i是n个二项式随机变量，N = 4且p = 1/3。这不是太糟糕;我只使用以下for循环：for(i in 1:100){yn[i] <- c(sum(rbinom(i, 4, (1/3))))}。然而，他然后希望我从1到10,000（即10,20,30，......，9990,10000）每十分之一的数字重新创建Y_n。我尝试使用此代码：yseq <- seq(10, 10000, by=10)
for(i in yseq){
Y2[i] <- c(sum(rbinom(i,4,(1/3))))}。它有点有效，但不是真的。它返回一个列表（我检查了它的类），看似正确的值，但是一堆NA。这给我带来了两个问题：1）R不会让我将列表重新分类为向量，2）R告诉我列表长度为1，这是一堆垃圾。
有些可以告诉我哪里出错了吗？我以前说过：编程不是我的强项，但我总是尽力学习！
谢谢！