* Stata
xtset id time
mixed y x##c.time || id: time, reml      // 随机截距 + 随机斜率
mixed y x##c.time || id:                 // 仅随机截距（供比较）
estat ic                                 // AIC/BIC 比较
* 二分类结局改用：melogit y x || id:

## R
library(lme4); library(lmerTest)         // lmerTest 提供 p 值
fit <- lmer(y ~ x*time + (1 + time | id), data = df)
summary(fit)
ranova(fit)                              // 随机效应结构的似然比检验