## Data originally from Davis 1991 Stat. Med., as packaged in geepack
## and transformed (center, id -> factor, idctr created, levels labeled)
library(lme4)
if (.Platform$OS.type != "windows") {
load(system.file("testdata","respiratory.RData",package="lme4"))
m_glmer_4.L <- glmer(outcome~center+treat+sex+age+baseline+(1|idctr),
family=binomial,data=respiratory)
m_glmer_4.GHQ5 <- glmer(outcome~center+treat+sex+age+baseline+(1|idctr),
family=binomial,data=respiratory,nAGQ=5)
m_glmer_4.GHQ8 <- glmer(outcome~center+treat+sex+age+baseline+(1|idctr),
family=binomial,data=respiratory,nAGQ=8)
m_glmer_4.GHQ16 <- glmer(outcome~center+treat+sex+age+baseline+(1|idctr),
family=binomial,data=respiratory,nAGQ=16)
} ## skip on windows (for speed)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.