}
ybarmean<-mean(ybar) ybarmean
error<-ybarmean-Ybar error
结论：当Iteration=1000,N=50,n=10时,的误差为：0.1576.
五、 编写循环语句,计算Iteration 从100变化到2000,步长为100时,的误差;
解：程序如下：
Y<-cars[,2] N<-length(Y) N
Ybar<-mean(Y) Ybar
Ite<-seq(from=100,to=2000,by=100) n<-10
times<-length(Ite)
6
ybarmean<-rep(0,times) for(j in Ite){ ybar<-rep(0,j) for(i in 1:j){ y<-sample(Y,n) ybar[i]<-mean(y) }
ybarmean[j/100]<-mean(ybar) }
error<-ybarmean-Ybar error
[1] -0.04400000 0.06950000 0.45166667 0.23850000 -0.09920000 0.20400000
[7] 0.23457143 0.21962500 0.06422222 0.45660000 -0.09890909 -0.02566667
[13] 0.16200000 -0.07285714 0.16560000 -0.02325000 0.09788235 -0.03666667
[19] -0.03973684 0.05570000
六、编写程序,画出(5)中,误差的变化情况;
7
解：程序如下：
输入程序：
plot(1:times,error,ylim=c(-1,1)) abline(h=0,col=\获得如下图形
七、 将(5)中结果记录至”record”文件
record<-error record
write.table(record,\
8
9