如何将metafor包中anova.rma的输出保存为CSV或数据框?
解决metafor包anova.rma结果导出为CSV的问题
当使用metafor包的anova.rma()比较多水平模型时,控制台会输出规整表格,但返回的对象是列表结构,直接导出会格式混乱。以下两种方法可以将结果整理为整洁格式并导出为CSV:
方法1:手动提取列表中的表格数据
anova.rma()返回的列表里,$table元素就是控制台显示的核心结果表格,直接将其转为数据框后调整格式即可:
# 假设你的anova结果存储在anova_res对象中 # 提取表格并转为数据框 anova_df <- as.data.frame(anova_res$table) # 添加模型名称列(从行名转换) anova_df$模型 <- rownames(anova_df) # 调整列顺序,让模型名称在前 anova_df <- anova_df[, c("模型", "df.resid", "logLik", "AIC", "BIC", "AICc", "LRT", "df", "pval")] # 导出为CSV write.csv(anova_df, "anova_rma结果.csv", row.names = FALSE, fileEncoding = "UTF-8")
方法2:用broom包自动整理为整洁数据框
broom包专门用于将统计模型结果转换为整洁的数据框,对anova.rma对象支持良好:
library(broom) # 整理结果 anova_tidy <- tidy(anova_res) # 导出CSV write.csv(anova_tidy, "anova_rma整洁结果.csv", row.names = FALSE, fileEncoding = "UTF-8")
如果遇到broom兼容性问题,优先用方法1手动提取,确保核心统计量(如LRT、p值、信息准则)完整保留。
内容的提问来源于stack exchange,提问作者jmmakoshark
相关产品推荐
相关产品推荐

