#----两水平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), data = data, family = binomial, control = glmerControl(optimizer = "bobyqa"))
summary(m_ri2)
performance::icc(m_ri2)
#与传统logistic模型比较
m2 <- glm(remission ~ cancertage+experience1, data = data, family = binomial)
summary(m2)
anova(m_ri2, m2)