第九章方差分析

第九章方差分析（ANOVA）本文檔共72頁；當前第1頁；編輯于星期三\10點56分方差分析又稱為變異分析（analysisofvariance，ANOVA），是由斯內(nèi)德克（GeorgeWaddelSnedecor）提出的一種方法。方差分析通過對多組平均數(shù)的差異進行顯著性檢驗，分析實驗數(shù)據(jù)中不同來源的變異對總變異影響的大小。

本文檔共72頁；當前第2頁；編輯于星期三\10點56分t檢驗法適用于樣本平均數(shù)與總體平均數(shù)及兩樣本平均數(shù)間的差異顯著性檢驗，但在生產(chǎn)和科學研究中經(jīng)常會遇到比較多個處理優(yōu)劣的問題，即需進行多個平均數(shù)間的差異顯著性檢驗(K≥3)。這時，若仍采用t檢驗法就不適宜了。這是因為：

本文檔共72頁；當前第3頁；編輯于星期三\10點56分

1、檢驗過程煩瑣例如，一試驗包含5個處理，采用t檢驗法要進行10次兩兩平均數(shù)的差異顯著性檢驗；若有k個處理，則要作k(k-1)/2次類似的檢驗。本文檔共72頁；當前第4頁；編輯于星期三\10點56分2、無統(tǒng)一的試驗誤差，誤差估計的精確性和檢驗的靈敏性低對同一試驗的多個處理進行比較時，應(yīng)該有一個統(tǒng)一的試驗誤差的估計值。若用t檢驗法作兩兩比較，由于每次比較需計算一個s，故使得各次比較誤差的估計不統(tǒng)一，同時沒有充分利用資料所提供的信息而使誤差估計的精確性降低，從而降低檢驗的靈敏性。本文檔共72頁；當前第5頁；編輯于星期三\10點56分

例如，試驗有5個處理，每個處理重復(fù)6次，共有30個觀測值。進行t檢驗時，每次只能利用兩個處理共12個觀測值估計試驗誤差，誤差自由度為2(6-1)=10；若利用整個試驗的30個觀測值估計試驗誤差，顯然估計的精確性高，且誤差自由度為5(6-1)=25?？梢姡谟胻檢法進行檢驗時，由于估計誤差的精確性低，誤差自由度小，使檢驗的靈敏性降低，容易掩蓋差異的顯著性。

本文檔共72頁；當前第6頁；編輯于星期三\10點56分3、這種兩兩比較會隨著樣本組數(shù)的增加而加大犯Ⅰ型錯誤的概率用t檢驗法進行4個處理平均數(shù)間的差異顯著性檢驗，若兩兩比較推斷正確的概率為95%,則所有比較都正確的概率為0.956=0.74,則降低了推斷的可靠性。

本文檔共72頁；當前第7頁；編輯于星期三\10點56分幾個常用術(shù)語:1、試驗指標(experimentalindex)為衡量試驗結(jié)果的好壞或處理效應(yīng)的高低，在試驗中具體測定的性狀或觀測的項目稱為試驗指標。由于試驗?zāi)康牟煌x擇的試驗指標也不相同。本文檔共72頁；當前第8頁；編輯于星期三\10點56分2、試驗因素(experimentalfactor)試驗中所研究的影響試驗指標的因素叫試驗因素。當試驗中考察的因素只有一個時，稱為單因素試驗；若同時研究兩個或兩個以上的因素對試驗指標的影響時，則稱為兩因素或多因素試驗。試驗因素常用大寫字母A、B、C、…等表示。

本文檔共72頁；當前第9頁；編輯于星期三\10點56分

