sas-方差分析.doc

-
. z.
四、统计推断Ⅱ〔方差分析——多个平均数的比拟〕〔1〕发布：admin 时间：2006-8-26
四、统计推断Ⅱ〔方差分析——多个平均数的比拟〕(1)
方差分析是关于程序的第一行称为显示管理命令，相当于
-
. z.
DOS中的批处理命令，它的作用是去除LOG窗和OUTPUT窗中的容；DROP r 用于删除变量r，因为r在后面的分析中用不到。当然，DM和DROP两句都可以不用。
3．系组分组资料
[例11] 。
DATA tb619;
DO trt=1 TO 4;
DO m=1 TO 3;
DO nop=1 TO 4;
INPUT y;
OUTPUT;
END;
END;
END;
CARDS;
50 55 40 35 35 35 30 40 45 40 40 50
50 45 50 45 55 60 50 50 55 45 65 55
85 60 90 85 65 70 80 65 70 70 70 70
60 55 35 70 60 85 45 75 65 65 85 75
;
PROC FORMAT;
VALUE trtv 1=〞A〞 2=〞B〞 3=〞C〞 4=〞D〞;
PROC ANOVA;
FORMAT trt trtv.;
CLASS trt m;
MODEL y=trt m(trt);
TEST H=trt E=m(trt);
MEANS trt / DUNCAN E=m(trt);
RUN;
这是一个系统分组资料，用ANOVA过程，其模型中效应包括培养液的效应〔trt〕、同一培养液盆间差异[m(trt)]，以及同一盆各植株间的差异〔即试验误差，在MODEL中不列出〕。MODEL这种定义方式称嵌套模型。
TEST语句指出测验trt的显著性时应使用m(trt)作误差项，假设不指明则以试验误差作为被比量进展测验。因此方差分析表中的trt项F测验不正确，必须使用TEST语句定义被比量进展测验。
对各培养液平均数的多重比拟用新复极差法〔DUNCAN〕，且平均数的标准误由m(trt)效应项方差计算。
本程序中还应用了FORMAT过程进展变量格式设定，VALUE则给出具体容，这种设定通过语句“FORMAT trt trtv.;〞指定输出格式。当然格式设定语句可以不要，这里只是为了使处理名称与相应效应直接对应。
五、统计推断Ⅱ〔方差分析——多个平均数的比拟〕〔2〕发布：admin 时间：2006-8-26
-
. z.
(二) 两向分组资料
1．组合只有单个观察值的两向分组资料
〔单因素随机区组或二因素完全随机设计〕
[例12] 。
DATA tb625;
DO szs=1 TO 5;
DO blk=1 TO 4;
INPUT y ;
OUTPUT;
END;
END;
CARDS;
60 62 61 60
65 65 68 65
63 61 61 60
64 67 63 61
62 65 62 64
61 62 62 65
;
PROC ANOVA;
CLASS szs blk;
MODEL y=szs blk;
MEANS szs / T;
MEANS szs/DUNNETT;
RUN;
本例用ANOVA过程，效应项包括生长素效应和区组效应以及省去的误差。对生长素效应分别用LSD法和DUNNETT氏最小显著差数法〔DLSD法，默认对照是处理项中第一项〕进展多重比拟。
2．组合有重复的两向分组资料的方差分析〔二因素完全随机设计〕
[例13] 。
DATA tb631;
DO a=1 TO 3;
DO b=1 TO 3;
DO r=1 to 3;
INPUT y ;
OUTPUT;
EN