R:如何从 coxme 包中的 lmekin 中提取 r 平方
R: how to extract r squared from lmekin in coxme package
我正在尝试了解我的具有家庭效应的混合模型与数据的拟合程度。是否可以从 lmekin 函数中提取 r 平方值?如果是这样,是否可以为每个协变量提取部分 r 平方值?
示例:
model= lmekin(formula = height ~ score + sex + age + (1 | IID), data = phenotype_df, varlist = kinship_matrix)
我试过 MuMin 包,但它似乎不适用于 lmekin 模型。谢谢
我可以使用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 包是一个很好的资源。
我正在尝试了解我的具有家庭效应的混合模型与数据的拟合程度。是否可以从 lmekin 函数中提取 r 平方值?如果是这样,是否可以为每个协变量提取部分 r 平方值? 示例:
model= lmekin(formula = height ~ score + sex + age + (1 | IID), data = phenotype_df, varlist = kinship_matrix)
我试过 MuMin 包,但它似乎不适用于 lmekin 模型。谢谢
我可以使用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 包是一个很好的资源。