[R] variance using lmer

From: Per Theilgaard <peth_at_dca.upv.es>
Date: Fri 28 Apr 2006 - 20:27:16 EST


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



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Apr 28 20:37:26 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 05 May 2006 - 14:09:58 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.