Inclusion of correlation structure glmer function of package lme4
Dear Marcos, You cannot add correlation structures to an lme4 model like you can with nlme models. Have a look at the INLA or glmmTMB packages instead. Best regards, Thierry ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance thierry.onkelinx at inbo.be Havenlaan 88 bus 73, 1000 Brussel www.inbo.be /////////////////////////////////////////////////////////////////////////////////////////// To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey /////////////////////////////////////////////////////////////////////////////////////////// <https://www.inbo.be> Op wo 19 feb. 2020 om 15:15 schreef Marcos Paulo <cjfr2020 at gmail.com>:
Dear,
Good night.
I am using the glmer function of the lme4 package of R. My answer is
Bernoulli and my data are longitudinal.
In the bild Package, I have already found article using correlation
structures.
Is there any way I can use the glmer function and be able to include a
correlation structure? If it is not possible, could you tell me which
default is?
The code in R is only to facilitate the help.
## generalized linear mixed model
library(lme4)
(gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
data = cbpp, family = binomial))
Thanks
Marcos Paulo
[[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models