概率论及sas作业.docVIP

概率论及SAS作业 P142 1．在校园中有3盏诱蛾灯，连续5天诱蛾记录如下： 地点 观测记录 1 12 46 30 38 17 2 11 21 24 5 22 3 13 4 16 11 8 试用SAS先做Levene的F检验，再作对数转换后的方差分析。 解：（1）检验方差是否同性质，用Levene F检验： H0：各方差具有相同性质 sas代码： data cap_fly ;do a=1 to 3; do i=1 to 5; input x@@;output ;end;end; cards; 12 46 30 38 17 11 21 24 5 22 13 4 16 11 8 ; proc anova;class a ;model x=a; means a/hovtest;run; 程序输出结果： Levenes Test for Homogeneity of x Variance ANOVA of Squared Deviations from Group Means Sum of Mean Source DF Squares Square F Value Pr F a 2 55635.7 27817.9 4.47 0.0353 Error 12 74596.9 6216.4 0.03530.05，拒绝原假设，认为方差不同质； （2）对数转换，再做方差分析 对数转换后，提出H0：各总体方差具有同质性 进行Levene F检验，结果见下。 方差分析：H0：各地点诱蛾数量不存在显著性差异 Sas代码： data cap_fly ;do a=1 to 3; do i=1 to 5; input x@@;y=log(x); output ;end;end; cards; 12 46 30 38 17 11 21 24 5 22 13 4 16 11 8 ; proc anova;class a ;model y=a; means a/hovtest;run; 程序输出结果： 方差来源 平方和 自由度 均方和 F值 显著性 因素A 2.50 2 1.25 3.57 N 误差 4.21 12 0.35 总和 6.71 14 可见转换后方差同质（0.80270.05），满足方差同质性假设，方差分析的结果是接受原假设，认为各地点诱蛾数量不存在显著差异（显著性水平0.05） P147 对3个玉米品种用3种方法进行管理，得到9个小区的产量（单位：kg）如下： 品种 方法 B1 B2 B3 A1 11 30 43 A2 31 29 45 A3 15 25 47 作方差分析，估计各总体的未知参数μi·和μ·j及μ·若有必要，求出两两均值差的双侧0.95置信区间。 解：H01：各μi·都相等，H02:各μ·j都相等 作双因素方差分析，并对均值差作双侧区间估计。 Sas代码： data corn;do a=1 to 3; do b= 1 to 3;input x@@;output;end;end; cards; 11 30 43 31 29 45 15 25 47 ; proc anova