资源描述:
《抽样调查实验二》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、浙述财鮭jL厝实验(实训)报告项目名称R语言初步所属课程名称抽样调查项目类型综合性实验实验(实训)日期2015年4月15日班级学号姓名指导教师罗季浙江财经大学教务处制实验名称:简单随机抽样(一)实验目的:进一步了解R语言应用;从实践中证明:E(y)=Y.实验要求:编写R程序,实现:(1)从R中装载数据集"cars";(2)当Iteration=100,N=50,n=10时,y的误差;(3)当Iteration=500,N=50,n=10时,y的误差;(4)当Iteration^1000,N=50,n=10时,y的误差;⑸编写
2、循环语句,计算Iteration从100变化到2000,步长为100时,y的误差;(6)编写程序,画出(5)中,误差的变化情况;(7)将⑸中结果记录至Second文件.作业:将上述任务逐个完成,并加以注释。实验过程:首先更改R的工作目录。本次更改在d:/r一、加载“cars”数据集。首先,在R软件里输入程序data("cars"),然后输入cars获得如下数据speeddist二、当Iteration=100,N=50,n=10时,y的误差解:程序如下Y<-cars[,2]N<-Iength(Y)NYbar<-mean(Y)
3、YbarIte<-100n<-10ybar<-rep(0Jte)for(iinl:Ite){y<-sample(Y,n)ybar[i]<-mean(y)}ybarmean<・mean(ybar)ybarmeanerrorv・ybarmcan・Ybarerror疳RConsole>Y<-cars[,2]R未命名-R编辑器>N<-length(Y)Y<-cars[,2]>NN<-length(Y)[1]50N>Ybar<-inean(Y)Ybar<-mean(Y)>YbarYbar[1]42・98Ite<-100>Ite<-10
4、0n<-10>n<-10ybar<-rep(0,Ite)>ybar<-rep(0^Ite)for(iin1:Ite){>for(iin1:Ite){y<-sample(Y#n)+y<-sample(Y,n)ybar[i]<-mean(y)+ybar[i]<-mean(y)}+}ybarmean<-mean(ybar)>ybarmsan<-mean(ybar)ybarmean>ybarmsanerror<-ybarmean-Ybar[1]43・207error>error<-ybarinean-Yb>error[1]0.227
5、>1结论:由R软件最后输出y的误差0227三、当Iteration=500,N=50,n=10时,y的误差;解:程序如下Y<-cars[,2]N<-length(Y)NYbar<-mean(Y)YbarIte<-500n<-10ybar<-rep(OJte)for(iinl:Ite){y<-sample(Y,n)ybar[i]<-mean(y)ybarmean<-mean(ybar)ybarmeanerror<-ybarmean-Ybarerror
6、RRConsole>Y<-cars[,2]>N<-length(Y)>>Yb
7、ar<-mean(Y)>Ybar[1]42・98>Ite<-500>n<-10>ybar<-rep(0,Ite)>for(iin1:Ite){+y<-sample(Yzn)+ybar[i]<-mean(y)+}>ybarmsan<-mean(ybar)>ybarmean[1]43・1708>error<-ybarmean-Ybar>error[1]0.1908结论:则当Iteration=500,N=50,n=10时,y的误差:0.1908四、当Iteration=1000,N=50,n=10时,y的误差解:程序如下:Y<-
8、cars[,2]N<-length(Y)Ybar<-mean(Y)YbarIte<-1000n<-10ybar<-rep(O,Ite)for(iinl:Ite){y<-sample(Y,n)ybar[i]<-mean(y)ybarmeanv・mean(ybar)ybarmcanerrorv・ybarmean・YbarerroronsoJ>Y<-cars[z2]>N<-length(Y)>>Ybar<-raean(Y)>Ybar[1]>Ite<-1000>n<-10>ybar<-rep(0,Ite)>for(iin1:Ite)
9、{+y<-sample(Y,n)+ybar[i]<-mean(y)+}>ybarmsan<-mean(ybar)>ybarmean[1]43・1376>error<-ybarinean-Ykiar>error[1]0.1576>I结论:当Iteration=1000,N=50,n=10时,