#----两水平logistic回归模型

#----两水平logistic回归模型-------------
install.packages("lmerTest")
install.packages("sjstats")
install.packages("Matrix")
install.packages("lme4", type = "source")
library(Matrix)
library(lmerTest)
library(sjstats)
library(lme4)

data <- read.csv("数据.csv",header = TRUE)
data <- read.csv("C:/Users/ASUS/Desktop/两水平logistics回归模型/两水平logistics回归模型/数据.csv")

head(data)
data$cancertage<-factor(data$cancertage,levels=c(1:4),labels=c("I","II","III","IV"))

#零模型，没有加入解释变量
m_ri <- glmer(remission ~ 1+
                (1 | DID), data = data, family = binomial, control = glmerControl(optimizer = "bobyqa"))
summary(m_ri)

#引入解释变量拟合两水平logistic回归模型
m_ri2 <- glmer(remission ~ cancertage +experience1+
                (1 | DID), da