我正在关注混合效果模型的教程。本教程使用来自 mlmRev 包的 egsingle 数据集。作为教程的一部分,作者使用 groupedData() 作为:
egsingle <- groupedData(math ~ year | schoolid/childid, data = egsingle)
有人可以帮我理解“schoolid/childid”是什么意思吗?
请注意,schoolid 和 childid 都是因素!
此外,在本教程的后面,作者采用大小为 50 的样本并使用 lmList() 来拟合每个主题的 OLS 回归,方法是:
egsingle <- groupedData(math ~ year | schoolid/childid, data = egsingle)
samp <- sample(levels(egsingle$childid), 50)
level2.subgroup <- subset(egsingle, childid %in% samp)
# fitting a separate OLS regression line to each student
level2 <- lmList(math ~ year | childid, data = level2.subgroup)
plot(augPred(level2))
当我运行上面的 lmList 命令时,我收到以下错误:
Error in eval(expr, envir, enclos) : object 'childid' not found
In addition: Warning messages:
1: In lmList(math ~ year | childid, data = level2.subgroup) :
lmList does not (yet) work correctly on groupedData objects
2: In Ops.factor(schoolid, childid) : ‘/’ not meaningful for factors
谁能帮我弄清楚为什么会出现这些错误?