抽样调查实验二

实 验（实训）报 告

项 目 名 称 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

}

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

本文来源：https://www.bwwdw.com/article/i8h3.html