生物統(tǒng)計學答案兩因素及多因素方差分析

1、第九章 兩因素及多因素方差分析9.1 雙菊飲具有很好的治療上呼吸道感染的功效，為便于飲用，制成泡袋劑。研究不同浸泡時間和不同的浸泡溫度對浸泡效果的影響，設(shè)計了一個兩因素交叉分組實驗，實驗結(jié)果(浸出率)見下表52：浸泡溫度/浸泡時間/min1015206023.7225.4223.588024.8428.3229.559530.6431.5832.21對以上結(jié)果做方差分析及Duncan檢驗。該設(shè)計已經(jīng)能充分說明問題了嗎？是否還有更能說明問題的設(shè)計方案？答：無重復二因素方差分析程序及結(jié)果如下：options linesize=76 nodate;data hermed;do temp=1 to 3

2、;do time=1 to 3;input effect ;output;end;end;cards;23.72 25.42 23.5824.84 28.32 29.5530.64 31.58 32.21;run;proc anova;class temp time;model effect=temp time;means temp time/duncan alpha=0.05;run;The SAS SystemAnalysis of Variance ProcedureClass Level InformationClass Levels ValuesTEMP 3 1 2 3TIME 3

3、1 2 3Number of observations in data set = 9The SAS SystemAnalysis of Variance ProcedureDependent Variable: EFFECTSum of MeanSource DF Squares Square F Value Pr > FModel 4 87.0707778 21.7676944 12.56 0.0155Error 4 6.9321778 1.7330444Corrected Total 8 94.0029556R-Square C.V. Root MSE EFFECT Mean0.9

4、26256 4.741881 1.31645 27.7622Source DF Anova SS Mean Square F Value Pr > FTEMP 2 78.7202889 39.3601444 22.71 0.0066TIME 2 8.3504889 4.1752444 2.41 0.2058The SAS SystemAnalysis of Variance ProcedureDuncan's Multiple Range Test for variable: EFFECTNOTE: This test controls the type I comparison

5、wise error rate, notthe experimentwise error rateAlpha= 0.05 df= 4 MSE= 1.733044Number of Means 2 3Critical Range 2.984 3.050Means with the same letter are not significantly different.Duncan Grouping Mean N TEMPA 31.477 3 3B 27.570 3 2C 24.240 3 1The SAS SystemAnalysis of Variance ProcedureDuncan

6、9;s Multiple Range Test for variable: EFFECTNOTE: This test controls the type I comparisonwise error rate, notthe experimentwise error rateAlpha= 0.05 df= 4 MSE= 1.733044Number of Means 2 3Critical Range 2.984 3.050Means with the same letter are not significantly different.Duncan Grouping Mean N TIM

7、EA 28.447 3 3AA 28.440 3 2AA 26.400 3 1從方差分析結(jié)果可以得知，溫度是極顯著的影響因素，時間是不顯著因素。在Duncan檢驗中，溫度的三個水平之間差異是顯著的。時間的三個水平間差異不顯著。本實驗是二因素固定模型設(shè)計，如果設(shè)置重復，會得到兩個因素之間的交互作用（如果存在的話），其結(jié)果能更好地說明問題。以上方差分析的結(jié)果可以歸納成下表：變差來源平方和自由度均方FP溫度(temp)78.720 288 9239.360 144 422.710.006 6時間(time)8.350 488 924.175 244 42.410.205 8誤差6.932 177 8

