[R-sig-ME] Model is nearly unidentifiable with glmer
David Duffy
David.Duffy at qimr.edu.au
Wed Feb 26 01:26:23 CET 2014
On Tue, 25 Feb 2014, Ramon Diaz-Uriarte wrote:
> http://ligarto.org/rdiaz/nearly-unidentifiable.tar.gz
>
> a small (230 KB) file with:
>
> - data
> - R code to show the problem
> - two runs of that code on different machines
I don't seem to get those warnings when I run your "f02" ("very large
eigenvalue") problem here - get the same parameter estimates (x86_64,
lme4_1.0-6, Matrix_1.1-0). The estimated RE variance is very small.
Your cut-down "singularity" example gives a variance for dataSetID with a
starting value for theta, eg glmer(nS ~ Method + (1|dataSetID),
family=poisson, data = d2, start=0.01) - again v. small.
| David Duffy (MBBS PhD)
| email: David.Duffy at qimrberghofer.edu.au ph: INT+61+7+3362-0217 fax: -0101
| Genetic Epidemiology, QIMR Berghofer Institute of Medical Research
| 300 Herston Rd, Brisbane, Queensland 4006, Australia GPG 4D0B994A
More information about the R-sig-mixed-models
mailing list