[R] variance using lmer

Per Theilgaard peth at dca.upv.es
Fri Apr 28 12:27:16 CEST 2006


Dear R help
I have a question on the variance of the binomial probit model.

I have fitted the following model :

> lmer1<-lmer(mp ~ l + op + l*op+ us_lev + bw_lev +(1|tatu) ,  
+                family = binomial(link="probit"), 
+                method = 'Laplace',
+                data = matings,
+                msVerbose= True)
> summary(lmer1)
Generalized linear mixed model fit using Laplace 
Formula: mp ~ l + op + l * op + us_lev + bw_lev + (1 | tatu) 
   Data: matings 
 Family: binomial(probit link)
      AIC      BIC    logLik deviance
 1455.121 1556.387 -706.5605 1413.121
Random effects:
     Groups        Name    Variance    Std.Dev. 
       tatu (Intercept)     0.42905     0.65502 
# of obs: 918, groups: tatu, 130

Estimated scale (compare to 1)  0.8521378 

I suppose that the dispersion parameter has been restricted to be 1.
Then the variance would be 1 + the variance of the random effect: 1+ 0.42905 ???

Thany you very much in advance
Per




More information about the R-help mailing list