多元統(tǒng)計(jì)分析第三章假設(shè)檢驗(yàn)與方差分析

1、第3章多元正態(tài)總體的假設(shè)檢驗(yàn)與方差分析從本章開始，我們開始轉(zhuǎn)入多元統(tǒng)計(jì)方法和統(tǒng)計(jì)模型的學(xué)習(xí)。統(tǒng)計(jì)學(xué)分析處理的對(duì)象是帶 有隨機(jī)性的數(shù)據(jù)。按照隨機(jī)排列、重復(fù)、局部控制、正交等原則設(shè)計(jì)一個(gè)試驗(yàn)，通過試驗(yàn)結(jié)果 形成樣本信息（通常以數(shù)據(jù)的形式），再根據(jù)樣本進(jìn)行統(tǒng)計(jì)推斷，是自然科學(xué)和工程技術(shù)領(lǐng)域 常用的一種研究方法。由于試驗(yàn)指標(biāo)常為多個(gè)數(shù)量指標(biāo)，故常設(shè)試驗(yàn)結(jié)果所形成的總體為多元 正態(tài)總體，這是本章理論方法研究的出發(fā)點(diǎn)。所謂統(tǒng)計(jì)推斷就是根據(jù)從總體中觀測(cè)到的部分?jǐn)?shù)據(jù)對(duì)總體中我們感興趣的未知部分作出推 測(cè)，這種推測(cè)必然伴有某種程度的不確定性，需要用概率來表明其可靠程度。統(tǒng)計(jì)推斷的任務(wù) 是“觀察現(xiàn)象，提取信息，建

2、立模型，作出推斷”。統(tǒng)計(jì)推斷有參數(shù)估計(jì)和假設(shè)檢驗(yàn)兩大類問題，其統(tǒng)計(jì)推斷目的不同。參數(shù)估計(jì)問題回答諸如“未知參數(shù) 二的值有多大？”之類的問題，而假設(shè)檢驗(yàn)回答諸如“未知參數(shù)二的值是“嗎？”之類的問題。本章主要討論多元正態(tài)總體的假設(shè)檢驗(yàn)方法及其實(shí)際應(yīng)用，我們將對(duì)一元正態(tài)總 體情形作一簡(jiǎn)單回顧， 然后將介紹單個(gè)總體均值的推斷，兩個(gè)總體均值的比較推斷， 多個(gè)總體均值的比較檢驗(yàn)和協(xié)方差陣的推斷等。3.1 一元正態(tài)總體情形的回顧一、假設(shè)檢驗(yàn)在假設(shè)檢驗(yàn)問題中通常有兩個(gè)統(tǒng)計(jì)假設(shè)（簡(jiǎn)稱假設(shè)），一個(gè)作為原假設(shè)（或稱零假設(shè)），另一個(gè)作為備擇假設(shè)（或稱對(duì)立假設(shè)），分別記為H。和H“。1、顯著性檢驗(yàn)為便于表述，假定考慮假

3、設(shè)檢驗(yàn)問題：設(shè)X1, X2，,Xn來自總體N（，二2）的樣本,我們要檢驗(yàn)假設(shè)(3.1 )Hi：原假設(shè)H。與備擇假設(shè)Hi應(yīng)相互排斥，兩者有且只有一個(gè)正確。備擇假設(shè)的意思是，一旦 否定原假設(shè)H。，我們就選擇已準(zhǔn)備的假設(shè)Hi。當(dāng)二2已知時(shí)，用統(tǒng)計(jì)量 zX在原假設(shè)Ho成立下，統(tǒng)計(jì)量z服從正態(tài)分布zn（0,1），通過查表，查得 N（0,1）的上分位點(diǎn)Z 一 2。對(duì)于檢驗(yàn)問題（3.1.1 ），我們制定這樣一個(gè)檢驗(yàn)規(guī)則（簡(jiǎn)稱檢驗(yàn)）：當(dāng)z z 2時(shí)，拒絕Ho ；當(dāng) Z WZa2時(shí)，接受 Ho。（3.2）我們稱z一很為臨界值，是N（0,1）的上分位點(diǎn)，不同的臨界值代表不同的檢驗(yàn)。稱拒絕原假設(shè)Ho的統(tǒng)計(jì)量z的范圍

