[R-sig-ME] how to extract the BIC value

[Andy Fugard (Work)]

Greetings,

Assuming you're using lmer, here's an example which does what you need:

> (fm1    <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
Linear mixed model fit by REML
Formula: Reaction ~ Days + (Days | Subject)
   Data: sleepstudy
  AIC  BIC logLik deviance REMLdev
 1756 1775 -871.8     1752    1744
Random effects:
 Groups   Name        Variance Std.Dev. Corr
 Subject  (Intercept) 612.092  24.7405
          Days         35.072   5.9221  0.066
 Residual             654.941  25.5918
Number of obs: 180, groups: Subject, 18

Fixed effects:
            Estimate Std. Error t value
(Intercept)  251.405      6.825   36.84
Days          10.467      1.546    6.77

Correlation of Fixed Effects:
     (Intr)
Days -0.138

> (fm1fit <- summary(fm1)@AICtab)
      AIC      BIC    logLik deviance  REMLdev
 1755.628 1774.786 -871.8141 1751.986 1743.628

> fm1fit$BIC
[1] 1774.786

Cheers,

Andy


Tahira Jamil wrote:
> Hi
> I can extract the AIC value of a model like this
> 
> AIC(logLik(fm0)
> 
> How can I extract the BIC value if I need!
> 
> Cheers
> Tahira 
> Biometris 
> Wageningen University
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models


-- 
Andy Fugard, Postdoctoral researcher, ESF LogICCC project
"Modeling human inference within the framework of probability logic"
Department of Psychology, University of Salzburg, Austria
http://www.andyfugard.info