实验12 统计推断资料.docxVIP

实验12 统计推断 【实验目的】 掌握数据的统计描述和参数估计、假设检验的基本概念与原理，及用MATLAB实现的方法； 练习用这些方法解决实际问题。 【实验内容】 【题目2】 据说某地汽油的价格是115美分/gal，为了验证这种说法，一位司机开车随机选择了一些加油站，得到某年1月和2月的数据如下： 1月 119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 114 109 109 118 2月 118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 125 （1）分别用两个月的数据验证这种说法的可靠性； （2）分别给出1月和2月汽油价格的置信区间（α=0.05） （3）如何给出1月和2月汽油价格差的置信区间（α=0.05） 2.1 汽油价格的置信区间 2.1.1 模型分析 本题是个假设检验问题。首先设该地的汽油价格服从正态分布。其原假设与备择假设分别为，。 此时总体方差未知，则~，取得分位数，记， ，从而可以据此选择接受或者拒绝H0。 2.1.2 matlab求解 用matlab运行如下程序 x1=[119 117 115 116 112 121 115 122 116 118 109 112 119 112 117 113 114 109 109 118]; x2=[118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 125]; m1=mean(x1) m2=mean(x2) [h1,sig1,ci1]=ttest(x1,115,0.05,0) [h2,sig2,ci2]=ttest(x2,115,0.05,0) 得到如下输出结果： m1 = 115.1500 m2 = 120.7500 h1 = 0 sig1 = 0.8642 ci1 = 113.3388 116.9612 h2 = 1 sig2 = 1.3241e-006 ci2 = 119.0129 122.4871 2.1.3 结果分析 月份均值检验结果价格的置信区间一月115.1500接受[113.3388 116.9612]二月120.7500拒绝[119.0129 122.4871]从一月份的数据来看，应接受平均有家为115美分/gal的说法，从二月份的数据看则不可以接受。 一月份汽油价格的置信区间为[113.3388 116.9612]，二月份汽油价格的置信区间为[119.0129 122.4871]。 2.2 汽油价格差的置信区间 对于此题有两种理解，一是两月的数据分别一一对应，来自于同一加油站；二是两组数据的来源均是随机的，来自于不同加油站。 2.2.1 两月数据一一对应来自同一加油站 对于此种理解，则两个月的差值及时一组样本。故课分别求出每个加油站的两个月份油价差，产生一组新的随机变量，直接对这组数据进行t检验，得到其差值的置信区间。 x3=x2-x1; [h3,sig3,ci3]=ttest(x3,0,0.05,0) 输出结果如下： h3 = 1 sig3 = 2.0582e-004 ci3 = 3.0393 8.1607 因此，置信区间为[3.0393 8.1607]。 2.2.2 两组数据来自不同加油站 在此种理解下，应构造T统计量由t分布检验，把一月与二月的数据看作两个总体，其均值的差值应该服从两个总体时的t 分布。取~，其中。。当t 满足处于置信区间内时我们可以认为一月和二月该地的汽油价格是相同的。该题要求确定置信区间，即为假定一月和二月的汽油价格不变，差值必定满足显著性水平，所以差值区间即为置信区间。 [h4,sig4,ci4]=ttest2(x2,x1,0.05,0) 得到如下结果： h4 = 1 sig4 = 3.6952e-005 ci4 = 3.1727 8.0273 此时置信区间为[3.1727 8.0273]。 2.3 小结 本题主要考察了 t 检验。虽然在编程时用到了ttest 的两种不同的检验，但是它们本质只是t 分布的几种变换形式。本题同时让我们复习了假设检验中的基本流程：先根据原题中的目标建立原假设，再写出备择假设，然后根据拒绝域的形式，利用正态分布、t 分布等构造出拒绝域形式，最终代入样本数据，对总体进行检验。MATLAB 利用此思想方便地利用函数进