3、因素水平(leveloffactor)試驗因素所處的某種特定狀態(tài)或數(shù)量等級稱為因素水平，簡稱水平。本文檔共72頁；當前第10頁；編輯于星期三\10點56分4、試驗處理(treatment)事先設(shè)計好的實施在試驗單位上的具體項目叫試驗處理，簡稱處理。在單因素試驗中，實施在試驗單位上的具體項目就是試驗因素的某一水平。進行單因素試驗時，試驗因素的一個水平就是一個處理。在多因素試驗時，試驗因素的一個水平組合就是一個處理。本文檔共72頁；當前第11頁；編輯于星期三\10點56分5、重復(fù)(repetition)在試驗中，將一個處理實施在兩個或兩個以上的試驗單位上，稱為處理的重復(fù)數(shù)。本文檔共72頁；當前第12頁；編輯于星期三\10點56分第一節(jié)方差分析的基本原理及步驟本文檔共72頁；當前第13頁；編輯于星期三\10點56分一、方差分析的基本原理：綜合的F檢驗（一）綜合的虛無假設(shè)與部分虛無假設(shè)1.綜合的虛無假設(shè)樣本所歸屬的總體的平均數(shù)相等，即

H0：μ1=μ2=μ32.此為部分虛無假設(shè)

組間的虛無假設(shè)本文檔共72頁；當前第14頁；編輯于星期三\10點56分（二）方差的可分解性（可加性）1.可加性方差分析作為一種統(tǒng)計方法，是把實驗數(shù)據(jù)的總變異分解為若干個不同來源的分量。因而它所依據(jù)的基本原理是變異的可加性。本文檔共72頁；當前第15頁；編輯于星期三\10點56分即每一個數(shù)據(jù)與數(shù)據(jù)的總體平均數(shù)差的平方和，可以分解為每一組數(shù)據(jù)各自的離差平方和與由各組數(shù)據(jù)的平均數(shù)組成的一組數(shù)據(jù)的離差平方和兩部分。前者表達的是組內(nèi)差異，即每組數(shù)據(jù)中各個數(shù)據(jù)之間的差異，也就是個體差異，表達的是抽樣誤差或隨機誤差程度；后者表達的是組間差異，即各組平均數(shù)之間的差異，表達的是實驗操縱的差異程度，實驗操縱即指自變量的操縱，這兩部分差異之間相互獨立。本文檔共72頁；當前第16頁；編輯于星期三\10點56分可用公式表示為：SST=SSB+SSw如：欲觀察某因素的三個水平對被試是否產(chǎn)生相同的影響：組一：A、A、A、A、A—水平一組二：B、B、B、B、B—水平二組三：C、C、C、C、C—水平三總組：A、A、A、A、A、B、B、B、B、B、C、C、C、C、C本文檔共72頁；當前第17頁；編輯于星期三\10點56分2.總體變異的構(gòu)成總體變異組間變異：組內(nèi)變異：組內(nèi)變異理論上要求齊性，實際計算取其均值本文檔共72頁；當前第18頁；編輯于星期三\10點56分3.方差的基本公式一般總體方差稱方差，樣本方差稱均方能使變量發(fā)生變異的原因很多，這些原因我們都將其稱為變異因素或變異來源。本文檔共72頁；當前第19頁；編輯于星期三\10點56分方差分析就是發(fā)現(xiàn)各類變異因素相對重要性的一種方法方差分析的思路就是：把整個試驗（設(shè)有k個總體）的樣本資料作為一個整體來考慮。把整個試驗的總變異按照變異的來源分解成不同因素的變異。由于方差等于平方和除以自由度，因此總方差分解成各因素的方差，就是將形成總方差的平方和和自由度分解為各因素的平方和和自由度。然后對各個因素的方差作出數(shù)量上的估計，從而發(fā)現(xiàn)各個因素的方差的相對重要程度。本文檔共72頁；當前第20頁；編輯于星期三\10點56分從總方差中除去各可控因素所引起的方差后，剩余方差又可以準確地估計試驗誤差，作為統(tǒng)計假設(shè)檢驗的依據(jù)因此，方差分析可以幫助我們抓住試驗的主要矛盾和技術(shù)關(guān)鍵，發(fā)現(xiàn)主要的變異來源，從而抓住主要的、實質(zhì)性的東西。本文檔共72頁；當前第21頁；編輯于星期三\10點56分4.平方和的剖分本文檔共72頁；當前第22頁；編輯于星期三\10點56分本文檔共72頁；當前第23頁；編輯于星期三\10點56分本文檔共72頁；當前第24頁；編輯于星期三\10點56分本文檔共72頁；當前第25頁；編輯于星期三\10點56分本文檔共72頁；當前第26頁；編輯于星期三\10點56分SS（sumofsquares）表示平方和SST（thesumofsquarestotal）總平方和，一個試驗中的總變異。SSB（sumofsquaresbetweengroups）組間平方和，表示由于不同的實驗處理而造成的變異。（主試）SSW（sumofsquareswithingroup）試驗誤差造成的變異。（被試）本文檔共72頁；當前第27頁；編輯于星期三\10點56分

