抽样技术第一次作业.docxVIP

抽样技术作业一学号：2013014850 姓名：金飞进 专业： 统计学2.5研究某小区家庭用于文化方面（报刊、电视、网络、书籍等）的支出，N=200，现抽取一个容量为20的样本，调查结果列于下表：编号文化支出1200215031704150516061307140810091101014011150121601318014130151501610017180181001917020120估计该小区平均的文化支出，并给出置信区间水平95%的置信区间。解：N：抽样总体 n：样本容量导入数据求求y的标准差、方差求的方差求的置信区间运行结果：得到的置信区间为[132.55,156.45]代码如下：N - 200n -20f - n/Nalpha - 0.05mydata - read.csv(C:\\Users\\Administrator\\Desktop\\ex2.5.csv)y - mydata$文化支出(ybar - mean(y))sy - sd(y)var_y - sy^2 var_ybar - (1-f)/n*sy^2radius - qnorm(1-alpha/2)* sqrt(var_ybar)ybar - radiusybar + radius2.9如果在解决习题2.5的问题时可以得到这些家庭的月总支出，见下表：编号文化支出总支出120023002150170031702000415015005160170061301400714015008100120091101200101401500111501600121601700131802000141301400151501600161001200171801900181001100191701800201201300而全部家庭的总支出平均为1600元，利用比估计的方法估计平均文化支出，给出置信水平95%的置信区间，并比较比估计和简单估计的效率。解：指标照旧由于有辅助变量x，所以还要求相应的均值，由公式求得求变量y、x的标准差、方差相关系数 套用书本所给公式，求取置信区间求两者之比作为效率运行结果：可以得到置信区间为[143.59,149.07],效率约为5.27%代码如下：N - 200n -20f - n/Nalpha - 0.05mydata - read.csv(C:\\Users\\Administrator\\Desktop\\ex2.9.csv)y - mydata$文化支出(ybar - mean(y))x - mydata$总支出(xbar - mean(x))R - ybar/xbarX - 1600yr - R*Xsy - sd(y)sx - sd(x)var_y - sy^2var_x - sx^2rho - cov(y,x)/(sy*sx) var_yr -(1-f)/n*(var_y-2*R*rho*sy*sx+R^2*var_x)radius - qnorm(1-alpha/2)* sqrt(var_yr)yr - radiusyr + radiusvar_ybar - (1-f)/n*var_yvar_yr/var_ybar2.10某养牛场购进了120头肉牛，购进时平均体重100千克。现从中抽取10头，记录重量，3个月后再次测量，结果如下表所示：编号原重量现重量1951502971553871404120180511017561151857103165810216099215010105170用回归估计法计算120头牛现在的平均重量，计算其方差的估计，并和简单估计的结果进行比较。解：和前面的题目一样，赋值指标数值，导入数据求两个变量的均值求，以及两个变量的标准差、方差相关系数简单估计（这里求的是现在的重量，所以用x）回归估计，先求，（由于这里求的是现在的重量，所以公式要变形下）运行结果：可知：简单估计的方差为19.45，而回归估计的方差为1.097，明显小于简单估计的，所以在这道题目里，用回归估计更好些代码如下：N - 120n -10f - n/Nalpha - .05mydata - read.csv(C:\\Users\\Administrator\\Desktop\\ex2.10.csv)y - mydata$原重量(ybar - mean(y))x - mydata$现重量(xbar - mean(x))R - ybar/xbarX - 100yr - R*Xsy - sd(y)sx - sd(x)var_y - sy^2var_x - sx^2rho - cov(y,x)/(sy*sx)var_ybar - (1-f)/n*var_xbeta -rho*sx/syylr - xbar+beta*(X-ybar)var_yl