1、manova 語(yǔ)句:當(dāng)效應(yīng)變量有多個(gè)時(shí)，此語(yǔ)句控制 anova 過程進(jìn)入多元方差分析模式，其后的選項(xiàng)用以指定多元方差分析時(shí)的各項(xiàng)指標(biāo)。第四題：data sas4;input groupvas1 vas2 vas3;cards;27 6626 6526 5325 4527 6626 5427 65645428 75;proc glm;class group;m vas1 vas2 vas3 =group ; manova h=group/prh pre summary;run;結(jié)果解釋：第一部分GLM 過程GLM 過程因變量: vas1源度平方和均方 F 值 Pr F模型1 3.20000000

2、 3.20000000 2.060.1686誤差18 28.00000000 1.55555556的觀測(cè)數(shù) 20使用的觀測(cè)數(shù) 20分類水平信息分類 水平 值group2 1 2GLM 過程因變量: vas2R 方 變異系數(shù)根 MSE vas2 均值0.001942 20.78081 1.1948975.750000源度平方和均方 F 值 Pr F模型1 0.05000000 0.05000000 0.040.8536誤差18 25.70000000 1.42777778校正合計(jì)19 25.75000000源度 III 型 SS均方 F 值 Pr Fgroup1 3.20000000 3.200

3、00000 2.060.1686源度I 型 SS均方 F 值 Pr Fgroup1 3.20000000 3.20000000 2.060.1686R 方 變異系數(shù)根 MSE vas1 均值0.102564 20.11644 1.2472196.200000源度平方和均方 F 值 Pr F校正合計(jì)19 31.20000000GLM 過程因變量: vas3源度I 型 SS均方 F 值 Pr Fgroup1 1.25000000 1.25000000 0.600.4474R 方 變異系數(shù)根 MSE vas3 均值0.032425 30.95745 1.4395224.650000源度平方和均方 F

4、 值 Pr F模型1 1.25000000 1.25000000 0.600.4474誤差18 37.30000000 2.07222222校正合計(jì)19 38.55000000源度 III 型 SS均方 F 值 Pr Fgroup1 0.05000000 0.05000000 0.040.8536源度I 型 SS均方 F 值 Pr Fgroup1 0.05000000 0.05000000 0.040.8536以上是各個(gè)變量組間比較的結(jié)果，相信大家都懂，就不做解釋了。第二部分：GLM 過程多元方差分析此部分是“E”矩陣，即三個(gè)指標(biāo)的誤差矩陣。給出了因變量對(duì)因子的偏相關(guān)系數(shù)及其顯著性。結(jié)果顯示v

5、as1 與vas2、vas1與 vas3、 vas2 與 vas3 的 P 值都小于 0.05，都顯示顯著正相關(guān)。GLM 過程誤差 SSCP 矩陣的偏相關(guān)系數(shù) / 概率 |r|度 = 18vas1vas2vas3 vas11.0000000.917043 0.53841.00010.0174vas20.9170431.000000 0.549070 Fgroup1 1.25000000 1.25000000 0.600.4474多元方差分析此部分顯示分組的三個(gè)指標(biāo)的平方和矩陣。顯示的是分組的差異。這部分三個(gè)指標(biāo)“H/E”矩陣，即組間平方和/向量積矩陣（H）差異除以組內(nèi)的誤差矩陣（E）得到的比例

6、。是檢驗(yàn)基于檢驗(yàn)組間的方差是否大于期望的組內(nèi)的方差。對(duì)角線就是組間平方和/單變量方差分析誤差平方和。假設(shè)無總體 group 效應(yīng)的 MANOVA 檢驗(yàn)準(zhǔn)則和精確 F 統(tǒng)計(jì)量H = III 型 SSCP 矩陣: group E = 誤差 SSCP 矩陣S=1 M=0.5 N=7以下對(duì)象的特性根和向量: E Inverse * H, where H = III 型 SSCP 矩陣: groupE = 誤差 SSCP 矩陣特性根 百分比 特性向量 VEV=1vas1vas2vas3 0.57156227 100.00 -0.46036699 0.44855242 -0.029325800.00000

7、0000.00 -0.11975321 -0.01595897 0.194796940.000000000.00 -0.02795959 0.22367671 0.00000000H = III 型 SSCP 矩陣: groupvas1vas2vas3vas13.20.42vas20.40.050.25vas320.251.25最主要看這部分結(jié)果：可以看出=0.63630950，轉(zhuǎn)化的近似F=3.05，其 P 值大于 0.05，無統(tǒng)計(jì)學(xué)意義。即 2 個(gè)處理組間無統(tǒng)計(jì)學(xué)意義，這是考慮了三個(gè)指標(biāo)的相關(guān)性在內(nèi)而得出來的結(jié)論。這是個(gè)整體結(jié)果。三個(gè)指標(biāo)整體存在組檢差異，具體是哪兩個(gè)有差異就看下面的結(jié)果第

8、三部分:因變量: vas1因變量: vas2源度 III 型 SS均方 F 值 Pr F源度III 型 SS均方 F 值 Pr F group1 3.20000000 3.20000000 2.060.1686誤差18 28.00000000 1.55555556統(tǒng)計(jì)量值F 值分子自分母Pr F由度度Wilks Lambda0.636309503.053160.0590Pillais Trace0.363690503.053160.0590Holing-Lawley Trace 0.571562273.053160.0590Roys Greatest Root0.571562273.05316

9、0.0590因變量: vas3結(jié)果顯示；三個(gè)指標(biāo)分別在組間均無統(tǒng)計(jì)學(xué)意義。GLM 過程最小二乘均值多重比較的調(diào)整: Bonferronigroup vas2 LSMEAN H0:LSMean1=LSMean2t 值 Pr |t| 15.70000000 -0.190.853625.80000000group vas1 LSMEAN H0:LSMean1=LSMean2t 值 Pr |t| 15.80000000 -1.430.168626.60000000源度III 型 SS均方 F 值Pr F group1 1.25000000 1.250000000.600.4474誤差18 37.30000000 2.07222222源度 III 型 SS均方 F 值 Pr Fgroup1 0.05000000 0.05000