4、為拒絕域，稱接受Ho的統(tǒng)計(jì)量z的范圍為接受域，因此給出一個(gè)檢驗(yàn)， 就是給出一個(gè)拒絕域。2、兩類錯(cuò)誤由于樣本具有隨機(jī)性，因此在根據(jù)樣本進(jìn)行判斷時(shí)，有可能犯兩種類型的錯(cuò)誤。一類錯(cuò)誤 是，原假設(shè)H。本來正確，但按檢驗(yàn)規(guī)則卻作出了拒絕H。的判斷，這類錯(cuò)誤稱為第一類錯(cuò)誤 （棄真錯(cuò)誤），其發(fā)生的概率 Pzza=稱為犯第一類錯(cuò)誤的概率； 另一類錯(cuò)誤時(shí)，原假設(shè)Ho 本來不正確，但按檢驗(yàn)規(guī)則卻作出了接收Ho的判斷，這類錯(cuò)誤稱為第二類錯(cuò)誤（存?zhèn)五e(cuò)誤），其發(fā)生的概率稱為犯第二類錯(cuò)誤的概率，記為一：。同時(shí)控制這兩類錯(cuò)誤是困難的，當(dāng)時(shí)在樣本容量n固定的條件下，要使:-和同時(shí)減小，通常是不可能的。在假設(shè)檢驗(yàn)的應(yīng)用中，由奈

5、曼（NEYMAN與皮爾遜（PEARSON提出了一個(gè)原則，即在控制犯第一類錯(cuò)誤的概率：條件下，盡量使犯第二類錯(cuò)誤的概率小,這種檢驗(yàn)問題，稱為顯著性檢驗(yàn)問題。根據(jù)這一原則，原假設(shè)受到保護(hù)，不至于被輕易拒絕，一旦檢驗(yàn)結(jié)果拒絕 了原假設(shè)，則表明拒絕的理由是充分的，如果接受了原假設(shè)，則只是表明拒絕的理由還不充分， 未必意味著原假設(shè)就是正確的。所以，在實(shí)際問題中，為了通過樣本觀測(cè)值對(duì)某一猜測(cè)取得強(qiáng) 有力的支持，通稱我們把這一猜測(cè)的否定作為原假設(shè)，而把猜測(cè)本身作為備擇假設(shè)。3、關(guān)于檢驗(yàn)的p值下面，我們?cè)俳榻B進(jìn)行檢驗(yàn)的另一種方式一一p值，我們就以（3.1.1）的檢驗(yàn)問題為例來X _ 4 a加以說明，對(duì)于樣本，

6、我們通過統(tǒng)計(jì)量，計(jì)算出zo =，是一確定值，這里的 X是CF /樣本觀測(cè)值的均值，再由統(tǒng)計(jì)量z服從正態(tài)分布zN（o ,1）,計(jì)算P zo為檢驗(yàn)的p值。由于z|等價(jià)于P = P z |zo蘭Pz，所以檢驗(yàn)規(guī)則可以表述為：當(dāng)p FWilks Lambda0.661127742.903170.0649Pillais Trace0.338872262.903170.0649Hotelling-Lawley Trace 0.512566992.903170.0649Roys Greatest Root 0.512566992.903170.0649可見P值為0.0649，所以否定原假設(shè)，即在0.10顯著