SST=SSB+SSW總變異=組間變異+組內(nèi)變異總變異（SST）是將所有被試的數(shù)值作為一個整體考慮到的結(jié)果，是用所有被試的因變量的值計算出來的。組間變異（SSB

）主要是因為接受不同的實驗處理而造成的各組之間的變異。用兩個平均數(shù)的離差表示。組內(nèi)變異（SSW）指組內(nèi)各被試因變量的差異范圍。本文檔共72頁；當前第28頁；編輯于星期三\10點56分5.組間方差（組間均方）與組內(nèi)方差（組內(nèi)均方）本文檔共72頁；當前第29頁；編輯于星期三\10點56分6.自由度的分解在計算處理間平方和時，各處理均數(shù)要受這一條件的約束，故處理間自由度為處理數(shù)減1，即k-1。處理間自由度記為dfB，即dfB=k-1。在計算處理內(nèi)平方和時，每組自由度為n-1,共有k組，故處理內(nèi)自由度k(n-1),

處理內(nèi)自由度記為dfW，即dfW=k(n-1)=kn-k。本文檔共72頁；當前第30頁；編輯于星期三\10點56分本文檔共72頁；當前第31頁；編輯于星期三\10點56分7、方差分析方差之間的差異分析用F檢驗，因此，組間與組內(nèi)方差的分析也用F檢驗。方差分析關(guān)注的是組間均方是否顯著大于組內(nèi)均方。因此，常用作單側(cè)檢驗。本文檔共72頁；當前第32頁；編輯于星期三\10點56分kj表明數(shù)據(jù)的總變異基本上是有不同的實驗處理造成的（不同的實驗處理間存在顯著差異）表明數(shù)據(jù)的總變異基本上是有實驗誤差和個體誤差造成的，與不同的實驗處理關(guān)系不大（不同的實驗處理間不存在顯著差異）。本文檔共72頁；當前第33頁；編輯于星期三\10點56分二、方差分析的基本假設(shè)1．總體正態(tài)分布2．各實驗處理是隨機的且相互獨立（一般情況下都能滿足）3．各實驗處理內(nèi)方差一致（需要進行檢驗）最為重要的假定本文檔共72頁；當前第34頁；編輯于星期三\10點56分三、方差齊性檢驗1.哈特萊最大F比率法找出要比較的組內(nèi)方差的最大值與最小值。最大方差與最小方差無顯著差異即為方差齊性。本文檔共72頁；當前第35頁；編輯于星期三\10點56分四、方差分析的基本步驟本文檔共72頁；當前第36頁；編輯于星期三\10點56分（一）．提出假設(shè)（二）．選擇檢驗統(tǒng)計量并計算

1．分解平方和總平方和SST

組間平方和SSB

組內(nèi)平方和SSW本文檔共72頁；當前第37頁；編輯于星期三\10點56分