8、41.733 044 4總和94.002 955 689.2 研究浙江蠟梅大苗移栽技術(shù)，處理方式包括移栽后的不同覆蓋方式和做床方法，統(tǒng)計每100株移栽苗的成活率，結(jié)果見下表53：做床方法精細作床僅挖穴覆蓋方法遮 陰9385未遮陰9081根據(jù)以往經(jīng)驗在覆蓋方法與作床方法之間不存在交互作用，對上述結(jié)果做方差分析。請注意，這里的結(jié)果是百分數(shù)。答：本例需對數(shù)據(jù)做反正弦變換，程序和結(jié)果如下：options linesize=76 nodate;data plum;do cover=1 to 2;do seedbed=1 to 2;input y ;surrate=arsin(sqrt(y/100)*18

9、0/3.14159265;output;end;end;cards;93 8590 81;run;proc anova;class cover seedbed;model surrate=cover seedbed;run;The SAS SystemAnalysis of Variance ProcedureClass Level InformationClass Levels ValuesCOVER 2 1 2SEEDBED 2 1 2Number of observations in data set = 4The SAS SystemAnalysis of Variance Proce

10、dureDependent Variable: SURRATESum of MeanSource DF Squares Square F Value Pr > FModel 2 64.5953445 32.2976722 90390.15 0.0024Error 1 0.0003573 0.0003573Corrected Total 3 64.5957018R-Square C.V. Root MSE SURRATE Mean0.999994 0.027238 0.01890 69.3987Source DF Anova SS Mean Square F Value Pr > F

11、COVER 1 9.4515512 9.4515512 26451.66 0.0039SEEDBED 1 55.1437933 55.1437933 99999.99 0.0001從結(jié)果可以看出，覆蓋方式和做床方式都是極顯著因素。以上結(jié)果可以歸納成下表：變差來源平方和自由度均方FP覆蓋方法(cover)9.451 551 219.451 551 226 451.660.003 9做床方法(seedbed)55.143 793 3155.143 793 399 999.990.000 1誤 差0.000 357 310.000 357 3總 和64.595 701 839.3 為了研究不同NaC

12、l質(zhì)量濃度對小麥愈傷組織生長的影響。 配制質(zhì)量濃度分別為0、0.1、0.3和0.5的NaCl MS培養(yǎng)基，接種15天后，測定每塊愈傷組織平均增重百分率，結(jié)果見下表54：材料名稱NaCl質(zhì)量濃度 / H870634/%G8901/%極早熟/%中國春/%0103.8063.7067.3267.100.199.3156.2752.2452.300.352.2645.0124.1734.300.518.3815.3720.4013.22對上述結(jié)果進行方差分析。作者已經(jīng)給出四種實驗材料都是鹽敏感型小麥，但是不同的小麥品種必定對鹽的抗性不同。也就是說，品種與鹽濃度之間存在交互作用，更完善的實驗應(yīng)當怎樣設(shè)計

13、？答：這是一個固定模型設(shè)計，程序不再給出，結(jié)果如下。The SAS SystemAnalysis of Variance ProcedureClass Level InformationClass Levels ValuesCONCEN 4 1 2 3 4MATERIAL 4 1 2 3 4Number of observations in data set = 16The SAS SystemAnalysis of Variance ProcedureDependent Variable: INCREASESum of MeanSource DF Squares Square F Value

14、 Pr > FModel 6 10411.9081 1735.3180 14.18 0.0004Error 9 1101.6046 122.4005Corrected Total 15 11513.5126R-Square C.V. Root MSE INCREASE Mean0.904321 22.54545 11.0635 49.0719Source DF Anova SS Mean Square F Value Pr > FCONCEN 3 8374.15867 2791.38622 22.81 0.0002MATERIAL 3 2037.74942 679.24981 5.

15、55 0.0196方差分析結(jié)果指出，鹽濃度是極顯著的影響因素，不同基因型的愈傷組織是顯著影響因素。根據(jù)以往的研究工作的經(jīng)驗，鹽濃度與基因型之間可能存在交互作用，最理想的設(shè)計應(yīng)當設(shè)置重復，從總平方和中分離出交互作用平方和，問題可以說明得更確切。以上結(jié)果可以歸納成下表。變差來源平方和自由度均方FP濃度間8 374.158 732 791.386 222.810.000 2品系間2 037.749 43679.249 85.550.019 6誤 差1 101.604 69122.400 5總 和11 513.512 7159.4 為了研究植物的光合作用，設(shè)計了一個實驗。將煙草的兩個變種種植在田間，利

16、用CO2-depletion技術(shù)檢測單位葉面積捕獲CO2的比率。實驗共涉及兩個因素：一個是變種，選用了兩個變種；另一個是抽樣時期，在整個生長季共進行40次田間抽樣。這是一個無重復兩因素實驗設(shè)計，方差分析表如下55：變差來源平方和自由度均方F抽樣時期3.356390.060 433.25*變 種0.015 710.015 78.66*誤 差0.071 1390.001 82總 和2.44379注：*P <0.01。根據(jù)實驗設(shè)計，該設(shè)計是一種什么模型？實驗所涉及的兩個因素屬于哪一種類型的因素？為什么？答：這是一個混合模型實驗。變種是固定因素，抽樣時期是隨機因素。因為實驗沒有設(shè)置重復，在無重復

17、的情況下，三種模型的檢驗統(tǒng)計量是一樣的，不知作者為什么不考慮設(shè)置重復。兩個變種是人為選定的，是固定因素。田間抽樣是隨機抽取的，是隨機因素。9.5 野生型C57BL/6及STAT1-/-型小鼠胰島，在移入四氧嘧啶糖尿病的BALB/c小鼠中之后的存活天數(shù)見下表56：實驗材料養(yǎng) 生 處 理未 處 理IL1ra*IL1ra+CsA野生型C57BL/66 11 11 11 12 1213 13 13 14 14 1415 15 17 10 1412 14 14 15 15 16 21STAT-1/型11 12 13 13 1310 1410 12 14 17 17 23注：*IL-1ra：interle

18、ukin1 receptor antagonist（白介素1受體拮抗物）。*CsA：cyclosporine A（環(huán)孢菌素A）。對上述結(jié)果進行方差分析，判斷兩種類型小鼠的胰島存活天數(shù)差異是否顯著？不同養(yǎng)生處理對移植的胰島存活天數(shù)的影響是否顯著？不同養(yǎng)生處理與不同型小鼠之間是否存在交互作用？答：這是一個重復數(shù)不等的兩因素固定模型實驗，所用程序及計算結(jié)果如下。options linesize=76 nodate;data mouse;infile 'e:dataexr9-5e.dat'do treat=1 to 3;do type=1 to 2;input n ;do repeti

19、t=1 to n;input days ;output;end;end;end;run;proc glm;class treat type;model days=treat type treat*type ;run;The SAS SystemGeneral Linear Models ProcedureClass Level InformationClass Levels ValuesTREAT 3 1 2 3TYPE 2 1 2Number of observations in data set = 37The SAS SystemGeneral Linear Models Procedu

20、reDependent Variable: DAYSSum of MeanSource DF Squares Square F Value Pr > FModel 5 70.1272844 14.0254569 1.67 0.1724Error 31 261.0619048 8.4213518Corrected Total 36 331.1891892R-Square C.V. Root MSE DAYS Mean0.211744 21.43162 2.90196 13.5405Source DF Type I SS Mean Square F Value Pr > FTREAT

21、2 26.5151696 13.2575848 1.57 0.2233TYPE 1 3.2191227 3.2191227 0.38 0.5409TREAT*TYPE 2 40.3929921 20.1964961 2.40 0.1075Source DF Type III SS Mean Square F Value Pr > FTREAT 2 10.5603197 5.2801598 0.63 0.5408TYPE 1 0.0089915 0.0089915 0.00 0.9741TREAT*TYPE 2 40.3929921 20.1964961 2.40 0.1075在方差分析表

22、中我們選用I型可估函數(shù)，從F的顯著性概率可以得出，不論是養(yǎng)生處理、小鼠類型還是兩者的交互作用都是不顯著因素。上述結(jié)果可以歸納成下表：變差來源平方和自由度均方FP處理間26.515 169 6213.257 584 81.570.223 3類型間3.219 122 713.219 122 70.380.540 9處理×類型40.392 992 1220.196 496 12.400.107 5誤 差261.061 904 8318.421 351 8總 和331.189 189 2369.6 野生型C57BL/6及STAT-1-/-型小鼠胰島，在移入自發(fā)糖尿病的NOD#小鼠中之后的存活

23、天數(shù)見下表56：實驗材料養(yǎng)生處理未處理IL1ra*CsA*IL1ra+CsA野生型C57BL/60 0 2 5 511 11 12 1313 15 17 0 5 8 12 12 158 8 8 10 10 11 185 10 11 11 12 16 20STAT-1/型6 10 10 1310 125 13 1410 11 11 1212 13注：# NOD：nonobese diabetic（非肥胖糖尿病）。* IL-1ra：interleukin1 receptor antagonist（白介素1受體拮抗物）。* CsA：cyclosporine A（環(huán)孢菌素 A）。對上述結(jié)果進行方差分析

24、，判斷兩種類型小鼠的胰島存活天數(shù)差異是否顯著？不同養(yǎng)生處理對移植的胰島存活天數(shù)的影響是否顯著？不同養(yǎng)生處理與不同型小鼠之間是否存在交互作用？答：本題與第5題的程序基本一樣，下面只給出計算的結(jié)果。The SAS SystemGeneral Linear Models ProcedureDependent Variable: DAYSSum of MeanSource DF Squares Square F Value Pr > FModel 7 81.7459473 11.6779925 0.55 0.7939Error 39 833.4880952 21.3714896Corrected

25、 Total 46 915.2340426R-Square C.V. Root MSE DAYS Mean0.089317 45.64659 4.62293 10.1277Source DF Type I SS Mean Square F Value Pr > FTREAT 3 67.9880108 22.6626703 1.06 0.3770TYPE 1 8.0410256 8.0410256 0.38 0.5432TREAT*TYPE 3 5.7169109 1.9056370 0.09 0.9656Source DF Type III SS Mean Square F Value

26、Pr > FTREAT 3 60.4851648 20.1617216 0.94 0.4290TYPE 1 8.0762698 8.0762698 0.38 0.5423TREAT*TYPE 3 5.7169109 1.9056370 0.09 0.9656本題的兩個主效應(yīng)和它們的交互作用都是不顯著因素。以上結(jié)果可以歸納成下表：變差來源平方和自由度均方FP處理間67.988 010 8322.662 670 31.060.337 0類型間8.041 025 618.041 025 60.380.543 2處理×類型5.716 910 931.905 637 00.090.965

27、 6誤 差833.488 095 23921.371 489 6總 和915.234 042 6469.7 一項音樂心理學研究，實驗是這樣設(shè)計的：為了避免熟悉的音樂環(huán)境，實驗安排在兩種非典型的音樂練習和演出環(huán)境中進行。一種環(huán)境是在劇場底層敞開的大廳中（環(huán)境A），另一種是在辦公室中（環(huán)境B）。要求實驗參與者學習并回憶所學習的練習曲。學習和回憶包括在相同環(huán)境中（AA，BB）和不同環(huán)境中（AB，BA），評判學習和回憶的得分，從而判斷得分與環(huán)境之間的關(guān)系57。該實驗是一個典型的兩因素交叉分組實驗設(shè)計，方差分析表如下：變差來源平方和自由度均方FP學習環(huán)境180.2671180.2671.4410.275

28、回憶環(huán)境640.2671640.2675.1200.064學習環(huán)境×回憶環(huán)境1 008.60011 008.6008.0650.030誤差750.3336125.056總和2 579.4679問：（1）本實驗共有幾次重復？為什么？（2）本實驗屬于哪一種模型？為什么？（3）本實驗的兩個因素中哪些因素是顯著因素？在本實驗中顯著因素的意義是什么？你可以得到什么結(jié)論？答：(1) 因為本實驗共有4種條件組合，df誤差(組合1重復數(shù)1)(組合2重復數(shù)1)( 組合3重復數(shù)1)( 組合4重復數(shù)1)重復數(shù)46。因此，重復數(shù)6410。(2) 屬固定模型。因為：根據(jù)作者所用的檢驗統(tǒng)計量，由推斷，環(huán)境的水平

29、是人為選定的。(3) 只有交互作用是顯著的。說明音樂的學習是與環(huán)境的兩種特定水平有關(guān)的。結(jié)論：音樂的學習屬本例的兩種特定環(huán)境依賴型記憶。作者的結(jié)論是：音樂的學習屬環(huán)境依賴型記憶。這樣的敘述不夠嚴格，只有隨機模型才能夠得到這樣的結(jié)論，固定模型只能說“音樂的學習屬本例的兩種特定環(huán)境依賴型記憶”。9.8 與上一實驗類似，這次是記憶一段16小節(jié)的鋼琴曲。在同一房間中放置兩臺鋼琴，一臺是平臺式大鋼琴（環(huán)境A），一臺是播音室鋼琴（環(huán)境B）。參與者在一臺鋼琴上學習這段曲譜之后，在同一臺鋼琴上（AA，BB）或不同鋼琴上（AB，BA），回憶這段曲子。根據(jù)回憶的正確性獲得評分57。變差來源平方和自由度均方FP學習

30、環(huán)境22.791122.7980.6520.426回憶環(huán)境0.28310.2830.0080.929學習環(huán)境×回憶環(huán)境1 188.15911 188.15933.9680.000誤差979.3962834.978總和2 190.63731問：（1）本實驗共有幾次重復？為什么？ （2）本實驗與上一實驗比較有什么不同，可以改變結(jié)論的性質(zhì)嗎？答：(1) 總的重復數(shù)為32次。(2) 結(jié)論與上一實驗結(jié)果類似，只能說不同鋼琴的這一環(huán)境所產(chǎn)生的交互作用更顯著。同樣不能把這一結(jié)論推廣到水平總體。9.9 研究318歲健康個體尿中Adrenarche標記物的值。其中兩性24小時尿樣中DHEA*的平均含量

31、*如下58：年齡/a性 別男孩/（mg ·d-1）女孩/（mg ·d-1）340.910.90560.900.99781.081.029101.531.4711121.901.5713142.271.8615162.092.1617182.552.31注：* DHEA：Dehydroepiandrosterone（脫氫表雄酮），是合成人體雌激素，雄激素，以及其他一些人體激素的最基本物質(zhì)。*該值已經(jīng)過對數(shù)變換。用兩因素方差分析判斷不同年齡組和不同性別的DHEA差異是否顯著？答：結(jié)果如下表：The SAS SystemAnalysis of Variance Procedure

32、Class Level InformationClass Levels ValuesAGE 8 1 2 3 4 5 6 7 8SEX 2 1 2Number of observations in data set = 16The SAS SystemAnalysis of Variance ProcedureDependent Variable: DHEASum of MeanSource DF Squares Square F Value Pr > FModel 8 4.90295000 0.61286875 35.44 0.0001Error 7 0.12104375 0.01729

33、196Corrected Total 15 5.02399375R-Square C.V. Root MSE DHEA Mean0.975907 8.247678 0.13150 1.59438Source DF Anova SS Mean Square F Value Pr > FAGE 7 4.84654375 0.69236339 40.04 0.0001SEX 1 0.05640625 0.05640625 3.26 0.1139從計算結(jié)果可以得知，年齡是極顯著因素，性別是不顯著因素。以上結(jié)果可以歸納成下表。變差來源平方和自由度均方FP年齡間4.846 543 7570.692

34、363 3940.040.000 1性別間0.056 406 2510.056 406 253.260.113 9誤 差0.121 043 7570.017 291 96總 和5.023 993 75159.10 嗜乳酸桿菌在體內(nèi)處于一種酸性環(huán)境，一項關(guān)于嗜乳酸桿菌(Lactobacillus acidophilus) IndI在體外模擬環(huán)境中，在不同pH和不同時間的活菌數(shù)（活菌數(shù)/mL）變化情況如下表59：時間/ hpH4.53.52.51.522.40×1091.34×1091.68×1081.08×1081.18×1081.02×

35、;1084.58×1043.24×10447.00×1094.14×1082.24×1091.48×1094.60×1073.98×1075.96×1032.36×10362.10×10102.38×1096.80×1084.88×1091.32×1071.30×1072.10×1031.92×103對表中的數(shù)據(jù)進行方差分析，數(shù)據(jù)是服從泊松分布的。答：對于服從泊松分布的數(shù)據(jù)，應(yīng)進行平方根變換。程序與結(jié)果如下：opti

36、ons linesize=76 nodate;data lacto;infile 'E:dataexr9-10e.dat'do time=1 to 3;do pH=1 to 4;do n=1 to 2;input y ;number=sqrt(y);output;end;end;end;run;proc anova;class time pH;model number=time pH time*pH;run;The SAS SystemAnalysis of Variance ProcedureClass Level InformationClass Levels Values

37、TIME 3 1 2 3PH 4 1 2 3 4Number of observations in data set = 24The SAS SystemAnalysis of Variance ProcedureDependent Variable: NUMBERSum of MeanSource DF Squares Square F Value Pr > FModel 11 2.0166E+10 1.8333E+09 2.86 0.0424Error 12 7.7030E+09 6.4191E+08Corrected Total 23 2.7869E+10R-Square C.V.

38、 Root MSE NUMBER Mean0.723601 96.47277 25336.0 26262.3Source DF Anova SS Mean Square F Value Pr > FTIME 2 1.7452E+09 8.7262E+08 1.36 0.2937PH 3 1.5229E+10 5.0763E+09 7.91 0.0036TIME*PH 6 3.1921E+09 5.3201E+08 0.83 0.5696只有“pH”是極顯著因素，“時間”和“時間×pH”都是不顯著因素。以上結(jié)果可以歸納成下表。變差來源平方和自由度均方FP時 間1.7452X109

39、28.7262X1081.360.2937pH1.5229X10935.0763X1097.910.0036時間×pH3.1921X10965.3201X1080.830.5696誤 差7.7030X109126.4191X108總 和2.7869X1010239.11 布氏輪藻（Chara braunii Gm.）的托葉長度與生態(tài)環(huán)境的狀況有密切關(guān)系。實驗選擇4種藥物（A：Cd2+，B：Hg2+，C：Cr6+，D：敵枯雙 ），每種藥物（因素）選擇4個水平，兩次重復。加藥培養(yǎng)5個月后，托葉的長度（mm）如下60：因 素水 平1234A6242603001 2076377021 300

40、722B1904007801 0922505008582 210C6506004101 0401 0404809101 300D7809476503008321 2481 8201 300這是一個有重復的兩因素交叉分組實驗設(shè)計，對上述數(shù)據(jù)進行分析，并解釋為什么會得到這樣的結(jié)果？ 答：結(jié)果如下：The SAS SystemAnalysis of Variance ProcedureClass Level InformationClass Levels ValuesDRUG 4 1 2 3 4LEVEL 4 1 2 3 4Number of observations in data set = 3

41、2The SAS SystemAnalysis of Variance ProcedureDependent Variable: LENGTHSum of MeanSource DF Squares Square F Value Pr > FModel 15 3672745.22 244849.68 1.39 0.2613Error 16 2823335.50 176458.47Corrected Total 31 6496080.72R-Square C.V. Root MSE LENGTH Mean0.565379 51.03545 420.070 823.094Source DF

42、Anova SS Mean Square F Value Pr > FDRUG 3 310025.34 103341.78 0.59 0.6331LEVEL 3 1435384.09 478461.36 2.71 0.0796DRUG*LEVEL 9 1927335.78 214148.42 1.21 0.3519根據(jù)以往的經(jīng)驗，重金屬和農(nóng)藥對植物的生長應(yīng)當有影響。然而，實驗結(jié)果卻是藥物、水平及藥物×水平三個因素都是不顯著因素。造成這種結(jié)果的原因是實驗的誤差平方和過大。我們知道，誤差平方和是重復間的平方和。在原始數(shù)據(jù)中，有些重復的數(shù)據(jù)相差甚大，例如，A3和D4的兩次重復間竟然相

43、差4倍有余，相差23倍的也有不少。重復間存在如此之大的偏差，說明實驗材料、實驗環(huán)境（條件）或?qū)嶒灢僮鞔嬖诓灰恢滦浴Ｖ貜烷g過大的偏差，造成過大的誤差均方，使本來存在的效應(yīng)被誤差掩蓋，而不能被檢驗出來。在設(shè)計實驗時，除所研究的因素外，一定要保證各方面的均一性。這一點在設(shè)計實驗和完成實驗的過程中是至關(guān)重要的，一定要特別注意。人們在接受誤差很大的背景下所得到的結(jié)論時，會持保留態(tài)度的。以上數(shù)據(jù)可以歸納成下表：變差來源平方和自由度均方FP藥 物310 025.343103 341.780.590.633 1水 平1 435 384.093478 461.362.710.079 6藥物×水平1 9

