系统分组或嵌套设计及其资料的统计分析
PROC SORT;
BY a b c r;
RUN;
PROC NESTED;
CLASS a b c;
VAR y;
RUN;
PROC GLM;
CLASS a b c;
MODEL y=a b(a) c(a b) / SS1;
TEST H=a E=b (a) ;
TEST H=b(a) E=c(a b);
RUN;
1)r代表试验批次;a,b,c代表3个因素,是分组变量;
2)用SAS提供的专门对系统分组(或嵌套)设计资料进行ANOVA的过程─NESTED时,虽过程步写起来很简单,但它受数据步中读取分组变量先后顺序的 ...
系统分组或嵌套设计及其资料的统计分析
PROC SORT;
BY a b c r;
RUN;
PROC NESTED;
CLASS a b c;
VAR y;
RUN;
PROC GLM;
CLASS a b c;
MODEL y=a b(a) c(a b) / SS1;
TEST H=a E=b (a) ;
TEST H=b(a) E=c(a b);
RUN;
1)r代表试验批次;a,b,c代表3个因素,是分组变量;
2)用SAS提供的专门对系统分组(或嵌套)设计资料进行ANOVA的过程─NESTED时,虽过程步写起来很简单,但它受数据步中读取分组变量先后顺序的影响很大,若不先用SORT过程对分组变量重新整理(因代表批次的变量r在最前面),则得出错误结果。
3)SORT过程对因素按主次重新排队; NESTED过程适用于对因素已按主次排列且来自系统分组设计的平衡资料的方差分析;第3个过程─GLM过程适用面很宽,资料可以是不平衡的, 也可不排序,仅取决于MODEL语句和TEST语句的写法,但写起来很复杂,需要的机器内存较大。对于系统分组设计资料,MODEL语句后只能选SS1(计算结果与因素在该语句中的顺序有关),不能选SS3。平衡资料只需用前2个过程,也可仅用最后1个过程;若不是平衡资料, 必须用最后1个过程。