7、水平下拒絕 Ho。在實(shí)際工作中，一元檢驗(yàn)與多元檢驗(yàn)可以聯(lián)合使用，多元的檢驗(yàn)具有概括和全面的優(yōu)點(diǎn)，而一 元的檢驗(yàn)容易發(fā)現(xiàn)各指標(biāo)之間的關(guān)系和差異，兩者的結(jié)合能給統(tǒng)計(jì)人員提供更多的統(tǒng)計(jì)分析信 息。3.3兩總體均值的比較檢驗(yàn)例3.2 為了研究日美兩國(guó)在華企業(yè)對(duì)中國(guó)經(jīng)營(yíng)環(huán)境的評(píng)價(jià)是否存在差異，從兩國(guó)在華企業(yè)對(duì)中國(guó)的政治、經(jīng)濟(jì)、法律、文化等環(huán)境打分，得表3-2。試分析日美兩國(guó)在華企業(yè)對(duì)中國(guó)經(jīng)營(yíng)環(huán)境的評(píng)價(jià)是否存在差異？表3-2日美兩國(guó)在華企業(yè)對(duì)中國(guó)經(jīng)營(yíng)環(huán)境的評(píng)價(jià)美國(guó)企業(yè)號(hào)政治環(huán)境X1經(jīng)濟(jì)環(huán)境X2法律環(huán)境X3文化環(huán)境X4美165352560美275502055美360453565美475404070美5703

8、03050美655403565美760453060美865402560美960503070美1055553575日本企業(yè)號(hào)政治環(huán)境Y1經(jīng)濟(jì)環(huán)境Y2法律環(huán)境Y3文化環(huán)境Y4日155554065日250604570日345453575日450505070日555503075日660404560日765554575日850653580日940453065日1045504570假設(shè) X =（Xi，X2，X3，X4）服從 Nx匸 x），Y = （丫1,丫2,丫3,丫4）服從 N（S，） 下，且有10對(duì)樣品，要做復(fù)合檢驗(yàn)bl1131 1一般情況下，我們考慮 Xi,X2,，Xn為取自P維正態(tài)總體Np(叫，三

9、)的一個(gè)樣本, 丫1,丫2,，Ym為取自P維正態(tài)總體Np(2，D的一個(gè)樣本。假定兩組樣本相互獨(dú)立，且_ 1 n - 1 nX = Xi，丫=Yin i amy一、有共同已知的協(xié)差陣時(shí)對(duì)于例3.2提出的問題，可歸類為假設(shè)檢驗(yàn)問題：叫：亠二叮Hi：h其中叫2為已知P維向量。在一維情形下，用了統(tǒng)計(jì)量 U = X丫 n m，與前面相似的思路，在 p維時(shí)，選用統(tǒng)ct0 T n + m計(jì)量nmn mT A(X -Y) Z (X -Y)當(dāng)H。成立時(shí)，2服從2(p)-分布。對(duì)給定的顯著性水平 ，從P 2 一 2(p :，求出:(p)。當(dāng)2 一(P)時(shí)，拒絕；當(dāng) FWilks LambdaPillais Tra

10、ce0.376077346.2240.623922666.224150.0037150.0037Hotelling-Lawley Trace 1.659027526.224150.0037Roys Greatest Root 1.659027526.224150.0037由此可見p值是0.0037，因而日美兩國(guó)在華企業(yè)對(duì)中國(guó)經(jīng)營(yíng)環(huán)境的評(píng)價(jià)存在顯著差異。3.4多個(gè)總體均值向量的比較檢驗(yàn)在研究作物栽培時(shí)，要考慮播種期、品種、土質(zhì)、施肥方式、灌溉方式對(duì)產(chǎn)量的影響；在 化學(xué)反應(yīng)中要觀察原料成分、劑量、催化劑、溫度、壓力，攪拌速度等對(duì)得率的影響。在很多 應(yīng)用領(lǐng)域尤其是科學(xué)研究中，都遇到過類似的問題，常涉

11、及許多因素，這類問題要分析出影響 最“大”的因素，就是比較各種因素對(duì)試驗(yàn)結(jié)果所起的作用問題。作為影響試驗(yàn)結(jié)果的每一因 素或因素的某一水平或某一方案，且試驗(yàn)結(jié)果都形成一個(gè)隨機(jī)總體。這樣，比較各種因素對(duì)試 驗(yàn)結(jié)果所起的作用問題就變成對(duì)各種因素的試驗(yàn)結(jié)果所形成的總體的比較問題。由于試驗(yàn)指標(biāo)常為多元指標(biāo)，故常設(shè)試驗(yàn)結(jié)果所形成的總體為多元正態(tài)總體。此外，我們 按照隨機(jī)排列、重復(fù)、局部控制、正交等原則設(shè)計(jì)一個(gè)試驗(yàn)，除要考察的因素外，其他試驗(yàn)條 件均要求一致，即要考察的試驗(yàn)因素的試驗(yàn)結(jié)果都是同協(xié)方差陣的且相互獨(dú)立的多元正態(tài)總體。 因而，各因素對(duì)試驗(yàn)結(jié)果影響的結(jié)果的比較，就變成了多個(gè)同協(xié)方差陣的多元正態(tài)總體均

12、值向 量的比較。統(tǒng)計(jì)上解決兩個(gè)以上同協(xié)方差陣多元正態(tài)總體均值向量比較的方法叫做多元方差分 析。多個(gè)總體均值向量的比較檢驗(yàn)，特別是多元方差分析正是本節(jié)的內(nèi)容，這類方法在經(jīng)濟(jì)管 理，系統(tǒng)控制，生物醫(yī)藥等許多領(lǐng)域有著廣泛的應(yīng)用。這里先看一個(gè)具體實(shí)例。341提出問題例3.3為了研究某種疾病，對(duì)三組人測(cè)量：第1組是20至35歲女性、第2組是20至25歲男性、第3組是30至55歲男性。每組取20個(gè)人，測(cè)量第I組的第J人4個(gè)指標(biāo)是：一：脂蛋白 （X；）、甘油三脂（X；2）、脂蛋白（X（3）、前0脂蛋白（X（4 ）。測(cè)量結(jié)果見表3-3。 問三組人的指標(biāo)間有沒有顯著差別？表3-3脂蛋白、甘油三脂、丄脂蛋白、前脂

13、蛋白數(shù)據(jù)(1)Xj；(1)Xj2(1)Xj3(1)Xj4(2) Xj1j2)Xj2j2)Xj3j2)Xj4j3)Xj1j3)Xj2j3)Xj3j3)Xj4260754018310122302132064391720072341731060351826059371124087451819040271536088282617065391722565341629510036122701103924170653716270653221205130342321082311738011436211906927 1528067 :371824055P 4210200464515210383617260553

