<<目次へ

最終更新日　：　2007.10.31

モデルの検査�T：独立性

　p141のBOX 8.4にあるMANOVA（多変量分散分析）をRで行うにはmanova()という関数を用いる。

# ExcelよりWEEKとLGWTの列のデータを取り込む
> dat <- excel.w(2)
Read 80 items

# seq関数を利用してLOGWT3から60までをベクトルにする
> LOGWT3 <- dat[seq(1,40,4), 2]
> LOGWT8 <- dat[seq(2,40,4), 2]
> LOGWT20 <- dat[seq(3,40,4), 2]
> LOGWT60 <- dat[seq(4,40,4), 2]
> DIET <- as.factor(rep(c(1,2), c(5,5)))

# manovaを用いて多変量分散分析を行う
> result <- manova(cbind(LOGWT3, LOGWT8, LOGWT20, LOGWT60) ~ DIET)
> summary(result, test="Wilks")
          Df   Wilks approx F num Df den Df Pr(>F)
DIET       1 0.58209  0.89745      4      5 0.5286
Residuals  8                                      
> summary(result, test="Hotelling-Lawley")
          Df Hotelling-Lawley approx F num Df den Df Pr(>F)
DIET       1          0.71796  0.89745      4      5 0.5286
Residuals  8                                               
> summary(result, test="Pillai")
          Df  Pillai approx F num Df den Df Pr(>F)
DIET       1 0.41791  0.89745      4      5 0.5286
Residuals  8                                      
> summary(result, test="Roy")
          Df     Roy approx F num Df den Df Pr(>F)
DIET       1 0.71796  0.89745      4      5 0.5286
Residuals  8