數(shù)學實驗-11：統(tǒng)計推斷

1、精選優(yōu)質(zhì)文檔-傾情為你奉上實驗 11：統(tǒng)計推斷習題2：據(jù)說某地汽油的價格是115美分/gal，為了驗證這種說法，一位司機開車隨機擇了一些加油站，得到某年1月和2月的數(shù)據(jù)如下（不完整）：1月1191171151161091091182月1181191151221261181251)分別用兩個月的數(shù)據(jù)檢驗這種說法的可靠性；2)分別給出1月和2月汽油價格的置信區(qū)間（=0.05）3)如何給出1月和2月汽油價格差的置信區(qū)間（=0.05）1 模型建立本題中應該假設兩個月該地區(qū)加油站的油價整體符合正態(tài)分布，可以進行正態(tài)性檢驗。隨后進行假設檢驗：為價格平均值，由此可以得到1月和2月汽油價格的置信區(qū)間；對于兩個

2、月的價格差的區(qū)間估計問題，應該分兩種情況進行討論，一種是認為兩個月的數(shù)據(jù)對應行是來自同一個加油站，則首先應求得兩個月差價的20組數(shù)據(jù)，然后在顯著性水平下確定的置信區(qū)間：如果認為兩個月的數(shù)據(jù)是完全隨機的，就不能采用1、2月對應數(shù)據(jù)的差作為新的數(shù)據(jù)，應采用兩總體的t分布檢驗，可以利用假設檢驗的命令實現(xiàn)區(qū)間估計。2 程序設計問題1和問題2：1）對1月和2月的數(shù)據(jù)分別進行正態(tài)性檢驗：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 1

3、20 122 128 116 120 123 121 119 117 119 128 126 118 125;mu=115;%原始數(shù)據(jù)h1,p1,lstat1,cv1=lillietest(x1)h2,p2,lstat2,cv2=lillietest(x2)%進行Lilliefors檢驗，適用于一般的正態(tài)性檢驗，即Ho:總體服從正態(tài)分布，其中，為樣本均值和方差估計。2）對1和2月的數(shù)據(jù)分別進行假設檢驗，并得到置信區(qū)間，由于總體方差未知，所以采用t檢驗法：mean(x1),mean(x2)%輸出均值h1,sig1,ci1,z1=ttest(x1,mu,0.05) h2,sig2,ci2,z2=t

4、test(x2,mu,0.05)3)針對數(shù)據(jù)的兩種可能采用不同的方法對兩個月的數(shù)據(jù)進行區(qū)估計，仍采用假設檢驗的命令：%假設兩個月的數(shù)據(jù)是對應的：mean(x1-x2)h3,sig3,ci3,z3=ttest(x2-x1,0,0.05);%假設兩個月的數(shù)據(jù)是完全隨機的：h4,sig4,ci4=ttest2(x1,x2,0.05);%兩個總體的均值的假設檢驗3 運行結(jié)果和分析1）對1月和2月的數(shù)據(jù)分別進行正態(tài)性檢驗：h1 = 0p1 = 0.0000lstat1= 0.1105cv1 = 0.0000h2 = 0p2 = 0.1399lstat2= 0.4013cv2 = 0.0000接受1、2月

5、數(shù)據(jù)負總正態(tài)分布的假設，可以進行假設檢驗。2）對1月和2月的數(shù)據(jù)分別進行假設檢驗 1.0000e+002 % mean(x1) 1.0000e+002 % mean(x2)h1 = 0%接受1月汽油價格為115美分/gal的假設sig1 = 0.9275% ci1 = 1.0e+002 * 1.9442 1.0558 %1月汽油價格的置信區(qū)間，可以看到115在區(qū)間內(nèi)。z1 = tstat: 0.3442%統(tǒng)計量 df: 19 sd: 3.7756h2 = 1 %拒絕接受2月汽油價格為115美分/gal的假設sig2 = 1.8971e-006ci2 = 1.0e+002 * 1.5904 1.

6、4096%2月汽油價格的置信區(qū)間，可以看到115不在區(qū)間內(nèi)z2 = tstat: 6.7383 df: 19 sd: 3.0572 結(jié)果：接受1月份價格為115美分/gal的假設，其置信區(qū)間為113.3388，116.9612;不接受2月份價格為115美分/gal的假設，其置信區(qū)間為119.0129，122.48713)針對數(shù)據(jù)的兩種可能采用不同的方法對兩個月的數(shù)據(jù)進行區(qū)估計，仍采用假設檢驗的命令：u 假設兩組數(shù)據(jù)來自同一個加油站 x1-x2=1 -2 0 -6 -6 0 -5 0 -12 2 -11 -11 -2 -7 0 -6 -14 -17 -9 -7;mean(x1-x2)=-5.6;