14、42025011721202806530232601102920200107282020076401729573332122513036 11120076 3920240114P 3818210125261728094261131010332181706431 11419060 13317330112211127076331329555301634512724201906034 16270125 :242125062P 22162808120182801203218260592119310119251524062 343027057318280692920345120

15、3618250673114370703020360107252326013539292804037172501173616問題中的3組人的測(cè)量值X、X、X，每個(gè)隨機(jī)向量有4個(gè)指標(biāo)，即4維隨機(jī) 向量。例3.3要從每個(gè)總體20個(gè)樣品值出發(fā)，檢驗(yàn) E(X)=E(X(2) = E(X(3)是否成立。3.4.2單因素方差分析的數(shù)學(xué)模型方差分析的目的在于找出自變量與因變量之間的線性關(guān)系，或自變量對(duì)因變量的實(shí)驗(yàn)效果。方差分析是一種處理實(shí)驗(yàn)數(shù)據(jù)的方法，考察一個(gè)被稱為因變量或相依變量(depe nde nt variable,)的連續(xù)響應(yīng)變量，又稱反應(yīng)變量(Response Variable)，其數(shù)值則是連續(xù)的

16、，它在由分類變量識(shí)別的幾種試驗(yàn)條件下被測(cè)量，這些分類變量被稱為自變量，獨(dú)立變量(independent variable)，定性變量(Qualitative Variable )或分類變量(Classification Variable)，其數(shù)值多半是不連續(xù)的。這些分 類變量的水平組合形成試驗(yàn)設(shè)計(jì)的單元。例如，某個(gè)試驗(yàn)要測(cè)量男人和女人的重量變化(因變 量)，他們采取了三種不同的減肥方法，這個(gè)設(shè)計(jì)的6個(gè)單元由性別(男、女)和減肥方法(A、B、C) 6種組合形成。一項(xiàng)試驗(yàn)有多個(gè)影響因素，因素也可以看成是一種變量，其取值不是數(shù)，而是水平。例如 “產(chǎn)地”是一個(gè)變量，它取的值是“北京”、“上?！?、“南京

17、”等。這種變量稱為屬性變量，定性變量或分類變量如果只有一個(gè)因素在發(fā)生變化，其他因素保持不變，貝U稱為單因素試驗(yàn)， 與之對(duì)應(yīng)的方差分析，稱為單因素方差分析。我們所考察的影響產(chǎn)品指標(biāo)的因素(如產(chǎn)地，溫度)也稱為因子，用大寫字母 A,B,C表示。因素所能處的狀況，如甲、乙、丙；60, 65, 70, 75,，稱為因素的水平，簡(jiǎn)稱為水平。水平常以A, A2, B1, B2,.表示。一般地，假設(shè)因素 A有k個(gè)水平：A,.,Ak。對(duì)第i個(gè)水平A進(jìn)行試驗(yàn)，獨(dú)立觀察 ni次， i=1,2，k，整個(gè)試驗(yàn)共作了 n1 n2 nk = n次，且完全隨機(jī)排列。設(shè)A的第j次觀察的試驗(yàn)指標(biāo)為p維向量X(i) =(x(1)

18、,x(2),，X(ip) Np(Ui2) i=1,2，k , j=1,2，ni 假設(shè)：(1) 同一個(gè)水平Aj下得到的觀測(cè)值X1(1), x2)xni ;X1(k)，X，由于實(shí)驗(yàn)過程中各種偶然因素的干擾及測(cè)量誤差所致，每次實(shí)驗(yàn)中這些偶然因素的總和稱為實(shí)驗(yàn)誤差，它們是方差相同的零均值正態(tài)隨機(jī)變量；(2) 所有誤差相互獨(dú)立；(3) 由于水平的不同，可能會(huì)給Xj一個(gè)定量的確定性的影響，其大小是未知的。假定=1 0譏令：.j =丄n i于是有模型：；ij1 Np(0j)且相互獨(dú)立i =1,2, ,k , j -1,2,其中卩稱為總體均值向量， 為A的主效應(yīng)向量， 為a的第j次觀察的隨機(jī)誤差向量，根 據(jù)假

19、設(shè)相互獨(dú)立且均服從 Np(0,匕)。判斷這個(gè)因素的影響是否顯著就是要檢驗(yàn)假設(shè):Ho ： 1 = 2 八八 k =0 Hi : 1 , 2 ,k 不全為 0( 3.7)1 ni設(shè)第I組樣本均值X 丄 x(i)ni j41 kni總均值Xx(i)n i 4 j4kni 樣本組內(nèi)差EM二(X； -茹)以- 丁)i 1 j 1k 樣本組間差 B = 7 n/X - X)(X - X(i),i 4k niA 八 、(Xj -XHXj - X)= B E，i =1 j 呂對(duì)于該檢驗(yàn)問題的統(tǒng)計(jì)量，取 WILKS統(tǒng)計(jì)量上 二 E/A定理3.3若打則上二E / A服從WILKS分布二p, n _k.1,k _2

20、x(2)x,x證明參見朱道元第177頁例3.3為了研究某種疾病，對(duì)三組人測(cè)量：第 1組是20至35歲女性、第2組是20至25歲男 性、第3組是30至55歲男性。每組取 20個(gè)人，測(cè)量第I組的第J人4個(gè)指標(biāo)是：一：脂蛋白(i)(i)(i)(x(7 )、甘油三脂(X；?)、脂蛋白(x；3 )、前：脂蛋白(Xj4 )。測(cè)量結(jié)果見表3.3。問三組人的指標(biāo)間有沒有顯著差別?x(2)x,xx(2)x,x解這兒有3個(gè)總體，建立假設(shè) H 0 ： S二篤八1不全相等x(2)x,xx(2)x,x計(jì)算三總體樣本均值231.0253.5292.7589.672.5590.232.932.4531.75.17.1 一.

21、17.9 一-18.4 一x(2)x,xx(2)x,x計(jì)算組內(nèi)差x(2)x,x305306298-107819515736.8 -796.8 1387.8 Ei =I .95590.2413.8517057021.5-1571.582712288.95 -807.95 321.1 E2 =.364.95-5.1133.843173.759959-1301.2572312441.2-333457.4E3 =.761.5-112476.812504.823278.5-395.75174840466.395 -1937.752166.3E =巳 +E2 + Er=.2082.5-26.9.1024.

22、2計(jì)算組間差39065.832307.92-724.08786 14017.23-35.82-26.9B =13.43-14.7J17.2 一計(jì)算總方差-164474.5825586.4244484.18-4674.83-1973.572534 12139.4A=B+E=.2095.93-41.61041.4 一計(jì)算上統(tǒng)計(jì)量一i = |E|/帆二0.6621 ，查得 A p,n4 4k/ (0.01) = 0.7090.6621 ；所以高度顯著否定 H。，故三組人身體指標(biāo)有顯著差異。3.5總體協(xié)差陣相等的檢驗(yàn)本章第三節(jié)和第四節(jié)中， 總假定不同總體的方差是相同的，這一假定是否合理？在一些問 題中

23、應(yīng)當(dāng)加以證明。3.5.1個(gè)正態(tài)總體協(xié)方差陣的檢驗(yàn)設(shè)Xi,X2,，Xn為取自P維正態(tài)總體Np(d2)的一個(gè)樣本，二未知，且三0。首先，我們考慮假設(shè)檢驗(yàn)問題:Ho : = I p， HiIp所構(gòu)造的檢驗(yàn)統(tǒng)計(jì)量為一 expQtrA A n2(np2其中nTA 八(Xj _X)(Xj _X)i A然后，我們考慮假設(shè)檢驗(yàn)問題：Ho :二匸0 = I p， H1 ： ：匸0 = I p因?yàn)閆o0，所以存在非奇異矩陣D，使得DoDT = Ip令丫 =DXj, i =1,2，,n則Y NP(D),D3Dt ) =NP(H )因此檢驗(yàn)二二70等價(jià)于匕=I p此時(shí)構(gòu)造檢驗(yàn)統(tǒng)計(jì)量為階沖* A*.苫嚴(yán)其中n*-、TA

24、=v (Yj -Y)(Yj -Y)i A給定檢驗(yàn)水平,因?yàn)橹苯佑械姆植加?jì)算臨界值很困難，所以通常采用的近似分布。在H 0成立時(shí)，- 21n 的極限分布是2 (p( p 1). 2)，因此當(dāng)n p ,由樣本值計(jì)算出值，若-21 n 2.，即 e_：2，則拒絕Ho，否則不能拒絕 Ho。3.5.2多個(gè)協(xié)方差陣相等檢驗(yàn)剛才討論的檢驗(yàn)、-lo是一個(gè)正態(tài)總體協(xié)方差陣的檢驗(yàn)，是檢驗(yàn)當(dāng)前協(xié)方差陣與過去是否一樣，在一些實(shí)際問題中，可能會(huì)遇到多個(gè)正態(tài)總體的協(xié)方差陣是否相等的問題。設(shè)有k個(gè)正態(tài)總體分別為 Np(叫,Zi ),，Np(k ,6)，匕i - 0且未知，i = 1,2，k從第i個(gè)總體中取ni個(gè)樣本Xj=

