SAS作业.doc

SAS作业 2011级栽培生理班 2011301110033 唐红玲SAS课后习题作业院系：植物科学技术学院学号：2011301110033姓名：唐红玲二O一一年十二月习题2.1.1观测200头某品种仔猪的一月窝重（单位：kg）得到频数表如下：组下限组上限组中值频数7.523.515.51023.539.531.51939.555.547.53055.571.563.56171.587.579.54987.5103.595.524103.5119.5111.57试用SAS作样本观测值整理后的描述性统计，并画出样本观测值的频率分布直方图。解：data ex; input x f ;cards;15.5 10 31.5 19 47.5 30 63.5 61 79.5 49 95.5 24 111.5 7;proc univariate; var x;freq f; run;Variable: xQuantiles (Definition 5)data hist01;input x ;cards;15.5 15.5 15.5 15.5 15.5 15.5 15.5 15.5 15.5 15.531.5 31.5 31.5 31.5 31.5 31.5 31.5 31.5 31.5 31.531.5 31.5 31.5 31.5 31.5 31.5 31.5 31.5 31.547.5 47.5 47.5 47.5 47.5 47.5 47.5 47.5 47.5 47.547.5 47.5 47.5 47.5 47.5 47.5 47.5 47.5 47.5 47.547.5 47.5 47.5 47.5 47.5 47.5 47.5 47.5 47.5 47.563.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.563.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.563.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.563.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.563.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.563.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.5 63.563.579.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 79.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.5 95.595.5 95.5 95.5 95.5 111.5 111.5 111.5 111.5 111.5 111.5 111.5 ;proc gchart;vbar x/type=pct space=0;run;习题2.4.1已知我国14岁女学生的平均体重为43.38kg，从该年龄的女学生中抽查10名运动员的体重，分别为39,36,43,43,40,46,45,45,42,41kg。试用SAS检验这些运动员的体重与上述平均体重是否有显著的差异（=0.05）解：假设H0：=43.38，运动员体重与平均体重无显著差异HA：43.38，运动员体重与平均体重有显著差异data ex; input x ; y= x-43.38;cards;39 36 43 43 40 46 45 45 42 41;proc means mean std t prt; var y ; run;SAS结果：因为P|X|T|0.05,所以接受H0：=43.38，即运动员体重与平均体重无显著差异.习题2.4.2某小麦品种经过4代选育，从第5代和第6代中分别抽出10株得到它们株高的观测值分别为66,65,66,68,62,65,63,66,68,62以及64,61,57,65,65,63,62,63,64,60.试用SAS检验这两代的株高是否有显著的差异（=0.05）解：data ex; do a = 1 to 2; do i = 1 to 10;input x ;output;end;end;cards;66 65 66 68 62 65 63 66 68 6264 61 57 65 65 63 62 63 64 60;proc ttest cochran;class a;var x;proc print;run; A N Mean Std Dev Std Error1 10 65.1 2.1833 0.69042 10 62.4 2.5033 0.7916Variable: x假设H0：第5代和第6代的株高没有显著差异，即H0：1=2,，对HA：12,显著水平=0.05如果认为方差相等，则DF=18，Prob|T|为0.0193；0.010.0193|T|为0.0194；0.010.0194|T|为0.0302；0.010.03020.05二者结果均表示：放弃H0综上可得：这两代的株高存在显著差异。习题2.5.1测定4种种植密度下金皇后玉米的千粒重（单位：g）如下：种植密度千粒重1234247238214210258244227204256246221200251236218210试用SAS对4种种植密度下金皇后玉米千粒重的观测值作正态性检验。解：假设H0：4种种植密度下金皇后玉米千粒重的观测值服从正态分布；H1：4种种植密度下金皇后玉米千粒重的观测值不服从正态分布data ex; do a = 1 to 4; do i= 1 to 4;input x ; output; end; end;cards ;247 258 256 251 238 244 246 236214 227 221 218 210 204 200 210;proc univariate normal ; run ;W检验的临界值0.05=0.859，PW0.05=0.859=0.05SAS结果表明PW0.867632=0.0250 FModel 2 934.722222 467.361111 31.10 0.0007A 2 934.722222 467.361111 31.10 0.0007Error 6 90.166667 15.027778C total 8 1024.888889SAS结果表明0.0007 |t|为0.0204 |t|为0.05160.05接受原假设，认为方差是同质的，即变换后的方差是齐性的。 Dunnetts t Tests for y根据上述结果，可知处理3与1有显著差异，处理2与1没有显著差异。各均值差的95%的区间估计为：Mu3-mu1 (0.00825，0.44087)Mu2-mu1 (-0.16223，0.27039)习题3.3.2有4个学生分析5个品种的稻米，得到含N量（单位：mg）的测量值如下：品种学生ABCDE甲2.42.53.23.42.0乙2.62.23.23.51.8丙2.12.73.53.81.8丁2.42.73.13.22.3试用SAS作方差分析，估计各个总体的未知参数i。和.j及。若有必要，试求出两两总体均值差的双侧0.95置信区间。解：设H01为各i。相等，也就是各学生分析的5个品种稻米的含氮量之间没有显著的差异，H02为.j相等，也就是品种之间的含氮量没有显著的差异。data ex; do a=1 to 4; do b= 1 to 5;input x ;output; end; end;cards;2.4 2.5 3.2 3.4 2.0 2.6 2.2 3.2 3.5 1.82.1 2.7 3.5 3.8 1.8 2.4 2.7 3.1 3.2 2.3;proc anova; class a b; model x= a b ;means a b/lsd cldiff; run;用SAS计算得到方差分析表为：Dependent Variable: x15因此接受H01，放弃H02，认为i。相等，.j不相等，即各学生分析的5个品种稻米的含氮量之间没有有显著的差异，品种之间的含氮量有显著的差异。t Tests (LSD) for x习题3.4.2某种化工过程在3种浓度、4种温度下成品的得率如下：温度浓度102438522%14,1011,1113,910,124%9,710,87,116,106%5,1113,1412,1314,10试先作得率数据的平方根反正弦变换后再用SAS作方差分析。解：设H01为各i.相等，即各种浓度对应的观测值之间没有显著差异，H02为各。j相等，即各种温度对应的观测值之间没有显著的差异，H03为各ij都等于0，即各种浓度与温度搭配试验的交互作用都等于0.data ex; do a=1 to 3; do b=1 to 4;do i= 1 to 2; input x ; y= arsin(sqrt(x/100);output; end; end; end;cards;14 10 11 11 13 9 10 12 9 7 10 87 11 6 10 5 11 13 14 12 13 14 10;proc anova; class a b; model y= a b a*b;means a b/lsd cldiff;run; 用SAS计算得到方差分析表：Dependent Variable: y结果显示：因素a*b的Pr |t|为0.57950.05，A与B交互不显著所以不考虑因素A与B的交互作用，即将上述程序中model y=a b a*b改为 model y=a b再次运行程序后结果如下：Dependent Variable: y Sum of Source DF Squares Mean Square F Value Pr F Model 5 0.01575797 0.00315159 2.07 0.1165 Error 18 0.02737766 0.00152098 Corrected Total 23 0.04313562 R-Square Coeff Var Root MSE y Mean 0.365312 11.95360 0.039000 0.326259 Source DF Anova SS Mean Square F Value Pr F a 2 0.01211420 0.00605710 3.98 0.0370 b 3 0.00364376 0.00121459 0.80 0.5107结果显示：拒绝H 01，也就是各种浓度(A)下成品的得率之间有显著的差异，接受H 02，也就是各种温度(B)下成品的得率之间没有显著的差异，所以对因素A作多重比较并求均值差的置信区间，即对上述程序增加means a /lsd cldiff; run;再次运行程序后结果如下： t Tests (LSD) for y Alpha 0.05 Error Degrees of Freedom 18 Error Mean Square 0.001521 Critical Value of t 2.10092 Least Significant Difference 0.041Comparisons significant at the 0.05 level are indicated by *. Difference a Between 95% Confidence Comparison Means Limits 3 - 1 0.00167 -0.03930 0.04264 3 - 2 0.04847 0.00750 0.08944 * 1 - 3 -0.00167 -0.04264 0.03930 1 - 2 0.04680 0.00583 0.08777 * 2 - 3 -0.04847 -0.08944 -0.00750 * 2 - 1 -0.04680 -0.08777 -0.00583 *根据上述结果，得出浓度3与2有显著差异，浓度1与2有显著差异，而浓度1与3没有显著差异。各均值差的95%的区间估计为：A3-A2 (0.00750，0.08944) A1-A2 (0.00583，0.08777) A1-A3 (-0.04264，0.03930) 习题3.5.2随机地抽取12个玉米自交系，以其中的3个为父本，9个为母本，每个父本又随机地与9个母本中的3个母本自交系杂交，共配成9个组合，每一个组合在田间种3个小区，共27个小区，收货时观测每穗行数的均值，结果如下表：父本母本每穗行数A12315.6 ， 16.4 ， 15.613.6 ， 15.6 ， 14.812.0 ， 12.0 ， 12.8B45616.0 ， 16.0 ， 15.615.6 ， 15.6 ， 16.814.4 ， 15.6 ， 16.8C78914.8 ， 15.6 ， 14.817.6 ， 18.8 ， 18.014.4 ， 15.2 ， 15.6试用SAS作方差分析，估计各个总体的未知参数i。，ij和。若有必要，试求出两两总体均值差的双侧0.95置信区间。解：设H01为各i。相等，也就是玉米自交系父本与母本杂交后代的每穗行数之间没有显著差异；H02为各ij与i。相等，也就是父本相同与各母本杂交后代的每穗行数之间没有显著差异。data ex;do a=1 to 3; do b=1 to 3;do i=1 to 3;input x ;output;end;end;end;cards;15.6 16.4 15.6 13.6 15.6 14.8 12.0 12.0 12.816.0 16.0 15.6 15.6 15.6 16.8 14.4 15.6 16.814.8 15.6 14.8 17.6 18.8 18.0 14.4 15.2 15.6; proc anova;class a b;model x=a b(a);means a b(a)/Duncan;means a b(a)/lsd cldiff;run;Dependent Variable: x因此放弃H01和H02，认为各i。不相等，各ij与i。也不相等，即玉米自交系父本与母本杂交后代的每穗行数之间有显著差异，父本相同与各母本杂交后代的每穗行数之间有显著差异。Duncans Multiple Range Test for x t Tests (LSD) for x 由以上结果可得如下方差分析表：方差来源平方和自由度均方和F值显著性AB(a)误差 总和16.65 39.887.01 63.552 6 18 268.326.650.3921.3717.06*根据上述结果，得出拒绝H 01，也就是各父本每穗行数之间有显著的差异，拒绝H 02，也就是同父本不同母本每穗行数之间有显著的差异。即父本2与1有显著差异，父本3与1有显著差异，而父本2与3没有显著差异。各均值差的95%的区间估计为：A2-A3 (-0.9960，0.2404) A2-A1 (0.8262，2.0626) A3-A1 (1.2040，2.4404) 习题4.1.2小麦基本苗数x及有效穗数Y(单位：万/亩)的5组观测数据如下: 基本苗数x15.025.830.036.644.4有效穗数y39.441.941.043.149.2试用SAS求线性回归方程并用三种方法作显著性检验，若x26，求Y的0.95预测区间。解：data ex;input x y ; cards;15.0 39.4 25.8 41.9 30.0 41.0 36.6 43.1 44.4 49.2 26. ;proc gplot;plot y*x;/* 以y为纵坐标，以x为横坐标*/symbol i=rl v=dot;/* i=rl表示画回归直线*/* v=dot表示观测值对应的点标记为小圆点*/proc reg;model y=x/cli;run;/*y=x表示以y为因变量，以x为自变量，*/*cli表示要求预测值的95%置信区间*/SAS程序结果表明所求线性回归方程为y=33.76085+0.30168x习题4.3.1三品种试验，每小区的株数x及产量Y(单位：Kg)的观测值如下，试用SAS以x为协变量作协方差分析。 品种 观 测 值(x,y) 1 3,10 2,8 1,8 2,11 2 4,12 3,12 3,10 5,13 3 1,6 2,5 3,8 1,7data ex;do a=1 to 3;do i=1 to 3;input x y ;output;end;end;cards;3 10 2 8 1 8 2 11 4 12 3 123 10 5 13 1 6 2 5 3 8 1 7;pro