我想在我构建的矩阵的每一行上运行t.test。然后我尝试拉出置信区间的值并将它们保存在单独的向量中,每次迭代一个输出。但是,在我运行代码后,我得到L1 = NA NA NA .... 8.155677。如果你能指出错误,我会很高兴。
(我知道有很多方法可以编写这个代码清理程序但是,我试着逐步编写它。)
set.seed(1234)
n= 24 # sample size or a number of RV's
N=100 # number of exrtractions or a number of sums for each rv
X=rnorm(N*n, 9, 1.5 ) # generate rv's
XMat=matrix(X,nrow=N)
#Problem部分:
L1=c()
L2=c()
for(i in N)
{
s=XMat[i,1:n]
K=t.test(s,conf.level=0.95)
M=K$conf.int
l1=M[1]
l2=M[2]
L1[i]=l1
L2[i]=l2
}
答案 0 :(得分:0)
将循环控制更改为:
for (i in seq(N))
您正在代码中运行单个i
值的循环。