抽样调查实验二.docVIP

实 验（实训）报 告 项 目 名 称 R语言初步 所属课程名称 抽样调查 项 目 类 型 综合性实验 实验(实训)日期 201年 月 日 班 级 学 号 姓 名 指导教师 罗 季 浙江财经教务处制 实验名称: 简单随机抽样（一） 实验目的: 进一步了解R语言应用；从实践中证明: . 实验要求: 编写R程序,实现： 从R中装载数据集”cars”; 当Iteration=100,N=50,n=10时,的误差; 当Iteration=500,N=50,n=10时,的误差; 当Iteration=1000,N=50,n=10时,的误差; 编写循环语句,计算Iteration 从100变化到2000,步长为100时,的误差; 编写程序,画出(5)中,误差的变化情况; 将(5)中结果记录至”record”文件. 作业: 将上述任务逐个完成,并加以注释。 实验过程: 加载“cars”数据集。首先，在R软件里输入程序data(cars)，然后输入cars获得如下数据speed dist 当Iteration=100,N=50,n=10时,的误差Y-cars[,2] N-length(Y) N Ybar-mean(Y) Ybar Ite-100 n-10 ybar-rep(0,Ite) for(i in 1:Ite){ y-sample(Y,n) ybar[i]-mean(y) } ybarmean-mean(ybar) ybarmean error-ybarmean-Ybar error 结论：由R软件最后输出的误差当Iteration=500,N=50,n=10时,的误差;Y-cars[,2] N-length(Y) N Ybar-mean(Y) Ybar Ite-500 n-10 ybar-rep(0,Ite) for(i in 1:Ite){ y-sample(Y,n) ybar[i]-mean(y) } ybarmean-mean(ybar) ybarmean error-ybarmean-Ybar error 结论：则当Iteration=500,N=50,n=10时,的误差当Iteration=1000,N=50,n=10时,的误差Y-cars[,2] N-length(Y) Ybar-mean(Y) Ybar Ite-1000 n-10 ybar-rep(0,Ite) for(i in 1:Ite){ y-sample(Y,n) ybar[i]-mean(y) } ybarmean-mean(ybar) ybarmean error-ybarmean-Ybar error 结论：当Iteration=1000,N=50,n=10时,的误差编写循环语句,计算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) 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 0 0 0-0 0 [7] 0 0 0 0-0-0[13] 0-0 0-0 0-0[19] -0 0六、编写程序,画出(5)中,误差的变化情况;plot(1:times,error,ylim=c(-1,1)) abline(h=0,col=red) 获得如下图形 将(5)中结果记录至”record”文件record-error record