2．分解自由度df總自由度：dfT=nk-1組間自由度：dfB=k-1組內(nèi)自由度：dfW=nk-k組內(nèi)自由度的計算（1）不同實驗處理人數(shù)相同時每組自由度n1-1，n2-1，n3-1…nk-1組內(nèi)自由度n1-1+n2-1+n3-1+…nk-1，因為n1=n2=n3=…nk，所以為K（n-1）（2）不同實驗處理人數(shù)不相同時每組自由度n1-1，n2-1，n3-1…nk-1組內(nèi)自由度n1-1+n2-1+n3-1+…nk-1本文檔共72頁；當前第38頁；編輯于星期三\10點56分

3．計算方差MSMSB=SSB/dfB和MSW=SSW/dfW4．計算F值F=MSB/MSW（三）．作出統(tǒng)計結(jié)論1.顯著性水平2.臨界水平F值3.比較推論本文檔共72頁；當前第39頁；編輯于星期三\10點56分（四）陳列方差分析表變異（差異）來源平方和（SS）自由度（df）

均方(MS) Fp組間SSBdfB=k-1MSBFp組內(nèi)SSWdfW=nk-kMSW總變異SSTdfT=nk-1MST本文檔共72頁；當前第40頁；編輯于星期三\10點56分例：研究人員采用四種不同的心理治療方案，對每個志愿參加治療的患者進行心理治療。他們用錄音機記錄了每個被試在一段時間中所講的詞數(shù)。由于錄音的困難每種方案記錄的人數(shù)各不相同，原始數(shù)據(jù)見表1。問這幾種方案是否有差異？本文檔共72頁；當前第41頁；編輯于星期三\10點56分序號治療方案X1X2X3X41309005025001832488774427454763814445631367860843462116664356341156603600458336462384424576765776562384444193666435663814445833645227047806400計算表1本文檔共72頁；當前第42頁；編輯于星期三\10點56分1．提出假設(shè)H0：μ1＝μ2＝μ3＝μ4H1：至少有兩個總體平均數(shù)不等2．選擇檢驗統(tǒng)計量并計算假定四組記錄數(shù)據(jù)是從四個正態(tài)總體中抽出的獨立樣本，對多個獨立樣本平均數(shù)的差異進行顯著性檢驗，應(yīng)采用完全隨機設(shè)計的方差分析。本文檔共72頁；當前第43頁；編輯于星期三\10點56分本文檔共72頁；當前第44頁；編輯于星期三\10點56分(1).計算平方和：組間平方和組內(nèi)平方和總平方和本文檔共72頁；當前第45頁；編輯于星期三\10點56分(2)．計算自由度組間自由度組內(nèi)自由度總自由度本文檔共72頁；當前第46頁；編輯于星期三\10點56分(3)．計算方差組間方差

組內(nèi)方差(4)．計算Ｆ值

本文檔共72頁；當前第47頁；編輯于星期三\10點56分3．做統(tǒng)計決斷,列方差分析表變異來源平方和自由度方差F值概率組間變異2850.43950.13.77*P<0.05組內(nèi)變異4786.519251.9總變異7636.922表9－3四組記錄數(shù)據(jù)的完全隨機設(shè)計方差分析表本文檔共72頁；當前第48頁；編輯于星期三\10點56分第二節(jié)完全隨機設(shè)計的方差分析本文檔共72頁；當前第49頁；編輯于星期三\10點56分單因素完全隨機設(shè)計的方差分析(completerandomizeddesign)把從同一個總體中隨機抽取的被試隨機地分配為幾個實驗組，每個實驗組分別接受某一因素的不同水平的操縱，檢驗這幾個獨立樣本平均數(shù)之間是否存在顯著差異本文檔共72頁；當前第50頁；編輯于星期三\10點56分一、各實驗處理組樣本容量相同例題：從五所中學同一個年級隨機抽取3名學生進行統(tǒng)一數(shù)學測驗，結(jié)果如下，問五所學校數(shù)學成績之間有無顯著差異？

