實驗四方差分析

1、實 驗 報 告課程名稱：_學(xué)院名稱：_班 級：_姓 名：_學(xué) 號：_ _學(xué)年第 _學(xué)期數(shù) 學(xué) 與 統(tǒng) 計 學(xué) 院 制實驗地點課程類別公共課 專業(yè)課實驗日期實驗編組第 組實驗所用時間 小時實驗名稱實驗四 方差分析實驗?zāi)康?.使學(xué)生了解SAS系統(tǒng)2.用SAS軟件進(jìn)行方差分析，并掌握其分析步驟。3.掌握Anova過程過程和GLM過程的應(yīng)用，以及結(jié)果的分析。實驗環(huán)境SAS9.1一、實驗準(zhǔn)備1、SAS過程的選擇：均衡設(shè)計、拉丁方設(shè)計、正交設(shè)計等的一元、多元方差分析和重復(fù)測量的方差分析用 Anova過程；不能用Anova過程分析的數(shù)據(jù)就用GLM過程；另外，若需要對各因素間的各水平進(jìn)行兩兩比較，則也用GLM

2、過程。2、結(jié)果分析的方法： 1) 單因素的方差分析： 先看該因素是否是顯著的，若是則看看哪些水平是有顯著性差異的，有必要的話還應(yīng)該找出一個最有利的水平。2) 雙因素或多因素的方差分析： 先看交互作用是否是顯著的，若是則找出一個最有利的組合；若交互作用不顯著，但存在主效應(yīng)是顯著的，則可以對顯著的主效應(yīng)分別按單因素的方法去分析。二、實驗舉例及內(nèi)容（一）實驗舉例例1、2002年12月進(jìn)行的青貯實驗，數(shù)據(jù)為7個不同品種青貯之間的可溶性有機(jī)物含量。玉米的7個不同品種variety分別為：a1高油玉米647a、a2高油玉米115a、a3高油玉米601、a4高油玉米289、a5農(nóng)大80、a6玉米3138、a

3、7特種玉米。每個品種有兩個重復(fù)，實驗數(shù)據(jù)一共有14個觀測值。分析品種玉米青貯之間的可溶性有機(jī)物WSC的含量是否存在顯著差異。varietywscvarietywscvarietywscvarietywsca120.19a34.08a539.38a716.04a121.03a34.40a538.95a715.86a236.06a417.23a630.04a238.97a416.85a628.68解： 程序如下（其中數(shù)據(jù)集為L3.anova5_1）： 主要輸出結(jié)果如下： (圖4.1)（圖4.2）（圖4.3）（圖4.4）（圖4.5）結(jié)果分析：(1) 用過程univariate判斷正態(tài)性 從圖4.1知

4、， Pr < W 的p值為 0.1694，可判斷wsc服從正態(tài)分布(2)用過程ANOVA ，檢驗是否存在顯著性差異 從圖4.2知 ，p <.0001, 可判斷存在顯著性差異(3) 用MEANS variety /duncan，可具體的看哪些水平存在差異從圖4.3知，(a2, a5), a6, a1, (a4,a7), a3 它們之間存在差異。(4) 用MEANS variety /Dunnett ('a1')，可分析它們與對照組a1的差異是否顯著。從圖4.4知，由于檢驗時都出現(xiàn)了*,說明它們與對照組a1的差異是顯著的。(5) 用過程N(yùn)PAR1WAY進(jìn)行非參數(shù)檢驗看是

5、否存在顯著性差異，以及哪個水平是最好的從圖4.5可看出，在0.05的水平下，各水平之間存在差異，且（a5, a2）的水平是最高的。例2 在組織培養(yǎng)過程中采用兩種激素（KT、2,4-D）的不同濃度配比，通過觀測材料的生根比例，選則合適某實驗品種的配比。設(shè)計濃度及觀測到的生根比例如下表： 2,4濃度（d24） 0 0.05 0.10KT濃度（kt）00.050.1000.050.1000.050.10 重復(fù)10.240.350.380.200.250.360.280.500.30 重復(fù)20.260.310.350.210.210.350.200.510.31 重復(fù)30.250.320.350.22

6、0.260.340.210.520.32解： 程序如下：主要輸出結(jié)果如下： （圖4.6） （圖4.7） 結(jié)果分析： (1)由 anlysis of variance中的輸出結(jié)果圖4.6，可得模型及d24, kt, d24*kt都是顯著的。 （2）由于交互作用是顯著的，所以選擇最優(yōu)的組合，結(jié)合圖4.7可知，d24=0.1, kt=0.05的組合。 注：若不存在交互作用，則分別選擇d24和kt的最優(yōu)者即可。（二）實驗題目1.2014年下半年同一苜蓿品種接種五種生防菌進(jìn)行比較實驗（數(shù)據(jù)如下表）。變量說明：deal為處理，期中CK為對照組，y代表根長。解決的問題：同一苜蓿品種接種不同生防菌后根長的差異

7、。dealy dealy dealy dealy dealy dealy XJ8307345NM024-238XJ8306040NMG048-239XJ8309737CK39XJ8307341NM024-233XJ8306038NMG048-238XJ8309736CK50XJ8307344NM024-225XJ8306037NMG048-242XJ8309736CK40XJ8307336NM024-235XJ8306032NMG048-235XJ8309737CK46XJ8307338NM024-240XJ8306042NMG048-247XJ8309740CK47XJ8307339NM02

8、4-241XJ8306040NMG048-238XJ8309735XJ8307341NM024-241XJ8306042NMG048-239XJ8309732XJ8307336NM024-235XJ8306036NMG048-236XJ8309737XJ8307338NM024-246XJ8306039NMG048-231XJ8309738XJ8307345NM024-238XJ8306040NMG048-240XJ8309736XJ8307337NM024-245XJ8306042NMG048-238XJ8309733XJ8307340NM024-240XJ8306039NMG048-237XJ83097422.為了研究某生長素不同濃度對豌豆生長的影響，從而選擇推薦使用濃度而進(jìn)行了此試驗（數(shù)據(jù)見下表）。m1、m2、m3、m4、m5分別代表不同的濃度，即5個處理（treat）。Weight為豌豆植株干重，單位為g/盆。每個處理有4個重復(fù)，共20個觀測值。選擇一