r - R：如何从 coxme 包中的 lmekin 中提取 r 平方

Question

我试图找出我的具有家庭效应的混合模型与数据的拟合程度。是否可以从 lmekin 函数中提取 r 平方值？如果是这样，是否可以为每个协变量提取部分 r 平方值？例子：

model= lmekin(formula = height ~ score + sex + age + (1 | IID), data = phenotype_df, varlist = kinship_matrix)

我已经尝试过 MuMin 包，但它似乎不适用于 lmekin 模型。谢谢。

Answer 1

我可以使用该r.squaredLR()功能，

library(coxme)
library(MuMIn)
data(ergoStool, package="nlme")  # use a data set from nlme 
fit1 <- lmekin(effort ~ Type + (1|Subject), data=ergoStool)
r.squaredLR(fit1)

（我很确定这行得通，但最好的一件事是创建一个可重现的示例，以便我可以运行您的代码进行仔细检查，例如，我不确定是什么phenotype_df样的，我无法运行你的代码是这样的，reprex包是一个很好的资源）。