0

lme4 ライブラリを使用して anova モデルを実行しようとしています。実際のデータがないので、MASS から Cars93 データをヤンクしました。lmer(...) を機能させることはできましたが、anova テーブルを印刷できないようです。

何が起こったかは次のとおりです。

> library(MASS)

> cars <- Cars93

> library(lme4)
Loading required package: Matrix
Loading required package: lattice

Attaching package: ‘lme4’

The following object(s) are masked from ‘package:stats’:

AIC, BIC

> fit <- lmer(Max.Price ~ 1+(1|AirBags)+(1|Man.trans.avail), data=cars)
> anova(fit)
Analysis of Variance Table
Df Sum Sq Mean Sq F value   
> summary(fit)
Linear mixed model fit by REML 
Formula: Max.Price ~ 1 + (1 | AirBags) + (1 | Man.trans.avail) 
Data: cars 
AIC   BIC logLik deviance REMLdev
688 698.1   -340    685.1     680
Random effects:
Groups          Name        Variance Std.Dev.
AirBags         (Intercept) 62.2339  7.8888  
Man.trans.avail (Intercept)  6.3022  2.5104  
Residual                    83.2786  9.1257  
Number of obs: 93, groups: AirBags, 3; Man.trans.avail, 2

Fixed effects:
            Estimate Std. Error t value
(Intercept)    23.81       5.00   4.761

どんな助けや提案も大歓迎です。

4

1 に答える 1

1

anovaが機能するには、式に変数が必要だと思います。例えば

fit <- lmer(Max.Price ~ Weight + (1 | AirBags) + (1 | Man.trans.avail),data=cars)

anova(fit)
Analysis of Variance Table
       Df Sum Sq Mean Sq F value
Weight  1 1932.1  1932.1   29.67
于 2013-03-01T07:42:30.520 に答える