【发布时间】:2016-12-11 23:08:47
【问题描述】:
我想在我构建的矩阵的每一行上运行一个 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)
#问题部分:
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
}
【问题讨论】: