软件包“ mitml”中的多级R平方错误

问题描述

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

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

nlme： 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平方统计的计算

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

解决方法

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

解决它

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

请注意，程序包的启动消息确实显示

***这是Beta版软件。请报告所有错误！

该软件包的“问题”（错误报告）列表为here ...