实 验(实训)报 告
项 目 名 称 R语言初步 所属课程名称 抽样调查 项 目 类 型 综合性实验 实验(实训)日期 2015年 4月 15 日
班 级 学 号 姓 名 指导教师 罗 季
浙江财经大学教务处制
1
实验名称: 简单随机抽样(一)
实验目的: 进一步了解R语言应用;从实践中证明: 实验要求: 编写R程序,实现: (1) 从R中装载数据集”cars”;
(2) 当Iteration=100,N=50,n=10时,的误差; (3) 当Iteration=500,N=50,n=10时,的误差; (4) 当Iteration=1000,N=50,n=10时,的误差;
(5) 编写循环语句,计算Iteration 从100变化到2000,步长为100时,的误差;
(6) 编写程序,画出(5)中,误差的变化情况; (7) 将(5)中结果记录至”record”文件.
作业: 将上述任务逐个完成,并加以注释。
.
2
实验过程:
首先更改R的工作目录。本次更改在d:/r
一、 加载“cars”数据集。首先,在R软件里输入程序data(\,然后
输入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
3
结论:由R软件最后输出的误差:0.227
三、 当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) }
4
ybarmean<-mean(ybar) ybarmean
error<-ybarmean-Ybar error
结论:则当Iteration=500,N=50,n=10时,的误差:0.1908
四、 当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)
5