抽样调查2
- 格式:doc
- 大小:1.25 MB
- 文档页数:8
实 验(实训)报 告 项 目 名 称 简单随机抽样(三),分层抽样 所属课程名称 抽样调查 项 目 类 型 综合性实验 实验(实训)日期 2015年 5月 13日
班 级 14应用统计C 学 号 姓 名 指导教师 罗 季
浙江财经大学教务处制 实验名称: 简单随机抽样(三),分层抽样 实验目的: 编写R程序解决简单随机抽样及分层抽样中的统计问题 实验要求: 解题,编写R程序,实现. 1. 调查某个地区的养牛头数,以村作为抽样单元。根据村的海拔高度和人口密度划分成四层,每层取10个村作为样本单元,经过调查获得下列数据
层 村总数 样本村养牛头数 1 2 3 4 5 6 7 8 9 10 1 2 3 4 1411 4705 2558 14997 43 84 98 0 10 44 0 124 13 0 50 147 62 87 84 158 170 104 56 160 228 262 110 232 139 178 334 0 63 220 17 34 25 34 36 0 25 7 15 31
要求: (1) 估计该地区养牛总头数Y~及其估计量的变异系数YYsˆ)ˆ(。 (2) 求分层随机抽样的设计效应。 (3) 若样本量不变采用Neyman分配可以减少方差多少?
某乡欲估计今年的小麦总产量,全县共实验过程: 一,(1) N1<-1411 N2<-4705 N3<-2558 N4<-14997 y1<-c(43,84,98,0,10,44,0,124,13,0) y2<-c(50,147,62,87,84,158,170,104,56,160) y3<-c(228,262,110,232,139,178,334,0,63,220) y4<-c(17,34,25,34,36,0,25,7,15,31) N<-N1+N2+N3+N4 W1<-N1/N W2<-N2/N W3<-N3/N W4<-N4/N ybar1<-mean(y1) ybar2<-mean(y2) ybar3<-mean(y3) ybar4<-mean(y4) ybar<-W1*ybar1+W2*ybar2+W3*ybar3+W4*ybar4 Y<-N*ybar Y [1] 1353572 VY<-(N1*(N1-n)*var(y1))/n+(N2*(N2-n)*var(y2))/n+(N3*(N3-n)*var(y3))/n+(N4*(N4-n)*var(y4))/n SY<-sqrt(VY)/Y SY [1] 0.09098019估计该地区养牛总头数1353572
及其估计量的变异系数YYsˆ)ˆ(0.09098019 (2) N1<-1411 N2<-4705 N3<-2558 N4<-14997 y1<-c(43,84,98,0,10,44,0,124,13,0) y2<-c(50,147,62,87,84,158,170,104,56,160) y3<-c(228,262,110,232,139,178,334,0,63,220) y4<-c(17,34,25,34,36,0,25,7,15,31) N<-N1+N2+N3+N4 W1<-N1/N W2<-N2/N W3<-N3/N W4<-N4/N ybar1<-mean(y1) ybar2<-mean(y2) ybar3<-mean(y3) ybar4<-mean(y4) ybar<-W1*ybar1+W2*ybar2+W3*ybar3+W4*ybar4 vybar<-(W1^2*var(y1)+W2^2*var(y2)+W3^2*var(y3)+W4^2*var(y4))/10-(W1*var(y1)+W2*var(y2)+W3*var(y3)+W4*var(y4))/N
vsrs=(N-10)/(10*(N-1))*(W1^2*var(y1)+W2^2*var(y2)+W3^2*var(y3)+W4^2*var(y4)-(W1^2*var(y1)+W2^2*var(y2)+W3^2*var(y3)+W4^2*var(y4))/10+W1*ybar^2+W2*ybar^2+W3*ybar^2+W4*ybar^2-ybar^2+vyabr) deff<-vybar/vsrs deff [1] 1.011586 (3) y1<-c(43,84,98,0,10,44,0,124,13,0) y2<-c(50,147,62,87,84,158,170,104,56,160) y3<-c(228,262,110,232,139,178,334,0,63,220) y4<-c(17,34,25,34,36,0,25,7,15,31) N1<-1411 N2<-4705 N3<-2558 N4<-14997 N<-N1+N2+N3+N4 n<-40 n1<-n*N1*sd(y1)/(N1*sd(y1)+N2*sd(y2)+N3*sd(y3)+N4*sd(y4)) n2<-n*N2*sd(y2)/(N1*sd(y1)+N2*sd(y2)+N3*sd(y3)+N4*sd(y4)) n3<-n*N3*sd(y3)/(N1*sd(y1)+N2*sd(y2)+N3*sd(y3)+N4*sd(y4)) n4<-n*N4*sd(y4)/(N1*sd(y1)+N2*sd(y2)+N3*sd(y3)+N4*sd(y4)) w1<-N1/N w2<-N2/N w3<-N3/N w4<-N4/N ybar1<-mean(y1) ybar2<-mean(y2) ybar3<-mean(y3) ybar4<-mean(y4) ybar<-w1*ybar1+w2*ybar2+w3*ybar3+w4*ybar4 Y<-N*ybar v<-(w1^2*var(y1)/4+w2^2*var(y2)/12+w3^2*var(y3)/14+w4^2*var(y4)/10)-(1/N)*(w1*var(y1)+w2*var(y2)+w3*var(y3)+w4*var(y4)) v [1] 23.40882
第一小题的方差是vybar<-(W1^2*var(y1)+W2^2*var(y2)+W3^2*var(y3)+W4^2*var(y4))/10-(W1*var(y1)+W2*var(y2)+W3*var(y3)+W4*var(y4))/N Vbar [1] 27.06595
2,方差可以减少:27.06595-23.40882=3.65713
第二题
2. 有123个村,按简单随机抽样抽取13个村作为样本,取得资料如下: 样本村 去年的小麦产量(百斤) 今年的小麦产量(百斤) 1 2 3 4 5 6 7 8 9 10 11 12 13 550 720 1500 1020 620 980 928 1200 1350 1750 670 729 1530 610 780 1600 1030 600 1050 977 1440 1570 2210 980 865 1710 (1) 若已知去年的小麦总产量为128200(百斤),采用比估计法估计今年的小麦总产量和置信度为95%的置信区间; (2) 采用回归估计法估计今年的小麦总产量和置信度为95%的置信区间; 比较两种估计方法的精确度。 Q<-function(Yr,sd,alpha){ u<-qnorm(1-alpha/2) q1<-Yr-u*sd q2<-Yr+u*sd Q<-c(q1,q2) } x<-c(550,720,1500,1020,620,980,928,1200,1350,1750,670,729,1530) y<-c(610,780,1600,1030,600,1050,977,1440,1570,2210,980,865,1710) xbar<-mean(x) ybar<-mean(y) N<-123 n<-13 X<-128200 varx<-var(x) vary<-var(y) R<-ybar/xbar Yrbar<-R*X sxy<-cov(x,y) f<-n/N varyrbar<-(1-f)*(varx+R^2*vary-2*R*sxy)/n varYrbar<-N^2*varyrbar sd<-sqrt(varYrbar) alpha<-0.05 u<-qnorm(1-alpha/2) C<-Q (Yrbar,sd,alpha) [1] 133950.5 157937.0
总量估计量为145943.8,置信区间为[133950.5 157937.0]
(2) Q<-function(Yr,sd,alpha){ u<-qnorm(1-alpha/2) q1<-Yr-u*sd q2<-Yr+u*sd Q<-c(q1,q2) x<-c(550,720,1500,1020,620,980,928,1200,1350,1750,670,729,1530) y<-c(610,780,1600,1030,600,1050,977,1440,1570,2210,980,865,1710)
N<-123 n<-13 f<-n/N X<-128200 xbar<-mean(x) ybar<-mean(y) Xbar<-X/N varx<-var(x) vary<-var(y) sxy<-cov(x,y)