4

Multivariate Tests: Site以下のMWEでfm1の出力からその部分を抽出する方法が気になります。

library(car)
fm1 <- summary(Anova(lm(cbind(Al, Fe, Mg, Ca, Na) ~ Site, data=Pottery)))
fm1

Type II MANOVA Tests:

Sum of squares and products for error:
           Al          Fe          Mg          Ca         Na
Al 48.2881429  7.08007143  0.60801429  0.10647143 0.58895714
Fe  7.0800714 10.95084571  0.52705714 -0.15519429 0.06675857
Mg  0.6080143  0.52705714 15.42961143  0.43537714 0.02761571
Ca  0.1064714 -0.15519429  0.43537714  0.05148571 0.01007857
Na  0.5889571  0.06675857  0.02761571  0.01007857 0.19929286

------------------------------------------

Term: Site 

Sum of squares and products for the hypothesis:
            Al          Fe          Mg         Ca         Na
Al  175.610319 -149.295533 -130.809707 -5.8891637 -5.3722648
Fe -149.295533  134.221616  117.745035  4.8217866  5.3259491
Mg -130.809707  117.745035  103.350527  4.2091613  4.7105458
Ca   -5.889164    4.821787    4.209161  0.2047027  0.1547830
Na   -5.372265    5.325949    4.710546  0.1547830  0.2582456

Multivariate Tests: Site
                 Df test stat  approx F num Df   den Df     Pr(>F)    
Pillai            3   1.55394   4.29839     15 60.00000 2.4129e-05 ***
Wilks             3   0.01230  13.08854     15 50.09147 1.8404e-12 ***
Hotelling-Lawley  3  35.43875  39.37639     15 50.00000 < 2.22e-16 ***
Roy               3  34.16111 136.64446      5 20.00000 9.4435e-15 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
4

2 に答える 2

3

fm1$multivariate.tests出力のSite一部に移動しfm1ます。

cat次に、との組み合わせを使用して、capture.output適切な印刷を行うかcapture.output、文字ベクトルのみを使用できます。

> cat(capture.output(fm1$multivariate.tests)[18:26], sep = "\n")
#
# Multivariate Tests: Site
#                  Df test stat  approx F num Df   den Df     Pr(>F)    
# Pillai            3   1.55394   4.29839     15 60.00000 2.4129e-05 ***
# Wilks             3   0.01230  13.08854     15 50.09147 1.8404e-12 ***
# Hotelling-Lawley  3  35.43875  39.37639     15 50.00000 < 2.22e-16 ***
# Roy               3  34.16111 136.64446      5 20.00000 9.4435e-15 ***
# ---
# Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

更新：の結果から

unlist(fm1$multivariate.tests, recursive = FALSE)

結果が数値として簡単にアクセスできるようには見えません。それで、あなたが要求したように、結果を行列に操作するために必要なものは次のとおりです。これを行ってからuser20650の回答を見たので、その提案に従い、ANOVAテーブルを介して値を取得することをお勧めします.

co <- capture.output(fm1$multivariate.tests)[20:24]
s <- strsplit(gsub("([*]+$)|[<]", "", co[-1]), "\\s+")
dc <- do.call(rbind, lapply(s, function(x) as.numeric(x[-1])))
row.names(dc) <- sapply(s, "[", 1)
s2 <- strsplit(co[1], " ")[[1]]
s2 <- s2[nzchar(s2)]
s3 <- s2[-c(1, length(s2))]
colnames(dc) <- c(s2[1], paste(s3[c(TRUE, FALSE)], s3[c(FALSE, TRUE)]), s2[10])
dc
#                  Df test stat  approx F num Df   den Df     Pr(>F)
# Pillai            3   1.55394   4.29839     15 60.00000 2.4129e-05
# Wilks             3   0.01230  13.08854     15 50.09147 1.8404e-12
# Hotelling-Lawley  3  35.43875  39.37639     15 50.00000 2.2200e-16
# Roy               3  34.16111 136.64446      5 20.00000 9.4435e-15

私の 2 番目のコード チャンクを改善したいと思う人がいたら、遠慮なく言ってください。

于 2014-09-17T19:39:32.907 に答える