Displaying 1 result from an estimated 1 matches for "myglm4".
Did you mean:
myglm
2006 Mar 08
1
Want to fit random intercept in logistic regression (testing lmer and glmmML)
...text ( 0.5*max(x), 0.5, expression( Prob( y[i] == 1) == frac( 1 , 1 +
exp(-eta[i] ))))
### Parameter estimates go to hell, as expected
myglm3 <- glm ( y ~ x, family=binomial(link="logit") )
summary(myglm3)
### Why doesn't quasibinomial show more evidence of the random intercept?
myglm4 <- glm(y~x, family=quasibinomial)
summary(myglm4)
# Can I estimate with lmer?
library(lme4)
### With no repeated observations, what does lmer want?
### Clusters of size 1 ?
### In lme, I'd need random= ~ 1
idnumber <- 1: length(y)
mylmer1 <- lmer( y ~ x + ( 1 | idnumber ), fami...