data<-rnorm(200,mean=15,sd=2)
y<-replicate(10,sample(data,100))
b=matrix(0,nrow=100,ncol=4)b
u<-qnorm(0.95)
for(i in 1:100)
{
b[i,1]=mean(y[i,])
b[i,2]=sd(y[i,])
b[i,3]=b[i,1]-u*sqrt(b[i,2])
b[i,4]=b[i,1]+u*sqrt(b[i,2])
}
这是我的代码。