如何获得鲁棒混合效应模型（rlmer命令；robustlmm）的R^2？

我使用

robustlmm

软件包中的

rlmer

命令估计了一个稳健的混合效应模型。有没有办法获得边际和条件R^2值？

我自己来回答这个问题。我在R中找不到与

lmerMod

对象中的

R.squaredGLMM

等效的包或函数，但我找到了一种快速解决方法，可用于

rlmerMod

对象。基本上，您只需提取固定效应、随机效应和残差的方差分量，然后根据Nakagawa&Schielzeth（2013）提供的公式手动计算边际和条件R^2

库（robustlmm）
图书馆（洞察）
图书馆（lme4）
数据（染料，包装=“lme4”）
robust.model我自己来回答这个问题。我在R中找不到与
lmerMod
对象中的
R.squaredGLMM
等效的包或函数，但我找到了一种快速解决方法，可用于
rlmerMod
对象。基本上，您只需提取固定效应、随机效应和残差的方差分量，然后根据Nakagawa&Schielzeth（2013）提供的公式手动计算边际和条件R^2

库（robustlmm）
图书馆（洞察）
图书馆（lme4）
数据（染料，包装=“lme4”）
稳健模型
library(robustlmm)
library(insight)
library(lme4)
data(Dyestuff, package = "lme4")
robust.model <- rlmer(Yield ~ 1|Batch, data=Dyestuff)
var.fix <- get_variance_fixed(robust.model)
var.ran <- get_variance_random(robust.model)
var.res <- get_variance_residual(robust.model)
R2m = var.fix/(var.fix+var.ran+var.res)
R2c = (var.fix+var.ran)/(var.fix+var.ran+var.res)