资源描述:
《用r软件进行区间估计》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数理统计上机报告上机实验题目:用R软件进行区间估计上机实验目的:为了了解并理解区间估计的概念,学会并且记住常见的区间估计的表示方法,用R软件实现区间估计。区间估计基本理论、方法:首先是利用假设检验判断接受还是拒绝原假设,在原假设为真的情况下利用样本数据对总体的区间估计,选用统计函数,利用样本数据()来对总体进行估计,再利用已知的置信度求置信水平α,再判断总体参数的区间。实验实例和数据资料:P377习题7.16随机的从一批钉子中抽取16枚,测得其长度为(单位:cm):2.142.102.132.152.132.122.132.102.152.122.142.
2、102.132.112.142.11设钉长的分布为正态分布,分别对下列两种情况求出总体均值μ的90%置信度的置信区间(i)已知σ=0.01cm;(ii)σ未知上机实验步骤:(i)x<-c(2.14,2.10,2.13,2.15,2.13,2.12,2.13,2.10,2.15,2.12,2.14,2.10,2.13,2.11,2.14,2.11)sigma<-0.01xbar<-mean(x)alpha<-0.10fws<-qnorm(1-alpha/2,lower.tail=TRUE)left<-Xbar-(fws*sigma)/(sqrt(length
3、(x)))right<-Xbar+(fws*sigma)/(sqrt(length(x)))实例计算结果及分析:3所以得到的置信区间为[2.121,2.129](ii)x<-c(2.14,2.10,2.13,2.15,2.13,2.12,2.13,2.10,2.15,2.12,2.14,2.10,2.13,2.11,2.14,2.11)xbar<-mean(x)S1<-var(x)n<-length(x)alpha<-0.10fws<-qt(1-alpha/2,Chang-1,lower.tail=TRUE)left<-xbar-(fws*sqrt(S1)
4、)/(sqrt(n))right<-xbar+(fws*sqrt(S1))/(sqrt(n))实例计算结果及分析:所以得到的置信区间为[2.1175,2.1325]P378习题7.20为了在正常条件下检验一种杂交作物的两种新处理方案,在同一地区随机的挑选8块地,在每一块实验地上按两种方案种植作物,这8块地的单位面积产量分别是:1号方案产量8687569384937579807958917782746632好方案产量设这两种方案的产量都服从正态分布,失球这两个平均产量之差的置信度为95%的一个置信区间。上机实验步骤:x<-c(86,87,56,93,84,9
5、3,75,79)y<-c(80,79,58,91,77,82,74,66)xbar<-mean(x)ybar<-mean(y)n1<-length(x)n2<-length(y)s1<-var(x)s2<-var(y)alpha<-0.05sw<-sqrt((n1-1)*s1+(n2-1)*s2)k<-sqrt(n1*n2*(n1+n2-2))fws<-qt(1-alpha/2,n1+n2-2,lower.tail=TRUE)left<--fws*sw*(sqrt(n1+n2))/k+(xbar-ybar)right<-fws*sw*(sqrt(n1+n2
6、))/k+(xbar-ybar)实例计算结果及分析:实验得到的置信区间为[-6.187,17.687]3