Texreg导出lmer模型结果到Latex时不显示R平方与调整R平方问题咨询

问题原因 #

texreg本身没有为lmer()生成的线性混合效应模型内置边际R²、条件R²的计算和输出逻辑，你看到的其他可输出该指标的示例，均是手动添加自定义拟合优度项实现的。而sjPlot::tab_model()默认集成了混合效应模型R²的计算逻辑，所以可以直接输出对应指标。

解决方法 #

你可以手动计算两个R²值后，通过texreg的custom.gof参数注入到输出结果中，具体操作如下：

1. 先安装加载依赖包MuMIn用于计算混合效应模型的R²

# 未安装包时先执行安装
install.packages("MuMIn")
# 加载所需包
library(MuMIn)
library(texreg)

2. 计算模型的边际R²和条件R²

r2_vals <- r.squaredGLMM(mixed.lmer.str.df)
# 提取两个R²值
marginal_r2 <- r2_vals[1, "R2m"]
conditional_r2 <- r2_vals[1, "R2c"]

3. 调用texreg时传入自定义拟合优度参数

texreg(mixed.lmer.str.df,
       booktabs = TRUE, 
       include.aic = TRUE, 
       include.bic = TRUE, 
       include.dic = TRUE,
       include.deviance = TRUE, 
       include.loglik = TRUE,
       include.nobs = TRUE,
       include.groups = TRUE,
       include.variance = TRUE,
       # 自定义添加R²指标行
       custom.gof.rows = list(
         "边际R²" = round(marginal_r2, 3),
         "条件R²" = round(conditional_r2, 3)
       ))

如果需要输出调整R²，你可以自行计算对应值后，按照同样的逻辑加入到custom.gof.rows列表中即可。

内容的提问来源于stack exchange，提问作者Sharareh Ghanbari Azar