X1 X2 X3是12名中学生的身高、体重、胸围,要检验现在中学生身体发育情况与10年前对比有无存在区别,使用多元方差分析的方法,将三者看成一个衡量身体素质的一个整体指标。
程序如下:
data ch2_2;
input x1 x2 x3 @@;
y1=x1-157.54;y2=x2-45.12;y3=x3-70.02;
cards;
171.0 58.5 81.0 154.8 44.5 74.0
175.0 65.0 87.0 164.0 51.0 72.0
159.0 38.0 71.0 165.2 55.0 79.0
155.3 45.0 74.0 164.5 46.0 71.0
152.0 35.0 63.0 159.1 48.0 72.5
158.3 44.5 75.0 164.2 46.5 73.0
;
run;
proc glm;
model y1 y2 y3= /nouni;
manova h=intercept/printe printh;
run;
问题本体:
根据Y1 Y2 Y3可以计算出一个3*3的协方差矩阵(或者SSCP),但是现在我只有一个矩阵,怎么去分解它得到ERROR的矩阵然后用一般方差分析的思想分析出是否存在区别呢。
我特别不理解H=intercepet这一句,我感觉这一句程序是解决我问题的关键。但是我看不懂这个intercept的矩阵从何而来。
请各位大大们指教。没有分组变量的情况下,SAS是如何完成这一次多元方差分析的呢?