sas第8周方差分析-作业

董大均书p286 第3,4,5题

image 

data test1;
   do i=1 to 8;
      do season=1 to 4;
         input x @@;output;
      end;
      end;
drop i;
cards;
22.6 19.1 18.9 19.0  22.8 22.8 13.6 16.9  21.0 24.5 17.2 17.6  16.9 18.0 15.0 14.8  
20.0 15.2 16.6 13.1  21.9 15.4 14.2 16.9  21.5 20.1 16.7 16.2  21.2 21.2 19.6 14.8
;
proc sort; by season;
proc univariate normal;
     var x;
       by season;
run;
proc anova ;
      class season;
        model x=season;
run;

 

 image

第三题:使用单因素方差分析。
1、正态性校验:
a=1(春)时,P=0.0829
a=2(夏)时,P=0.8013
a=3(秋)时,P=0.7348
a=4(冬)时,P=0.9216
结论:P值皆大于0.05,不拒绝来自正态的假设,说明数据为正态分布。
2、方差分析:
结论:F=7.88  P=0.0006<0.01 说明有高度统计意义,即不同季节氯化物含量有差别。

 

image

DATA test2; 
  DO a=1 TO 2; /*a=1:治疗前;a=2:治疗后*/
DO i=1 TO 9;
     INPUT x @@;
     OUTPUT;
   END; 
  END; 
DROP i; 
CARDS;
68 65 55 75 50 70 76 65 72
128 82 80 112 125 110 85 80 105
;
/*正态检验*/
PROC SORT;
BY a;
PROC UNIVARIATE NORMAL;
VAR x;
BY a;
/*方差分析*/
PROC ANOVA;
CLASS a;
MODEL x = a;
RUN;
QUIT;

 

第四题:使用单因素方差分析。
1、正态性校验:
a=1(治疗前)时,P=0.3085
a=2(治疗后)时,P=0.1131
结论:P值皆大于0.05,不拒绝来自正态的假设,说明数据为正态分布。
2、方差分析:
结论:F=23.64  P=0.0002<0.01 说明有高度统计意义,即治疗前后血红蛋白含量有显著性差别。

 

image

data test3;
do a =1 to 8;
do b=1 to 4;
input x @@; output;
end; 
end;
cards;
5.27 5.27 4.94 4.61
5.27 5.22 4.88 4.66
5.88 5.83 5.38 5.00
5.44 5.38 5.38 5.00
5.66 5.44 5.38 4.88
6.22 6.22 5.61 5.22
5.83 5.72 5.38 4.88
5.27 5.11 5.00 4.44
;
proc sort;
by a;
proc univariate normal;
var x;
by a;
proc sort;
by b;
proc   univariate normal;
var x;
by b;
proc anova;
class i;
model x=i;
run;
proc anova;
class a b;
model x =a b;
means a b /snk;
run;
quit;

 

      Source                      DF        Anova SS     Mean Square    F Value    Pr > F
                      a                            7      2.49487188      0.35641027      26.22    <.0001
                      b                            3      2.89038438      0.96346146      70.87    <.0001
0.001《0.05
不同时间的血糖 浓度有区别。

 

 

薛毅书p358(纸介质) 例7.11,要求使用SAS方法完成。

image

data test4;
do a =1 to 3;
do b=1 to 4;
do i=1 to 5;
input x @@; output;
end;  
end;
end;
cards;
23 25 21 14 15 20 17 11 26 21
16 19 13 16 24 20 21 18 27 24
28 30 19 17 22 26 24 21 25 26
19 18 19 20 25 26 26 28 29 23
18 15 23 18 10 21 25 12 12 22
19 23 22 14 13 22 13 12 22 19
;
drop i;
proc sort;
by a;
proc univariate normal;
var x;
by a;
proc sort;
by b;
proc   univariate normal;
var x;
by b;
proc anova;
class i;
model x=i;
run;
proc anova;
class a b;
model x =a b;
means a b /snk;
run;
quit;

image

结论:F=9.37  P=0.0003<0.01 说明有高度统计意义

即 a 是表不同的树种  是有差别
b 是表示不同区域  不同区域 是没有差别的

您可以选择一种方式赞助本站