Displaying 1 result from an estimated 1 matches for "myglmm1".
Did you mean:
mygamm1
2006 Mar 08
1
Want to fit random intercept in logistic regression (testing lmer and glmmML)
...f size 1 ?
### In lme, I'd need random= ~ 1
idnumber <- 1: length(y)
mylmer1 <- lmer( y ~ x + ( 1 | idnumber ), family=binomial, method="Laplace" )
summary(mylmer1)
### Glmm wants clusters, and I don't have any clusters with more than
1 observation
###
library(glmmML)
myglmm1 <- glmmML(y~x, family=binomial, cluster = idnumber )
summary(myglmm1)
--
Paul E. Johnson
Professor, Political Science
1541 Lilac Lane, Room 504
University of Kansas