我正在使用R的lme4包的glmer函数。我的答案是伯努利,我的数据是纵向的。
在bild包中,我已经找到了使用Mc1或MC2作为相关结构的文章。
有什么办法可以使用glmer函数并且能够包含相关结构吗?如果不可能,你能告诉我哪个是默认值吗?
R中的代码只是为了方便帮助。
library(dplyr)
n = 300
xx<-c("r1","r2","r3","r4","r5")
xxx<-c("e1","e2","e3")
p=0.3
df1 <- data_frame(
xx1 = runif(n, min = 0, max = 10),
xx2 = runif(n, min = 0, max = 10),
xx3 = runif(n, min = 0, max = 10),
School = factor(sample(xxx, n,re=TRUE)),
yx = as.factor(rbinom(n, size = 1, prob = p))
)
df1
library(lme4)
mm2 <- glmer(yx ~ xx1 + xx2 + xx3 + (1 | School), data = df1,
family = "binomial",control = glmerControl(calc.derivs = FALSE))