数学实验——统不计推断.docxVIP

实验9 统计推断分1 黄浩 2011011743实验目的掌握数据的参数估计、假设检验的基本原理、算法，及用MATLAB 实现的方法练习用这些方法解决实际问题实验内容1.《数学实验》（第一版）习题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 1182月 118 119 115 122 118 121 120 122 128 116 120 123 121 119 117 119 128 126 118 1251）分别用两个月的数据验证这种说法的可靠性；2）分别给出1月和2月汽油价格的置信区间；3）如何给出1月和2月汽油价格差的置信区间。实验过程：（1）（2）小问：第一问验证可靠性，既可以用1、2月的均值来比较，也可以用假设检验的方法来进行推断，当然后者对于问题的分析更有信度，因为假设检验不仅利用了均值，还利用了样本方差这一统计量。因此，可以将1、2小问合并，使用代码：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];[mu1 sigma1 muci1 sigmaci1]=normfit(x1,0.05)[mu2 sigma2 muci2 sigmaci2]=normfit(x2,0.05)所得结果如下：（显著性水平为0.05）样本均值样本方差总体均值区间估计1月115.153.87[113.34,116.96]2月120.753.712[119.01,122.49]对以上表格进行分析，得：1月数据对总体均值的点估计为115.15，而2月数据对总体均值的点估计为120.75，即1月的数据可以支持题中说法，而相对而言2月的数据不能支持该种说法。当然，仅仅从点估计不能很好的说明问题，因为它没有给出到底样本均值偏离多少的时候，就不能接受该说法，也无法了解1月份到底有多大程度支持该说法、2月份有多大程度不支持该说法。因此，这就要用到假设检验，记油价总体均值为，并设：使用代码：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];[h1 sig1 ci1]=ttest(x1,115,0.05,0)[h2 sig2 ci2]=ttest(x2,115,0.05,0)结果如下：H0P值总体均值置信区间1月接受0.8642[113.34,116.96]2月拒绝1.324*10^-6[119.01,122.49]由上表可见，在显著性水平为0.05的条件下，1月份的数据支持原假设，而2月份的数据不支持原假设。而且，在原假设成立的条件下，1月份数据的P值为0.8642，很接近1，说明可信度比较高，而2月份的P值的数量级在10^-6，这说明在H0成立的条件下，2月份的数据只可能是小概率事件，因而通过反证法可知，2月份的数据是反对原假设的。此外，使用假设检验给出的置信区间与用区间估计给出的区间是一致的，这是由于二者都是用了相同的统计量和相同的分位数，实际上，在进行上一步的区间估计的时候，就已经能做出接受/拒绝的结论了。但是，在matlab的具体实现中，使用ttext函数进行假设检验还可以直接得到接受和拒绝的结论、以及相应的P值，既方便又全面，比使用normfit进行区间估计更强大。（3）小问：对于这一问，我首先想到的是将两行数据相减，然后进行区间估计，直接给出置信区间，使用代码：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];[mu sigma muci sigmaci]=normfit(x1-x2,0.05)结果如下：样本均值样本方差总体均值区间估计价格差-5.