如何解决软件包“ 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 ...
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。