方差分析包括三因素課件

1如某種農(nóng)作物的收獲量受作物品種、肥料種類及數(shù)量等的影響；選擇不同的品種、肥料種類及數(shù)量進行試驗，日常生活中經(jīng)常發(fā)現(xiàn)，影響一個事物的因素很多，希望找到影響最顯著的因素2看哪一個影響大？并需要知道起顯著作用的因素在什么時候起最好的影響作用。方差分析就是解決這些問題的一種有效方法。3例如：某廠對某種晴棉漂白工藝中酸液濃度（g/k）進行試驗，以觀察酸液濃度對汗布沖擊強力有無顯著影響。序號沖擊強力濃度123456A116.215.115.814.817.115.0

A216.817.517.115.918.417.7A319.020.118.918.220.519.7方差分析就是把總的試驗數(shù)據(jù)的波動分成1、反映因素水平改變引起的波動。2、反映隨機因素所引起的波動。然后加以比較進行統(tǒng)計判斷，得出結(jié)論。5方差分析的基本思想：把全部數(shù)據(jù)關(guān)于總均值的離差平方和分解成幾部分，每一部分表示某因素諸水平交互作用所產(chǎn)生的效應，將各部分均方與誤差均方相比較，從而確認或否認某些因素或交互作用的重要性。用公式概括為：總變異=組間變異+組內(nèi)變異各因素引起由個體差異引起（誤差）種類：常用方差分析法有以下4種1、完全隨機設計資料的方差分析（單因素方差分析）2、隨機區(qū)組設計資料的方差分析（二因素方差分析）3、拉丁方設計資料的方差分析（三因素方差分析）4、R*C析因設計資料的方差分析（有交互因素方差分析）6第二節(jié)單因素方差分析一、假設檢驗設：A1、A2、A3、為三個總體X1、X2、X3，每個總體有6個樣本Xi1、Xi2、…、Xi6(i=1,2,3)。注：要判斷酸液濃度的3種水平對汗布的沖擊強力是否有顯著影響，實質(zhì)上就是檢驗3種不同水平所對應的3個總體是否有顯著差異的問題。即檢驗3個總體數(shù)學期望是否相等。以后就是求解問題，為了說明一般解的公式（方法），如下作一般分析。7二、離差平方和的分解與顯著檢驗記：將Q進行分解：由于9故：下面通過比較QE和QA來檢驗假設H0。在假設H0成立的條件下，可以證明：相互獨立10方法：（檢驗方法）（1）當H0：1=2=…=m成立時。（2）統(tǒng)計量：即：11同樣可推出：2、數(shù)據(jù)的簡化：試驗數(shù)據(jù)經(jīng)過變換數(shù)據(jù)簡化后對F值的計算沒有影響，不會影響檢驗的結(jié)果四、方差分析表方差來源離差平方和自由度F值F0.05F0.01顯著性因素AQAm-1試驗誤差QEm(n-1)總誤差Qmn-113例：前例題1、對數(shù)據(jù)的簡化得下表：序號沖擊強力濃度123456A1-8-19-12-221-20-801454

