孙振球主编的《医学统计学》(人民卫生出版社, 2002:211)上有这样一道题,我用sas计算的的总方差要大于总均数的离均差平方和,请教高手这是为什么?按我的理解总方差(即总均数的离均差平方和)应该等于下图中红色部分之和(sas没有直接给出)。结果为60121.2035。但是你将原始数据在excel中计算总均数的离均差平方和(可以用DEVSQ函数求),结果为60006.6885。随后我按书上的公式在excel中分别计算各项,发现差别在time项。sas的SStime为27893.95516,按公式excel中的SStime为27779.440142。excel中的SS总为60006.6885。sas程序如下:
data aa;
input group $ T1-T4;
cards;
A 90.53 142.12 65.54 73.28
A 88.43 163.17 48.95 71.77
A 100.01 144.75 86.06 80.01
A 46.32 126.33 48.95 39.54
A 73.69 138.96 70.02 60.89
A 105.27 126.33 75.01 83.66
A 86.32 121.06 78.95 70.24
B 70.53 97.38 112.12 28.5
B 68.43 95.27 133.17 56.9
B 57.37 78.43 83.16 48.34
B 105.8 120.54 136.33 84.03
B 80.01 104.75 114.75 65.61
B 56.32 75.27 96.33 47.52
B 53.69 110.02 138.96 45.44
B 85.27 110.01 126.33 69.47
B 66.32 115.27 129.06 55.29
;
proc glm;
class group;
model T1-T4=group/nouni;
repeated time 4 contrast(1);
run;
quit;