代码之家  ›  专栏  ›  技术社区  ›  b82mo

R: 如何从coxme包中的lmekin中提取R平方

  •  0
  • b82mo  · 技术社区  · 7 年前

    我试图找出我的带有家庭效应的混合模型与数据的吻合程度。可以从lmekin函数中提取r平方值吗?如果是这样,是否有可能提取每个协变量的部分r平方值? 例子:

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

    我试过MuMin软件包,但它似乎不适用于lmekin模型。谢谢

    1 回复  |  直到 7 年前
        1
  •  0
  •   Lucy    7 年前

    我可以使用 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 包装)。