我正在尝试报告我们在 R 中拟合的线性混合模型的效果大小。现在我正在研究报告部分 eta 平方或 eta 平方。但是,为此我需要计算平方误差之和。我们正在使用该lme()函数,它不报告 MSE 或效果大小。我不是进行主要分析的人,所以我不知道我们是否可以切换到使用ez包,如Omega squared for measure of effect in R 中所述?.
lmeObject函数返回的有lme()一些信息,但我不确定哪个最合适。在 上使用anova()时lmeObject,它将 F 检验的分母自由度报告为 56683,因此使用 MSE 的值,我可以计算 SSE 并对 F 检验进行逆向工程,以获得部分 eta 平方所需的 SS 处理。我有 2 个固定效应和一个随机效应(用于重复测量设计)。
我已经查看lmeObject$sigma并计算了 的平方和lmeObject$residuals[,1],但他们不同意(我sigma将 SS 除以自由度)。
那里有任何 R 大师可以告诉如何从 计算 MSE 或 SSE lmeObject?