[R-sig-ME] False convergence in null model
Iker Vaquero Alba
karraspito at yahoo.es
Tue Aug 23 19:44:22 CEST 2011
Dear list:
Sometimes I have had "false convergence" problems while performing some analyses with lmer. Usually, I was told that the problem could be an overparameterization. But this time, I've got the false convergence message just at the end of a stepwise model simplification, after removing the last remaining term. Obviously, I would like to know the reason for this, and I would appreciate any ideas about it. I attach the data set I am using. The analysis was performed after standardizing the continuous independent variables (tlength, brithr, briventral, inslarge and cond).
This is the step where I get the error message:
s.clutchmodel1.57<-update(s.clutchmodel1.56,~.-brithrst)
Mensajes de aviso perdidos
In mer_finalize(ans) : singular convergence (7)
> summary(s.clutchmodel1.57)
Generalized linear mixed model fit by the Laplace approximation
Formula: clsize1 ~ (1 | site/pair) + (1 | year)
AIC BIC logLik deviance
30 39.68 -11 22
Random effects:
Groups Name Variance Std.Dev.
pair:site (Intercept) 0 0
site (Intercept) 0 0
year (Intercept) 0 0
Number of obs: 83, groups: pair:site, 37; site, 12; year, 2
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.60702 0.04915 32.7 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> summary(s.clutchmodel1.56)
Generalized linear mixed model fit by the Laplace approximation
Formula: clsize1 ~ brithrst + (1 | site/pair) + (1 | year)
AIC BIC logLik deviance
31.52 43.62 -10.76 21.52
Random effects:
Groups Name Variance Std.Dev.
pair:site (Intercept) 0.0000e+00 0.0000e+00
site (Intercept) 7.6128e-17 8.7251e-09
year (Intercept) 1.0638e-14 1.0314e-07
Number of obs: 83, groups: pair:site, 37; site, 12; year, 2
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.60645 0.04918 32.67 <2e-16 ***
brithrst 0.03418 0.04930 0.69 0.488
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
(Intr)
brithrst -0.034
Do you think the problem could be in the fact that I am including random factors which can explain none of the variance? Should I do the analysis again without including those factors?
Thank you very much.
Iker
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Clutch Size1 2009+2010.txt
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20110823/fffe3b54/attachment-0003.txt>
More information about the R-sig-mixed-models
mailing list