我这样理解你的问题,在体重、种类一致的前提下,不同剂量(A-E)对胰腺重量的影响。
代码如下:
data a1;
input class I $ II $ III $ IIII $ IIIII $ @@;
cards;
1 C E A B D
1 65 85 57 49 79
2 E B D C A
2 82 63 77 70 46
3 A D C E B
3 73 68 51 76 52
4 D C B A E
4 92 67 63 41 68
5 B A E D C
5 81 56 99 75 66
;
run;
proc transpose data=a1 out=a2(RENAME=(_NAME_=weight col1=dose));
by class;
var I II III IIII IIIII;
run;
data a3;
set a2;
pancreas=input( col2,10.);
drop col2;
run;
proc glm data=a3;
class class weight dose;
model pancreas= class weight dose;
lsmeans dose/pdiff;
run;
quit;
结果如下:
Output1:
pancreas LSMEAN
dose LSMEAN Number
A 54.6000000 1
B 61.6000000 2
C 63.8000000 3
D 78.2000000 4
E 82.0000000 5
说明:这个表是每个剂量组的胰腺重均值,你的问题就是比较这几个值是不是有显著差异吧?。
Output2:
Least Squares Means for effect dose
Pr > |t| for H0: LSMean(i)=LSMean(j)
Dependent Variable: pancreas
i/j 1 2 3 4 5
1 0.2505 0.1385
0.0016 0.0005
2 0.2505 0.7109 0.0143 0.0042
3 0.1385 0.7109 0.0287 0.0085
4 0.0016 0.0143 0.0287 0.5245
5 0.0005 0.0042 0.0085 0.5245
说明:这个表代表剂量A-E组两两进行均值比较的显著性,低于0.05可以认为差别是显著的。
例如红色0.0016,表示A组和D组均值比较,A组的54.6要显著低于D组78.2。
我也不是很懂,求高手来解释,哈。
我个人感觉,ABC可以归为一类,DE可以归为一类