7、h3 = 1%拒絕x1-x2=0（、月份價格不變）的假設，實際上這一問是做區(qū)間估計，假設檢驗的結(jié)果無所謂sig3 = 2.5117e-004ci3 = 3.2156 8.7844%x1-x2的置信區(qū)間z3 = tstat: 4.4384 df: 19 sd: 5.3948結(jié)果：，即、月份汽油價格差在.03938.1607區(qū)間內(nèi)的概率是。u 假設兩組數(shù)據(jù)完全獨立h4 = 1拒絕、月份價格不變的假設sig4 = 3.9776e-005ci4 = -8.3219 -3.6770置信區(qū)間結(jié)果：，即、月份汽油價格差在3.17278.0274區(qū)間內(nèi)的概率是。兩種計算得到的結(jié)果略有不同。4 SIG參數(shù)的意義

8、進行假設檢驗時得到的參數(shù)，其意義為，可以自編程序?qū)崿F(xiàn)對的重現(xiàn)：t=(mean(x2)-mu)/(std(x2)/sqrt(20)%統(tǒng)計量（月）tt=tcdf(t,19);%，ttt=tpdf(t,19);%值對應的概率密度函數(shù)值 sig22=2*(1-tt)%x=-5:0.01:5;plot(x,tpdf(x,19),t,ttt,'r*'); %繪制（）分布的圖形，并標記出統(tǒng)計%量的位置結(jié)果：t = 6.7383sig22 = 1.8579e-006同之前直接利用ttest命令得到的結(jié)果相同（藍色標注）?？梢钥吹浇y(tǒng)計量的位置已經(jīng)大大超過了的范圍，sig 的值是上圖中紅點和其關于

9、點對稱的綠點外側(cè)對應的概率，根據(jù)其求解過程，可以將其定義為：在一組統(tǒng)計數(shù)據(jù)下，錯誤的否定正確的Ho假設的概率。這是建立在假設Ho恒成立的基本前提下的一種判別方法，事實上體現(xiàn)了假設檢驗的基本思想。對1月數(shù)據(jù)進行同樣的處理：t = 6.7383sig11 = 0.9275習題8：20名學生參加了某課程進行的、考察同樣知識的兩次測驗，成績?nèi)缦卤?，根?jù)這些數(shù)據(jù)判斷兩次檢驗的難度是否相同。表：20名學生的兩次測驗成績（每列是同一名學生的兩次成績）第一次938579 90838385第二次888986858888901. 程序設計1）本題的兩個樣本不是獨立的。首先應該檢驗兩次測驗成績的差是否符合總體正態(tài)分

10、布，如果符合，可以對差進行假設檢驗的方法進行處理。x=93 85 79 90 78 76 81 85 88 68 92 73 88 84 90 70 69 83 83 85;. 88 89 86 85 87 88 75 93 88 78 86 86 80 89 85 79 78 88 88 90;h0,p0=lillietest(x(1,:)- x(2,:)2）在正態(tài)檢驗合格的前提下，對兩組數(shù)據(jù)的差，進行假設檢驗,其中，Ho: x1-x2=0:mean(x(1,:)-x(2,:)mu=0;alpha=0.05;h1,sig1,ci1,z1=ttest(x(1,:)-x(2,:),mu,alpha)2.運行結(jié)果和分析1) 正態(tài)性檢驗的結(jié)果：p0 = 0.57060.050.0350.0360.037h01011根據(jù)參數(shù)的意義，p0表示：x(1,:)-x(2,:),即成績之差符合正態(tài)分布但是被錯誤決絕的概率；作為顯著性水平，其越小允許錯誤篩去符合正態(tài)的數(shù)據(jù)的概率越小。對于默認的=0.05,不能接受數(shù)據(jù)為正態(tài)的假設，但是只要略微放寬約束，=0.35，即如果成績之差的確為正態(tài)，僅有3.5%的概率被篩去，此時可以接受成績之差符合正態(tài)分布的假設。由于值仍然不是太小，所以可以接受并繼續(xù)運算。2) 假設檢驗的結(jié)果：h1 = 1%拒絕接受兩次考試成績相同(難度相同)