[R] GLMM(..., family=binomial(link="cloglog"))?

[Peter Dalgaard]

Spencer Graves <spencer.graves at pdf.com> writes:

> Data: DF
>   log-likelihood:  -55.8861
> Random effects:
>   Groups Name        Variance   Std.Dev.
>   smpl   (Intercept) 1.7500e-12 1.3229e-06
> 
> Estimated scale (compare to 1)  3.280753
> 
> Fixed effects:
>              Estimate Std. Error z value Pr(>|z|)
> (Intercept) 0.148271   0.063419  2.3379  0.01939
> 
> Number of Observations: 10
> Number of Groups: 10

The only information you have about sigma is in the overdispersion
of y, so you probably cannot both have a scale parameter and a random
effect of sample in there. Doug did say something about fixing the
scale in GLMM, but I forget whether he had implemented it or planned
to...

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907