《第十六章-SAS在统计分析中的应用课件》由会员分享，可在线阅读，更多相关《第十六章-SAS在统计分析中的应用课件（111页珍藏版）》请在金锄头文库上搜索。
1、第十六章第十六章 SAS SAS在统计分析中的应用在统计分析中的应用 资料整理与描述资料整理与描述 n例例 用LI6400光合测定仪测定植物叶片速率，共测定6次，分别为12.34，12.82，12.46，12.62，12.52，12.33。求该组数据的算术平均值、几何平均值、调和平均值、中值和极差。 noptions ps=60 ls=80 nodate;ndata photosynthesis;ninput value;ncards;n12.34 12.82 12.46 12.62 12.52 12.33n;nproc univariate;nvar value;noutput;noptio
2、ns ps=60 ls=80 nodate;ndata photosynthesis;ninput value;ncards;n12.34 12.82 12.46 12.62 12.52 12.33n;nproc means var std range cv;nvar value;noutput;频率分布频率分布 n例例 调查100个小区水稻产量的数据如下表（小区计产面积1m2，单位kg）。编制次数分布表并绘制直方图。 3436393634353331383446353933413332344132383842333939303839333834333541313435393039353634
3、3635373536323537362835353633382735373830263637323330333234333437353234323536353535343230363036353836313332333634noptions ps=60 ls=80 nodate;ndata rice;ninput zhugao;ncards;n34363936343533313834n46353933413332344132n38384233393930383933n38343335413134353930n39353634363537353632n35373628353536333827n3
4、5373830263637323330n33323433343735323432n35363535353432303630n36353836313332333634n;nproc freq;ntables zhugao;ntitlefrequency table of zhugao;noptions ps=30;nproc chart;nvbar zhugao;ntitlefrequency chart of zhugao;nrun;n例例 一批种子发芽率为一批种子发芽率为70%70%，每穴播种，每穴播种6 6粒种子计算每穴出粒种子计算每穴出6 6棵苗、棵苗、5 5、4 4、3 3、2 2、1
5、 1、0 0棵苗的概率各为多少？棵苗的概率各为多少？n解：设解：设x x表示每穴出苗数表示每穴出苗数n P(y=6)= C P(y=6)= C6 66 60.70.76 6 0.3 0.30 0=0.117 649=0.117 649n P(y=5)=C P(y=5)=C6 65 50.70.75 5 0.3 0.31 1=0.302 526 =0.302 526 n P(y=4)=C P(y=4)=C6 64 40.70.74 4 0.3 0.32 2=0.324 135=0.324 135n P(y=3)=C P(y=3)=C6 63 30.70.73 3 0.3 0.33 3=0.185
6、 220=0.185 220n P(y=2)=C P(y=2)=C6 62 20.70.72 2 0.3 0.34 4=0.059 535 =0.059 535 n P(y=1)=C P(y=1)=C6 61 10.70.71 1 0.3 0.35 5=0.010 206=0.010 206n P(y=0)=C P(y=0)=C6 60 00.70.70 0 0.3 0.36 6=0.000 729 =0.000 729 data germination；p6=probbnml(0.7,6,6)-probbnml(0.7,6,5); p5=probbnml(0.7,6,5)-probbnml(
7、0.7,6,4); p4=probbnml(0.7,6,4)-probbnml(0.7,6,3); p3=probbnml(0.7,6,3)-probbnml(0.7,6,2); p2=probbnml(0.7,6,2)-probbnml(0.7,6,1); p1=probbnml(0.7,6,1)-probbnml(0.7,6,0); p0=probbnml(0.7,6,0);put p6-p0;proc print;run; 例例 麦田内，平均每麦田内，平均每10m10m2 2有有1 1株杂草，现在要问每株杂草，现在要问每100m100m2 2麦田中，有麦田中，有0 0株杂草，株杂草，1
8、1株杂草，株杂草，2 2株杂草，株杂草，.的概率是多少？的概率是多少？ m=100/10=10 m=100/10=10株株 p p（y y）=10=10y ye e-10-10/y/y！ 将将y=0y=0，1 1，2 2，代入公式即可得出相应概率，代入公式即可得出相应概率data weed;p0=poisson(10,0);p1=poisson(10,1)-poisson(10,0);p2=poisson(10,2)-poisson(10,1);p3=poisson(10,3)-poisson(10,2);p4=poisson(10,4)-poisson(10,3);put p0-p4;pro
9、c print;run;先将先将y y转换为转换为u u ，转化的公式为：，转化的公式为：例例例例 假定假定假定假定y y y y是一随机变数，具有正态分布，平均数是一随机变数，具有正态分布，平均数是一随机变数，具有正态分布，平均数是一随机变数，具有正态分布，平均数 =30=30=30=30，标准差，标准差，标准差，标准差 =5=5=5=5，试计算小于，试计算小于，试计算小于，试计算小于26262626，小于，小于，小于，小于40404040的概率，的概率，的概率，的概率，或者介于或者介于或者介于或者介于26262626和和和和40404040之间的概率以及大于之间的概率以及大于之间的概率以及
10、大于之间的概率以及大于40404040的概率。的概率。的概率。的概率。查附表查附表2 2，当，当 =-0.8=-0.8时，时，F FN N（2626）=0.2119=0.2119，说明这一分，说明这一分布从布从- - 到到2626范围内的变量占全部变量数的范围内的变量占全部变量数的21.19%21.19%，或者说，或者说y y 2626的概率为的概率为0.21190.2119。同样计算同样计算（4040）查附表查附表2 2，当，当 =+2=+2时，时， F FN N（4040）=0.9773=0.9773，这指出从，这指出从- - 到到4040范围内的变量数占全部变量数的范围内的变量数占全部变
11、量数的97.73%97.73%，或者说，或者说，y y 4040的概率为的概率为0.97730.9773。计算计算P（26F(F=5.91F.01.01(4.72)(4.72)，其中在，其中在.05.05和和.01.01水平水平上第上第1 1块田与第块田与第3 3、4 4田的虫口密度有显著差异，田的虫口密度有显著差异，其他处理间差异均不显著。其他处理间差异均不显著。data Cnaphalocrocis;do field=1 to 4;input N;do repeat=1 to N;input pest_density;output;end;end;cards;712 13 14 15 15
12、 16 17614 10 11 13 14 1189 2 10 11 12 13 12 11712 11 10 9 8 10 12;proc glm;class field;model pest_density=field;means field/lSD duncan snk; run;n例 在温室内以4种培养液(l=4)培养某种作物，每种3盆(m=3)，每盆4株(n=4)，一个月后测定其株高生长量(mm),结果如下表，试作方差分析。培养液培养液A AB BC CD D总和总和盆号盆号A A1 1A A2 2A A3 3B B1 1B B2 2B B3 3C C1 1C C2 2C C3 3D
13、 D1 1D D2 2D D3 3生生长长量量505035354545505055555555858565657070606060606565555535354040454560604545606070707070555585856565404030304040505050506565909080807070353545458585353540405050454550505555858565657070707075757575盆总和盆总和T Tijij180180140140175175190190215215220220320320280280280280220220265265290290
14、T=2725T=2725培养液培养液总和总和T Ti i495495625625880880775775培养液培养液平均平均41.341.352.152.173.373.364.664.6n一、自由度和平方和分解一、自由度和平方和分解n总自由度总自由度DFDFT T= =lmnlmn-1=(434)-1=47-1=(434)-1=47n培养液间自由度培养液间自由度DFDFt t= =l l-1=4-1=3-1=4-1=3n培养液内盆间自由度培养液内盆间自由度nDFDFe1e1= =l l( (m m-1)=4(3-1)=8-1)=4(3-1)=8n盆内株间自由度盆内株间自由度nDFDFe2e2
15、= =lmlm( (n n-1)=43(4-1)=36-1)=43(4-1)=36n2 2、培养液间差异、培养液间差异二、二、F F 测验测验1 1、盆间差异、盆间差异假设假设H H0 0：，求得：求得：F F=157.81/89.06=1.77=157.81/89.06=1.77此此F F值小于值小于1 1=8,=8,2 2=36 =36 F F0.050.05=2.22,=2.22,所以接受所以接受H H0 0假设假设，求得：，求得：n推断：该试验同一培养液内盆间的生长量无推断：该试验同一培养液内盆间的生长量无显著差异；而不同培养液间的生长量有极显显著差异；而不同培养液间的生长量有极显著的
16、差异。著的差异。 F F=2375.25/157.81=15.05=2375.25/157.81=15.05此此F F值大于值大于1 1=3,=3,2 2=8=8F F0.010.01=7.59,=7.59,故否定故否定，接受接受变异来源变异来源DFDFSSSSMSMSF FF F0.050.05F F0.010.01培养液间培养液间3 37126.567126.562375.522375.5215.0515.054.074.077.597.59培养液内盆间培养液内盆间8 81262.501262.50157.81157.811.771.772.222.223.043.04盆内株间盆内株间36
17、363206.253206.2589.0689.06总变异总变异474711595.3111595.31方差分析表方差分析表n三、各培养液平均数间的比较三、各培养液平均数间的比较p pSSRSSR0.050.05SSRSSR0.010.01LSRLSR0.050.05LSRLSR0.010.012 23.263.264.744.7411.8311.8317.2117.213 33.393.395.005.0012.3112.3118.1518.154 43.473.475.145.1412.6012.6018.6618.664 4种培养液的种培养液的LSRLSR值值( (新复极差测验新复极差测
18、验) )n4 4种培养液植株生长量种培养液植株生长量(mm)(mm)的差异显著性的差异显著性培养液平均生长量差异显著性0.050.01C73.3aAD64.6a ABB52.1 b BC A41.3 b Cdata growth;do i=1 to 4; input solution$; do j=1 to 3; input pot$; do k=1 to 4; input height; output; end; end;end;cards;A A1 50 55 40 35 A2 35 35 30 40 A3 45 40 40 50B B1 50 45 50 45 B2 55 60 50 5
19、0 B3 55 45 65 55C C1 85 60 90 85 C2 65 70 80 65 C3 70 70 70 70D D1 60 55 35 70 D2 60 85 45 75 D3 65 65 85 75;proc glm;class solution pot;model height=solution pot;means solution pot/duncan; run;例例 采用采用5 5种生长素处理豌豆，未处理为对照，待种种生长素处理豌豆，未处理为对照，待种子发芽后，分别每盆中移植子发芽后，分别每盆中移植4 4株，每组株，每组6 6盆，每盆一盆，每盆一个处理，试验共有个处理，
20、试验共有4 4组组2424盆，并按组排列于温室中，盆，并按组排列于温室中，使同组各盆的环境条件一致。当各盆见第一朵花时使同组各盆的环境条件一致。当各盆见第一朵花时记录记录4 4株豌豆的总节间数，结果见下表，试作方差分株豌豆的总节间数，结果见下表，试作方差分析。析。处理（处理（A A）组（组（B B）总和总和T Ti.i.平均平均I IIIIIIIIIIIIVIV未处理未处理(CK)(CK)606062626161606024324360.860.8赤霉素赤霉素656565656868656526326365.865.8动力精动力精636361616161606024524561.361.3吲哚
21、乙酸吲哚乙酸硫酸腺嘌呤硫酸腺嘌呤马来酸马来酸64646262616167676565626263636262626261616464656525525525325325025063.863.863.363.362.562.5总和总和T.T.j j375375382382377377375375T T15091509(1)(1)自由度和平方和的分解自由度和平方和的分解变异来源变异来源DFDFSSSSMSMSF FF F0.050.05F F0.010.01组间组间3 35.455.451.821.821 1处理间处理间5 565.8765.8713.1713.174.564.562.902.90
22、4.564.56误差误差151543.3043.302.892.89总变异总变异2323114.62114.62方差分析表方差分析表查得查得1515时，时，t t0.050.05=2.131=2.131，t t0.010.01=2.947=2.947LSDLSD0.050.05=1.2022.131=2.56, =1.2022.131=2.56, LSDLSD0.010.01=1.202 2.947=3.54=1.202 2.947=3.54处理处理平均数平均数与对照的差数与对照的差数对照（对照（CKCK）赤霉素赤霉素动力精动力精吲哚乙酸吲哚乙酸硫酸腺嘌呤硫酸腺嘌呤马来酸马来酸60.860.8
23、65.865.861.361.363.863.863.363.362.562.55.0*5.0*0.50.53.0*3.0*2.52.51.71.7处理间比较处理间比较此例有预先指定的对照，故用此例有预先指定的对照，故用LSD LSD 法。法。data pea;do hormone=1 to 6; do group=1 to 4; input node; output; end;end;cards;606261606565686563616160646763616265626461626265;proc glm;class hormone group;model node=hormone gr
24、oup;means hormone group/lsd; run;n例 施用A1、A2、A33种肥料于B1、B2、B33种土壤，以小麦为指示作物，每处理组合种3盆，得产量结果(g)如下表，试作方差分析。肥料种类肥料种类(A)(A)盆盆土壤种类土壤种类(B)(B)总和总和T Ti.i.平均平均B B1 1( (油沙油沙) )B B2 2( (二合二合) )B B3 3( (白僵白僵) )A1A11 121.421.419.619.617.617.6169.2169.218.818.82 221.221.218.818.816.616.63 320.120.116.416.417.517.5T T
25、ij.ij.62.762.754.854.851.751.7A2A21 112.012.013.013.013.313.3118.2118.213.113.12 214.214.213.713.714.014.03 312.112.112.012.013.913.9T Tij.ij.38.338.338.738.741.241.2A3A31 112.812.814.214.212.012.0122.0122.013.613.62 213.813.813.613.614.614.63 313.713.713.313.314.014.0T Tij.ij.40.340.341.141.140.640
26、.6总和总和T T.j.j.141.3141.3134.6134.6133.5133.5T T=409.4=409.4平均平均15.715.715.015.014.814.8n1 1、自由度和平方和的分解、自由度和平方和的分解2 2、F F 测验测验将上述结果及自由度录于方差分析表中，将上述结果及自由度录于方差分析表中，以固定模型作以固定模型作F F测验测验变异来源变异来源DFDFSSSSMSMSF FF F0.010.01处理处理8 8202.58202.5825.2225.2227.2827.283.713.71肥料间肥料间2 2179.38179.3889.6989.6996.6596.
27、656.016.01土类间土类间2 23.963.961.981.982.132.136.016.01 肥料肥料土类土类4 419.2419.244.814.815.185.184.584.58试验误差试验误差181816.7016.700.9280.928总变异总变异2626219.28219.28n3 3、平均数的比较、平均数的比较n(1)(1)各处理组合平均数的比较各处理组合平均数的比较p p2 23 34 45 56 67 78 89 9SSRSSR0.050.052.972.973.123.123.213.213.273.273.323.323.353.353.373.373.393
28、.39SSRSSR0.010.014.074.074.274.274.384.384.464.464.534.534.594.594.644.644.684.68LSRLSR0.050.051.651.651.731.731.781.781.821.821.851.851.861.861.871.871.881.88LSRLSR0.010.012.262.262.372.372.442.442.482.482.522.522.552.552.582.582.602.60各处理组合平均数的各处理组合平均数的LSRLSR值值各处理平均数的新复极差测验各处理平均数的新复极差测验处理组合平均数(g)差
29、异显著性0.050.01A1B120.9aAA1B218.3 b BA1B317.2 b BA2B313.7 c CA3B213.7 c CA3B313.5 c CA3B113.4 c CA2B212.9 c CA2B112.8 c C (2) (2)各肥类平均数的比较各肥类平均数的比较p pSSRSSR 0.050.05SSRSSR 0.010.01LSRLSR 0.050.05LSRLSR 0.010.012 22.972.974.074.070.950.951.301.303 33.123.124.274.271.001.001.371.37肥料种类肥料种类平均数平均数差异显著性差异显著
30、性0.050.050.010.01A A1 118.818.8a aA AA A3 313.613.6 b b B BA A2 213.113.1 b b B B肥类平均数的肥类平均数的LSRLSR测验测验data wheat;do fertilizer=1 to 3; do soil=1 to 3; do repeat=1 to 3; input yield; output; end; end;end;cards;21.4 21.2 20.119.6 18.8 16.417.6 16.6 17.512.0 14.2 12.113.0 13.7 12.013.3 14.0 13.912.8 1
31、3.8 13.714.2 13.6 13.312.0 14.6 14.0;proc glm;class fertilizer soil;model yield=fertilizer soil fertilizer*soil;means fertilizer soil/duncan; run;例例 调查经过种子灭菌处理与未经种子灭菌处理的小麦调查经过种子灭菌处理与未经种子灭菌处理的小麦发生散黑穗病的穗数，得相依表如下，试分析种子灭发生散黑穗病的穗数，得相依表如下，试分析种子灭菌与否和散黑穗病穗多少是否有关。菌与否和散黑穗病穗多少是否有关。处理项目处理项目发病穗数发病穗数未发病穗数未发病穗数总数总
32、数种子灭菌种子灭菌26(34.7)26(34.7)50(41.3)50(41.3)7676种子未灭菌种子未灭菌184(175.3)184(175.3)200(208.7)200(208.7)384384总数总数210210250250460460查附表查附表6 6， 现实得现实得 故故P0.05Pt 0.01，故接受HA。data Scirpophaga;input tem time;cards;35.5 12 34.1 16 31.7 9 40.3 2 36.8 7 40.2 3 31.7 13 39.2 9 44.2 -1;proc corr nosimple;var tem time;r
33、un;data wheat;input x1 x2 y;cards;18.5132.066818.3132.064017.2127.859719.9110.358119.8114.560718.8122.561822.3114.766321.2 98.155117.9115.454418.9129.064120.2102.452120.8108.157321.4113.760821.3110.958318.0128.1576;proc reg;model y=x1 x2;run;data wheat;input x1 x2 x3 x4 y;cards;30.8 33.0 50.0 90 520
34、.823.6 33.6 28.0 64 195.031.5 34.0 36.6 82 424.019.8 32.0 36.0 70 213.527.7 26.0 47.2 74 403.327.7 39.0 41.8 83 461.716.2 43.7 44.1 83 248.031.2 33.7 47.5 80 410.023.9 34.0 45.3 75 378.330.3 38.9 36.5 78 400.835.0 32.5 36.0 90 395.033.3 37.2 35.9 85 400.027.0 32.8 35.4 70 267.525.2 36.2 42.9 70 361.
35、323.6 34.0 33.5 82 233.821.3 32.9 38.6 80 210.021.1 42.0 23.1 81 168.319.6 50.0 40.3 77 400.021.6 45.1 39.3 80 319.432.3 25.6 39.8 71 376.2;proc stepwise;model y=x1-x4;run;F 测验 ：查表法 ：data wheat;input x1 x2 y;cards;18.5132.066818.3132.064017.2127.859719.9110.358119.8114.560718.8122.561822.3114.766321
36、.2 98.155117.9115.454418.9129.064120.2102.452120.8108.157321.4113.760821.3110.958318.0128.1576;proc cancorr short;var x1 x2;with y;run;t 测验 ：查表法 ：data wheat;input x1 x2 y;cards;18.5 132.0 668 18.3 132.0 64017.2 127.8 597 19.9 110.3 58119.8 114.5 607 18.8 122.5 61822.3 114.7 663 21.2 98.1 55117.9 115
37、.4 544 18.9 129.0 64120.2 102.4 521 20.8 108.1 57321.4 113.7 608 21.3 110.9 58318.0 128.1 576;proc corr nosimple;var x1 y;partial x2;title Correlation Matrix for x1 and y;proc corr nosimple;var x2 y;partial x1;title Correlation Matrix for x2 and y;proc corr nosimple;var x1 x2;partial y;title Correla
38、tion Matrix for x1 and x2;run;例例例例 为研究为研究为研究为研究4 4种不同肥料种不同肥料种不同肥料种不同肥料A1A1、A2A2、A3A3、A4A4对梨树单株产对梨树单株产对梨树单株产对梨树单株产量的影响，选择量的影响，选择量的影响，选择量的影响，选择4040株梨树做试验，把株梨树做试验，把株梨树做试验，把株梨树做试验，把4040株梨树完全随株梨树完全随株梨树完全随株梨树完全随机分为机分为机分为机分为4 4组，每组包含组，每组包含组，每组包含组，每组包含1010株梨树，每组施用株梨树，每组施用株梨树，每组施用株梨树，每组施用1 1种肥料。种肥料。种肥料。种肥料。各
39、株梨树的起始干周和单株产量列表各株梨树的起始干周和单株产量列表各株梨树的起始干周和单株产量列表各株梨树的起始干周和单株产量列表10-210-2，检验，检验，检验，检验4 4种肥种肥种肥种肥料的单株产量是否有差异。料的单株产量是否有差异。料的单株产量是否有差异。料的单株产量是否有差异。 表表表表10-2 10-2 梨树梨树梨树梨树4 4种肥料比较试验的起始干周（种肥料比较试验的起始干周（种肥料比较试验的起始干周（种肥料比较试验的起始干周（x, cmx, cm）与单株产量（与单株产量（与单株产量（与单株产量（y, kgy, kg）变异来源变异来源SSSSdfdfMSMSF FF F0.010.01
40、误差回归误差回归475.993475.9931 1475.993475.99311.29 *11.29 *7.427.42误差离回归误差离回归1475.0071475.007353542.14342.143误差总和误差总和195119513636表表表表10-5 10-5 单株产量与起始干周的回归关系显著性检验表单株产量与起始干周的回归关系显著性检验表单株产量与起始干周的回归关系显著性检验表单株产量与起始干周的回归关系显著性检验表表表10-6 矫正梨树单株产量的方差分析表矫正梨树单株产量的方差分析表 表表表表10-7 10-7 各处理的校正平均单株产量计算表各处理的校正平均单株产量计算表各处理
41、的校正平均单株产量计算表各处理的校正平均单株产量计算表 肥料肥料实际平均单实际平均单株产量株产量校正平均单校正平均单株产量株产量A1A124.6-23.67524.6-23.6750.7359*0.7359*（ 24.6- 24.6-23.67523.675）75.575.574.81974.819A2A222.9-23.67522.9-23.6750.7359*0.7359*（22.9-22.9-23.67523.675）67.467.467.9767.97A3A323.2-23.67523.2-23.6750.7359*0.7359*（23.2-23.2-23.67523.675）58.8
42、58.859.1559.15A4A424.0-23.67524.0-23.6750.7359*0.7359*（ 24.0- 24.0-23.67523.675）54.554.554.26154.261肥料肥料校正平均单校正平均单株产量株产量A1A174.81974.81920.558*20.558*15.668*15.668*6.849*6.849*A2A267.9767.9713.709*13.709*8.82*8.82*A3A359.1559.154.8894.889A4A454.26154.261表表表表10-10 10-10 各处理校正平均单株产量多重比较表（各处理校正平均单株产量多重
43、比较表（各处理校正平均单株产量多重比较表（各处理校正平均单株产量多重比较表（SSRSSR法）法）法）法） data pear;do i=1 to 4; do j=1 to 10; input x y; output; end;end;cards;36 89 30 80 26 74 23 80 26 85 30 68 20 73 19 68 20 80 16 5828 64 27 81 27 73 24 67 25 77 23 67 20 64 18 65 17 59 20 5728 55 33 62 26 58 22 58 23 66 20 55 22 60 23 71 18 55 17 48
44、32 52 23 58 27 64 23 62 27 54 28 54 20 55 24 44 19 51 17 51;proc glm;class i;model y=x i/solution;lsmeans i/stderr pdiff;run;data bean;do i=1 to 6; do j=1 to 4; input x y; output; end;end;cards;34.0 93.0 33.4 94.8 34.7 91.7 38.9 80.839.6 47.3 39.8 51.5 51.2 33.3 52.0 27.231.7 81.4 30.1 109 33.8 71.6
45、 39.6 57.537.7 66.9 38.2 74.1 40.3 64.7 39.4 69.324.9 119.5 24.0 128.5 24.9 125.6 23.5 129.030.3 106.6 29.1 111.4 31.7 99.0 28.3 126.1;proc glm;class i j;model y=x i j/solution;lsmeans i/stderr pdiff;run;例例 水稻品比水稻品比5X5拉丁方试验的产量结果拉丁方试验的产量结果横行区组横行区组 纵行区组纵行区组 Tr D(37)A(38)C(38) B(44) E(38) 195 D(37)A(38
46、)C(38) B(44) E(38) 195 B(48)E(40)D(36) C(32) A(35) 191 B(48)E(40)D(36) C(32) A(35) 191 C(27)B(32)A(32) E(30) D(26) 147 C(27)B(32)A(32) E(30) D(26) 147 E(28)D(37)B(43) A(38) C(41) 187 E(28)D(37)B(43) A(38) C(41) 187 A(34)C(30)E(27) D(30) B(41) 162 A(34)C(30)E(27) D(30) B(41) 162Tc 174 177 176 174 181
47、 T=882Tc 174 177 176 174 181 T=882 方差分析表方差分析表变异来源变异来源 DF SS MS F F0.05横行区组横行区组 4 348.64 87.16 纵行区组纵行区组 4 6.64 1.66 品种品种 4 271.44 67.86 4.33 * 3.26 试验误差试验误差 12 188.32 15.69 总变异总变异 24 815.04 LSR法法水稻品比试验的新复极差测验水稻品比试验的新复极差测验品种品种 平均产量平均产量 差异显著性差异显著性 5% 1%B 41.6 a AA 35.4 b ABC 33.6 b ABD 33.2 bE 32.6 b B
48、data rice;input row col entry$ yield;cards;1 1 D 37 1 2 A 38 1 3 C 38 1 4 B 44 1 5 E 382 1 B 48 2 2 E 40 2 3 D 36 2 4 C 32 2 5 A 353 1 C 27 3 2 B 32 3 3 A 32 3 4 E 30 3 5 D 264 1 E 28 4 2 D 37 4 3 B 43 4 4 A 38 4 5 C 415 1 A 34 5 2 C 30 5 3 E 27 5 4 D 30 5 5 B 41;proc anova;class row col entry;model
49、 yield=row col entry;means entry/duncan;run;A1B1 A2B2 A3B3 A3B2 A1B3 A3B1 A1B2 A2B1 8 7 10 8 6 7 7 9A2B3 A3B2 A1B2 A3B1 A1B3 A2B A2B2 A3B3 A1B17 7 7 7 5 9 9 9 8A3B1 A1B3 A2B1 A1B2 A2B2 A3B3 A1B1 A2B3 A3B2 6 6 8 6 6 9 8 6 8例例 区组和处理两向表处理 总和TABA1B1 8 8 8 24A1B2 7 7 6 20A1B3 6 5 6 17A2B1 9 9 8 26A2B2 7
50、9 6 22A2B3 8 7 6 21A3B1 7 7 6 20A3B2 8 7 8 23A3B3 10 9 9 28 总和Tr 70 68 63 T=201变异来源变异来源F0.05区组区组处理处理282.8930.001.453.722.967.65*3.632.59品种品种26.233.126.33*3.63密度密度21.560.781.593.63品种品种 密度密度422.215.5511.33*3.01误差误差167.780.49总变异总变异2640.67二因素随机区组设计早稻栽培实验方差分析表二因素随机区组设计早稻栽培实验方差分析表 （1）品种间的比较品种间的比较 SE=0.233
51、 新复极差测验新复极差测验品种 产量 差异显著性 5% 1%A3 7.9 a AA2 7.7 a ABA1 6.8 b Bdata rice;do A=1 to 3; do B=1 to 3; do block=1 to 3; input yield; output; end; end;end;cards;8 8 8 7 7 6 6 5 6 9 9 8 7 9 6 8 7 6 7 7 6 8 7 8 10 9 9;proc anova;class A B block;model yield=A B block A*B;means A B/duncan;run;例例 用三种不同方法从植物中提取有
52、效成分，按用三种不同方法从植物中提取有效成分，按4种不同浓度种不同浓度添加到培养基中，观察对培养基植株生长的促进作用，记录添加到培养基中，观察对培养基植株生长的促进作用，记录培养一个月后植株重量。（裂区设计培养一个月后植株重量。（裂区设计A因素为主处理，因素为主处理，B因素因素为副处理）为副处理） 重复重复 重复重复 A1 A3 A2 A3 A2 A1B237B129B315B231B413B313B318B417B416B130B128B231B127B314B412B313B232B314B415B228B228B129B416B128B415B317B231B413B125B229B23
53、1B132B126B311B310B412重复重复重复重复A1 A3 A2A1 A3 A21、结果整理 区组和处理两向表主处理A 副处理B 区组 TAB TA A1 B1 29 28 32 89 B2 37 32 31 100 B3 18 14 17 49 B4 17 16 15 48 Tm 101 90 95 286A2 B1 28 29 25 82 B2 31 28 29 88 B3 13 13 10 36 B4 13 12 12 37 Tm 85 82 76 243A3 B1 30 27 26 83 B2 31 28 31 90 B3 15 14 11 40 B4 16 15 13 44
54、 Tm 92 84 81 257Tr 278 256 252 T=786裂区方差分析表变异来源变异来源 DF SS MS F F0.05主区部分主区部分 区组区组 2 32.67 16.34 7.14* 6.94 A 2 80.17 40.09 17.51* 6.94 Ea 4 9.16 2.29 总变异总变异 8 122 副区部分副区部分 B 3 2179.67 726.56 282.71* 3.16 AxB 6 7.16 1.19 1 Eb 18 46.17 2.57总变异总变异 35 2355data ingredient;do A=1 to 3; do B=1 to 4; do block=1 to 3; input height; output; end; end;end;cards;29 28 32 37 32 31 18 14 17 17 16 15 28 29 25 31 28 29 13 13 10 13 12 1230 27 26 31 28 31 15 14 11 16 15 13;proc anova;class A B block;model height=A B block A*B A*block;test h=A e=A*block;means A A*B/duncan;run;
