【数据求助】请问，有没有人用R做过多层线性模型 （Hierarchical Linear Model）

本人刚开始学R，我想用R软件做一个4层 HLM，请大家给提提建议。
这个模型是不同省市的6年的学生成绩，第一层是学生，第二层是学校，第三层是一个growth model， 第四层是省市。CITY_ID是省市的编号，SCH_ID是学校的编号， YEAR是测量的年份。


我主要关注的是第3层的intercept（Π000k)和时间的slope（Π001k)，所以只在第一二层加两个变量，公式如下：



以下是我写的R的代码，请问这样写对不对？

#intall package////
library(Matrix)
library(lme4)
library(lmerTest)
library(dplyr)
library(reghelper)

#open file//////
data <- read.csv(file="D:stuschr.csv")
summary(data)

#centering quantitative variable/////
data$AGE_cen<- data$AGE-mean(data$AGE, na.rm=TRUE)
data$SCHSES_cen<- data$SCHSES-mean(data$SCHSES, na.rm=TRUE)

#model ////
conditional <- lmer(READING~ TIME + SCHSES_cen + SCHTYPE  +GENDER +AGE_cen  +(1+TIME|CITY_ID)+(1|CITY_ID:YEAR)+(1|CITY_ID:YEAR:SCH_ID), data=data,REML = FALSE, control=lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=1e5)))
summary(conditional)

数据结构和公式如下：

# 安装必要的包（如果尚未安装）
install.packages("lme4")    # 多层模型核心包
install.packages("lmerTest") # 为lme4提供p值
install.packages("ggplot2")  # 数据可视化
install.packages("effects")  # 模型效果可视化
install.packages("performance") # 模型诊断和性能评估

# 加载这些包
library(lme4)
library(lmerTest)
library(ggplot2)
library(effects)
library(performance)

hub957880 发表于 2025-9-10 17:45
# 安装必要的包（如果尚未安装）
install.packages("lme4")    # 多层模型核心包
install.packages("lme ...

已经安装了

还是路过。

为啥不用python？

Killua609 发表于 2025-9-14 10:36
为啥不用python？

零基础，先学了R，发现可以出结果，但是不知道模型写得对不对