1

我正在尝试使用 R 中的包 'mitml' 来计算多级模型的 R 平方度量。我尝试使用 lme4 和 nlme 来指定我的模型。但是，当我使用 lme4 指定模型时，所有 4 次计算中的 R 平方值都是相同的（不是我所期望的），当我使用 nlme 时出现错误。这是我用来跨两个包指定模型的代码：

lme4： H1 <- lmer(PAA_groupmc ~ Velocity.difference+(1|mydata$ID), data=mydata)

名词： h1.1 <- lme(PAA_groupmc ~ Velocity.difference, data=mydata, random = ~1|ID, method="REML")

PAA_groupmc 和 Velocity.difference 都是连续变量，ID 是代表每个人的因素，因为我有一个重复测量数据集。我允许按人进行随机拦截。

当我运行时，multilevelR2(H1)我得到以下结果： RB1：0.1004472
RB2：0.1004472
SB：0.1004472
MVP：0.1013596

问题1：我认为结果如此相似很奇怪，因为我没想到会是这样。有人可以解释为什么会发生这种情况或我可能做错了什么吗？

问题2：运行时multilevelR2(h1.1)出现以下错误：

multilevelR2(h1.1) 中的错误：类模型不支持计算 R 平方统计量

这个错误是什么意思，我该如何解决？

4

1 回答 1

1

对于您的第二个问题，该错误看起来像 中的错误mitml::multilevelR2，它将对象误解lme为模型列表（然后无法为列表的第一个元素找到合适的类）；您可以通过以下方式解决它

mitml:::.getRsquared(<your_model>,
                     print=c("RB1", "RB2", "SB", "MVP"), 
                     method="nlme")

请注意，包的启动消息确实说

*** 这是测试版软件。请报告任何错误！

包的“问题”（错误报告）列表在这里...

于 2020-09-13T15:47:56.983 回答