44、27 335.789214 148.421.210.351 9誤 差2 823 335.5016176 458.47總 和6 496 080.72319.12 六味木香袋泡劑是一種中藥新劑型。藥物的浸出率與粒度的大小、浸泡時間、浸泡水溫等因素有關(guān)。以下數(shù)據(jù)是不同粒度及不同水量的浸出率（）61，對這些數(shù)據(jù)進行分析，推斷因素的顯著性。粒 度/目10203040加水量/ mL10041.8339.1034.9334.8840.1438.2135.7932.6615041.1840.3035.3334.4037.8938.0834.6832.0520033.8532.9027.2334.2735.27

45、31.0331.0031.36答：這里的因變量是浸出率，它不是二項分布數(shù)據(jù)，不需做變換。結(jié)果如下：The SAS SystemAnalysis of Variance ProcedureClass Level InformationClass Levels ValuesWATER 3 1 2 3GRANULE 4 1 2 3 4Number of observations in data set = 24The SAS SystemAnalysis of Variance ProcedureDependent Variable: PERCENTSum of MeanSource DF Squa

46、res Square F Value Pr > FModel 11 278.901933 25.354721 10.28 0.0002Error 12 29.604000 2.467000Corrected Total 23 308.505933R-Square C.V. Root MSE PERCENT Mean0.904041 4.443402 1.57067 35.3483Source DF Anova SS Mean Square F Value Pr > FWATER 2 126.373908 63.186954 25.61 0.0001GRANULE 3 118.520

