两因素随机区组实验的区组如何分组?

如题所述

一个两因素交叉分组实验,若每一处理重复n次,全部实验共abn次(见课本9.1.1)。这abn次实验的实验条件或实验材料必须具有同质性。否则,由于实验材料或实验条件的差异所引起的误差会混杂于实验误差中,影响试验结果的可靠性。为避免这种情况的发生,与随机化完全区组的做法一样,将每一套水平组合,安排在一个区组中,n次重复构成了n个区组。这样一种设计称为两因素随机化区组设计。统计模型为:

其中ai 、bj 和(ab)ij分别为A因素、B因素和AB交互作用效应,dk 是第k区组效应。设A因素为固定因素,B和区组为随机因素,模型中各分量的均方期望可由下表推演出。

因素
F R R R 均方期望
F R R R 均方期望

abn1
abn1

ijkl
ijkl

αi
0bn1

0bn1

βj
a1n1

a1n1

(αβ)ij
01n1

01n1

δk
ab11

ab11

(αδ)ik
0b11

1111

(βδ)jk
a111

(αβδ)ijk
0111

ε(ijk)l
1111

上表中的l为区组内的重复,因l = 1,这时s2无法估计(误差自由度dfe = 0)。假设区组与因素间不存在交互作用,即将上表左半部的后四行合并,作为误差估计,得到表的右半部。由均方期望可以得到检验统计量,FA = MSA / MSAB,FB = MSB / MSe,FAB = MSAB / MSe。两因素随机区组实验的方差分析与三因素交叉分组实验的方差分析程序基本相同。

例 2.13 课本上表9-11中的实验,共需32名同质受试者,因32名同质受试者很难找到,因此将实验的两个重复安排为两个区组,每一区组只要16名同质受试者即可。

解: 先创建一个名为a:\2-7data.dat的外部数据文件。SAS程序为:

options linesize = 76;

data work;

infile ‘a:\2-7data.dat’;

input block a b energy @@;

run;

proc anova;

class block a b;

model energy = block a b a*b;

test h = a e = a*b;

means a / duncan e = a*b;

run;
温馨提示:答案为网友推荐,仅供参考