4

を使用してネストされたランダム効果の2つのモデルを比較しanova()たとすると、結果は次のようになります。

new.model: new
current.model: new 
              Df    AIC    BIC  logLik  Chisq Chi Df Pr(>Chisq)    
new.model      8 299196 299259 -149590                             
current.model  9 299083 299154 -149533 115.19      1  < 2.2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

テーブル部分のみを使用したい（下記参照）： 

              Df    AIC    BIC  logLik  Chisq Chi Df Pr(>Chisq)    
new.model      8 299196 299259 -149590                             
current.model  9 299083 299154 -149533 115.19      1  < 2.2e-16 ***

attribute（anova.object）$ heading = NULLを使用して見出しをnullに設定することで、見出し部分（blowを参照）を削除できることはわかっていますが、下部部分を削除する方法がわかりません：Signif 。コード：....。

new.model: new
current.model: new

空白のセルをNAに変更するため、data.frame（以下を参照）は使用したくありません。

data.frame(anova(new.model, current.model))
              Df      AIC      BIC    logLik    Chisq Chi.Df   Pr..Chisq.
new.model      8 299196.4 299258.9 -149590.2       NA     NA           NA
current.model  9 299083.2 299153.6 -149532.6 115.1851      1 7.168247e-27

皆さんはこの状況に対処する方法を知っているのだろうか。

[更新]：print.anovaを使用してラッパーを作成することになりました： 

anova.print = function(object, signif.stars = TRUE, heading = TRUE){
    if(!heading)
         attributes(object)$heading = NULL
    print.anova(object, signif.stars = signif.stars)
}

例：

dv = c(rnorm(20), rnorm(20, mean=2), rnorm(20))
iv = factor(rep(letters[1:3], each=20))
anova.object = anova(lm(dv~iv))

Analysis of Variance Table
Response: dv 
          Df Sum Sq Mean Sq F value    Pr(>F)    
iv         2 46.360 23.1798  29.534 1.578e-09 ***
Residuals 57 44.737  0.7849                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 


anova.print(anova.object, F, F)
          Df Sum Sq Mean Sq F value    Pr(>F)
iv         2 46.360 23.1798  29.534 1.578e-09
Residuals 57 44.737  0.7849
4

2 に答える 2

10

編集：anovaパラメータとしてsignif.starsを使用したprintメソッドがあります

anova(new.model, current.model, signif.stars=FALSE)

> x <- anova(lm(hp~mpg+am, data=mtcars))
> print(x, signif.stars=F)
Analysis of Variance Table

Response: hp
          Df Sum Sq Mean Sq F value    Pr(>F)
mpg        1  87791   87791 54.5403 3.888e-08
am         1  11255   11255  6.9924   0.01307
Residuals 29  46680    1610
于 2012-05-17T15:50:58.853 に答える
3

先日、NAを表示しないことについて同様の投稿がありました。あなたができること：

x <- as.matrix(anova(new.model, current.model))
print(x, na.print="", quote=FALSE)

mtcarsデータセットを使用したより再現性の高い例：

x <- as.matrix(anova(lm(hp~mpg+am, data=mtcars)))
print(x, na.print="", quote=FALSE)
于 2012-05-17T16:04:13.637 に答える