47、600 39.506867 16.01 0.0002WATER*GRANULE 6 34.007425 5.667904 2.30 0.1037結(jié)果指出，加水量和粒度都是極顯著因素，但兩者的交互作用是不顯著的。交互作用不顯著的含義是，不是只有在特定的加水量和特定的粒度下才有最佳的浸出率。以上結(jié)果可以歸納成下表：變差來源平方和自由度均方FP加水量126.373 908263.186 95425.610.000 1粒 度118.520 600339.506 86716.010.000 2加水量×粒度34.007 42565.667 9042.300.103 7誤 差29.604 0001

48、22.467 000總 和308.505 933239.13 長沙市2005年7月份不同地點、不同日期和每一天三個時間的空氣溫度測量結(jié)果列在下表中62：地點日期/日溫度/8:0014:0020:00長沙汽車西站429.434.132.5629.735.932.91630.436.333.71831.536.033.22729.236.033.42930.536.633.5岳麓金峰小區(qū)429.039.232.0630.335.832.61630.536.933.81832.137.133.62729.836.833.82931.237.033.7五一中路袁家?guī)X430.537.532.8630.6

49、37.433.11631.038.634.21830.939.134.82729.837.834.22930.638.834.6解放中路瀏城橋430.437.432.4630.537.231.91630.838.534.01830.739.034.62729.537.634.02930.438.534.2馬坡嶺430.235.832.3630.335.932.71630.437.133.71830.737.534.02729.236.533.62930.037.434.2首先判斷這是一個什么模型，然后對上述記錄結(jié)果，做無重復三因素交叉分組方差分析。答：程序和結(jié)果如下：options linesize=76 nodate;data changsha;infile 'e:dataexr9-13e.dat'do place=1 to 5;do date=1 to 6;do time=1 to 3;input temp ;output;end;end;end;run;proc anova;class place date time;model temp=plac