2

私は多くの ANOVA テーブルを実行していますが、最終的には xtable パッケージの 'xtable' 関数を使用してそれらを実行し、LaTeX にエクスポートしたいと考えています。ただし、テーブル オブジェクトを使用して xtable を提供する必要があります。ANOVA オブジェクトは受け入れられません。私は基本的にANOVAテーブルをテーブルオブジェクトとして作りたいと思っています。再現可能なコードを次に示します。

utils::data(npk, package="MASS")
npk.aovE <- aov(yield ~  N*P*K + Error(block), npk)
summary(npk.aovE) ## THIS IS THE TABLE I WANT AS A TABLE OBJECT

通常の容疑者(as.table、print、およびxtable(summary(npk.aoE)))をすべて試しましたが、成功しませんでした。どんな助けでも大歓迎です

4

1 に答える 1

1

一般に、人々が望むのは、次のように取得された行列です。

 coef( summary(npk.aovE) )  # which returns NULL

ヘルプ ページにあるように、「関数coefは、標準誤差、t 統計量、および p 値を含む係数の行列を抽出します。」残念ながら、理論と実践は必ずしも一致しません。その概要オブジェクトは、実際には 2 つのデータフレームであり、その動作は で説明されてい?summary.aovlistます。

> summary(npk.aovE)[[2]]
          Df Sum Sq Mean Sq F value  Pr(>F)   
N          1 189.28  189.28  12.259 0.00437 **
P          1   8.40    8.40   0.544 0.47490   
K          1  95.20   95.20   6.166 0.02880 * 
N:P        1  21.28   21.28   1.378 0.26317   
N:K        1  33.14   33.14   2.146 0.16865   
P:K        1   0.48    0.48   0.031 0.86275   
Residuals 12 185.29   15.44                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
> summary(npk.aovE)[[1]]
          Df Sum Sq Mean Sq F value Pr(>F)
N:P:K      1   37.0   37.00   0.483  0.525
Residuals  4  306.3   76.57               
于 2012-09-21T06:07:12.553 に答える