ABCDE176788683732738184827437081858778本文檔共72頁；當前第51頁；編輯于星期三\10點56分實驗分析：自變量：學校類型自變量的不同水平：不同類型的學校因變量：數(shù)學成績實驗?zāi)康募邦A(yù)學校類型對數(shù)學成績有無影響？重點中學應(yīng)該好于一般中學本文檔共72頁；當前第52頁；編輯于星期三\10點56分1.提出假設(shè)：虛無假設(shè)：u1=u2……=u5；備選假設(shè)：至少兩個總體的平均數(shù)不相等；2.選擇并計算統(tǒng)計量（1）計算平方和容量平均數(shù)總和平方和均方組一3732191600515987組二380240920619200組三3852552166721168組四3842522118216899組五3752252116816875合計1511919495994905本文檔共72頁；當前第53頁；編輯于星期三\10點56分組間平方和：SSB=94905-11912/3*5=339.6組內(nèi)平方和：SSW=94959-94905=54（2）計算自由度：組間自由度：dfB=k-1=-5-1=4組內(nèi)自由度：dfW

=k(n-1)=5×2=10本文檔共72頁；當前第54頁；編輯于星期三\10點56分（3）計算方差組間方差：MSB=SSB/(k-1)=339.6/4=84.9組內(nèi)方差：MSW

=SSW

/(nk-n)=54/10=5.4（4）計算F值F=MSB/MSW=84.9/5.4=15.723.查F值表做出統(tǒng)計推斷F0.05（dfB，dfW

）=F0.05（4，10）

=3.48統(tǒng)計推斷：F=15.72>F0.05（dfB，dfW）=3.48,p<0.05。存在顯著差異。本文檔共72頁；當前第55頁；編輯于星期三\10點56分4.列出方差分析表本文檔共72頁；當前第56頁；編輯于星期三\10點56分二、各實驗處理組樣本容量不同例：研究人員采用四種不同的心理治療方案，對每個志愿參加治療的患者進行心理治療。他們用錄音機記錄了每個被試在一段時間中所講的詞數(shù)。由于錄音的困難每種方案記錄的人數(shù)各不相同，原始數(shù)據(jù)見表1。問這幾種方案是否有差異？本文檔共72頁；當前第57頁；編輯于星期三\10點56分1．提出假設(shè)H0：μ1＝μ2＝μ3＝μ4H1：至少有兩個總體平均數(shù)不等2．選擇檢驗統(tǒng)計量并計算假定四組記錄數(shù)據(jù)是從四個正態(tài)總體中抽出的獨立樣本，對多個獨立樣本平均數(shù)的差異進行顯著性檢驗，應(yīng)采用完全隨機設(shè)計的方差分析。本文檔共72頁；當前第58頁；編輯于星期三\10點56分本文檔共72頁；當前第59頁；編輯于星期三\10點56分(1).計算平方和：組間平方和組內(nèi)平方和總平方和本文檔共72頁；當前第60頁；編輯于星期三\10點56分(2)．計算自由度組間自由度組內(nèi)自由度總自由度本文檔共72頁；當前第61頁；編輯于星期三\10點56分(3)．計算方差組間方差

組內(nèi)方差(4)．計算Ｆ值

本文檔共72頁；當前第62頁；編輯于星期三\10點56分3．做統(tǒng)計決斷,列方差分析表變異來源平方和自由度方差F值概率組間變異2850.43950.13.77*P<0.05組內(nèi)變異4786.519251.9總變異7636.922表9－3四組記錄數(shù)據(jù)的完全隨機設(shè)計方差分析表本文檔共72頁；當前第63頁；編輯于星期三\10點56分三、利用樣本統(tǒng)計量進行方差分析適用范圍：分析資料沒有原始數(shù)據(jù)，只有si2、ni，平均數(shù)等。適用原始