A2-251-1114714396A32031191235271443820由表中數(shù)據(jù)可算出14計算計算出F值：15五、各水平下試驗次數(shù)不等時的方差分析設第i個水平試驗次數(shù)為ni,則有自由度分別為f=n-1,fE=n-m,fA=m-1.檢驗統(tǒng)計量為式中：17第三節(jié)雙因素方差分析例如：某廠對生產(chǎn)的高速鋼銑刀進行淬火工藝試驗，考察回火溫度A和淬火溫度B兩個因素對強度的影響。今對兩個因素各3個水平進行試驗，得平均硬度見表：BjAi試驗結(jié)果B1（1210‘C）B2（1235’C）B3（1250‘C）A1（280’C）646668A2（300‘C）666867A3（320’C）656768假設：美中不足組合水平下服從正態(tài)分布、互相獨立、方差相等。所需要解決的問題是：所有Xij的均值是否相等。18假設檢驗：1）在假設H0成立的條件下。2）統(tǒng)計量3）給定顯著水平，查表得臨界值4）由樣本觀察值計算FA、FB5）若時，接受H0，因素的影響不顯著。若時，拒絕H0。對因素B同理說明。19SAS系統(tǒng)中區(qū)分兩種情況：1、每組觀測數(shù)據(jù)相等，可用ANOVA過程處理以上四種情形的方差分析。2、若每組觀測數(shù)據(jù)不相等，可用GLM過程處理以上四種情形的方差分析。21均衡數(shù)據(jù)的方差分析（ANOVA過程）過程說明：1、PROCANOVA；2、CLASS變量表；3、MODEL因變量表=效應；4、MEANS效應[/選擇項]；5、ALPHA=p顯著性水平（缺省值為0.05）CLASS和MODEL是必需的，CLASS必須的MODEL之前?？梢允菙?shù)值型和字符型是指因變量與自變量效應，模型如下：1、主效應模型MODELy=abc;(abc是主效應，y是因變量）2、交互模型MODELy=abca*ba*cb*ca*b*c;3、嵌套效應MODELy=abc(ab);4、混合效應模型號MODELy=ab(a)c(a)b*c(a);輸出因變量均數(shù)，對主效應均數(shù)間的檢驗。22例：1、單因素方差分析某勞動衛(wèi)生組織研究棉布、府綢、的確涼、尼龍四種衣料內(nèi)棉花吸附十硼氫量。每種衣料做五次測量，所得數(shù)據(jù)如下。試檢驗各種衣料見棉花吸附十硼氫量有沒有顯著差別？棉布府綢的確涼尼龍2.332.483.064.002.002.343.065.132.932.683.004.612.732.342.662.802.332.223.063.6023組間占總的比例單因素方差分析SAS程序輸出結(jié)果：結(jié)論：在CLASS語句中指出的P值。P《0.003,可得出各衣料組間有非常顯著差異。252、兩因素方差分析例：用4種不同方法治療8名病者，其血漿凝固時間的資料如表，試分析影響血漿凝固的因素。受試者編號（區(qū)組）處理組1234

8.49.49.812.212.815.212.914.49.69.111.29.89.88.89.912.08.48.28.58.58.69.99.810.98.99.09.210.47.98.18.210.0263、三因素方差分析（交互作用不存在）例：五種防護服，由五人各在不同的五天中穿著測定脈搏數(shù)，如表。試比較五種防護服對脈搏數(shù)有無不同。試驗日期受試者甲乙丙丁戊A129.8B116.2C114.8D104.0E100.6B144.4C119.2D113.2E132.8A115.2C143.0D118.0E115.8A123.0B103.8D133.4E110.8A114.0B98.0C110.6E142.8A110.6B105.8C120.0D109.829三因素方差分析SAS程序的輸入：30三因素方差分析SAS程序輸出結(jié)果：結(jié)論：因F=6.80,P《0.0011,故總體有非常顯著差異。其中K因素影響極大F=16.27,P《0.001。因素P、C都無顯著差異。314、有交互因素的方差分析

例：治療缺鐵性貧血病人12例，分為4組給予不同治療，一個月后觀察紅細胞增加（百萬/mm），資料如表。試分析兩種藥物對紅細胞增加的影響。甲藥（A）不用（A0）用（A1）乙藥（B）不用B0用B10.80.90.71.31.21.10.91.11.02.12.22.032有交互因素方差分析SAS程序的輸入：33有交互因素方差分析SAS程序輸出結(jié)果：結(jié)論：因F=98.75,P=0.0001<0.01,故總體有非常顯著的差異，因素A、B、A*B都對紅細胞增加數(shù)有非常大的影響。34非平衡數(shù)據(jù)的方差分析（GLM過程）在SAS/STAT中GLM（GeneralLinearModels)過程分析功能最多。包括：1、簡單回歸（一元）2、加權(quán)回歸3、多重回歸及多元回歸4、多項式回歸5、方差分析（尤其不平衡分析）6、偏相關(guān)分析7、協(xié)方差分析8、多元方差分析9、反應面模型分析10、重復測量方差分析35GLM過程在方差分析中的應用：MODEL語句反映因變量與自變量的模型，其形式：模型說明模型類型MODELY=ABＣ；主效應ＭODELY=ABA*B；交互效應MODELY=ABA(B)；嵌套效應ＭODELY1Y2=AB；多元方差分析MODELY=AX；協(xié)方差分析Ａ，Ｂ，Ｃ是分類變量，Ｘ，Ｙ是連續(xù)型變量。36１、不平衡單因素方差分析例：健康男子各年齡組淋巴細胞轉(zhuǎn)化率（％）如表，問各組淋