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)

（我很确定这样做有效，但有一点很重要，那就是创建一个可重现的示例，这样我就可以运行你的代码进行双重检查，例如我不确定是什么{{1}看起来像，我无法按原样运行你的代码，这是一个很好的资源reprex包。）