25、(x(1),x(2)，x(p) Np(Ui Ji ), j這里n1 n2川卷nk =門為總樣本容量。=1 ,2, ni我們考慮假設(shè)檢驗(yàn)問題為H o : 1=匸2=_=k，構(gòu)造檢驗(yàn)統(tǒng)計(jì)量為Hi仁i 不全相等np 2i AAn2k p：2丨ni 4山A (Xj(i)_X(i)(X(i)_X(i)T，j 二(i)X1ni二丄V X(i) j ?ni j m按照Bartlett的建議，記 Nj二山kni 二血 = NANi 2NNp2k廠 NiNip2i=1-1, N = n - k得到修正的檢驗(yàn)統(tǒng)計(jì)量則在Ho成立時(shí),2(f)，其中1-2lnr的極限分布是一1 d1f =2 p(p 1)(k T)汁+

26、3p-1眩1 1 I Nk環(huán)必相等; d =6(p+1)(k-1) #Ni N!(2p3p-1)(k+1),也廂等I6(p+1)Nk例3.4有甲、乙兩品種，取得如表3-4所示的兩個(gè)二元正態(tài)樣本，試檢驗(yàn) H。： 1=匕2觀察值和Z Xi2工 X1X23002322171002863201455385109甲X23525431017123341726085 (2)200 :150r 333150283383350P 3002149635167乙X1616381* Z v5043834173808610055642044表3-4方差陣檢驗(yàn)數(shù)據(jù)&32271.5-3742.51解：A =-3742.58

27、95.5 一90163.4 8540.0 A 8540.04297.5A =14892822A2 =460883525A =314545504.17 8991 1 2 2 852A I2|1 2 2 852 3 4 0.2In Aj =16.5164In A2 =17.6461In A =19.5666由于 p =2，k =2，-2In 丨 N In A-p I nN 亠 p N k I nN k N k IniTiTN -12，故Ak=175.1614 43.3371 -206.0957 =7.815d =0.18741匕p(p皿仆311 -d0.052(3) = 9.6176由于-2ln

28、6176，故應(yīng)拒絕H0 :72，即認(rèn)為11 ,匕2有顯著差異。3.6獨(dú)立性檢驗(yàn)X (1) | q一個(gè)隨機(jī)向量 X = | (2)，若其中兩子向量相互獨(dú)立，則可化為兩個(gè)低維隨即向量x(2處理，給統(tǒng)計(jì)分析帶來極大的便利，因此檢驗(yàn)一個(gè)隨機(jī)向量的子向量之間是否獨(dú)立是參數(shù)假設(shè)檢驗(yàn)中的重大課題，而當(dāng)陟Np隠冒-X 一 I 尸一二 21-12 時(shí)，X,x(2)相互獨(dú)立口 X,x(2)互不相22J關(guān)二712 =0(二 二21 =0)。這時(shí)，X (1) , X的獨(dú)立性檢驗(yàn)可歸結(jié)為參數(shù)假設(shè)檢驗(yàn)。般情況下，設(shè) XNp(1) ,3正定，將X分割成k個(gè)子向量:X =(x(1),x(2),，x(k) )T,其中X的維數(shù)為

29、Pi , ik= 1,2，kPi = p，將與三也作相應(yīng)的剖分:i =4亠工11V己121 k卩(2)9,龍=21S22a2kmHkl1 1k1Ik2kk _必卩=檢驗(yàn)子向量x(1)，,X(k)之間的相互獨(dú)立的假設(shè)問題可寫成:i , j 2, kH1 : j =0至少有一對(duì)i, j也就是說，如果H。成立，則00-22-kk現(xiàn)從總體X Np(丄，0抽取容量為n的樣本X1 ,X2，Xn ,將樣本的總離差陣nA (x X)(Xi X)T，剖分成i 4a11a12a1 pA11A-m12Aika21a22sa2pa=A21A22A2ka.a p1a p2app 1Ak1A-Mk2Akk _A 二P1P2Pk也可以計(jì)算樣本相關(guān)陣，并作相應(yīng)剖分112r1 pR11R12R1kr21i1 亠ar2p費(fèi)R21aR22aR2karp1rp2 1px：pRk1 a SRk2Rkk _PlP2Pk其中rij= aij 、aiiajj檢驗(yàn)問題(3.6.1)式的似然比統(tǒng)計(jì)量為IA通常取v|AlR 作為統(tǒng)計(jì)量n Aiii z!n Riii z1關(guān)于在Ho成立