how to extract the BIC value

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

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