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

Question

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

Answer 1

我自己来回答这个问题。我在 R 中找不到等同于例如的包或函数r.squaredGLMM 在 lmerMod 个对象的情况下，但我发现了一个适用于 rlmerMod 个对象的快速解决方法。基本上，您只需提取固定效应、随机效应和残差的方差分量，然后根据 Nakagawa & Schielzeth (2013) 提供的公式手动计算边际和条件 R^2。

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)

文学：
Nakagawa, S. 和 Schielzeth, H. (2013)，一种从广义线性混合效应模型中获取 R2 的通用且简单的方法。方法 Ecol Evol，4：133-142。 doi:10.1111/j.2041-210x.2012.00261.x

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

How to obtain R^2 for robust mixed effect model (rlmer command; robustlmm)?

r

robust

mixed-models