今天群内小伙伴们介绍rempsyc包可以快速制作word文档。小编测试后发现非常不错,rempsyc包提供SCI论文类似统计表格,一般格式为APA的第7版表格(Times New Roman字体,size12,三线表水平线,双倍行距,正确的小数位数等)。
当然rempsyc包还可轻松定制各种统计图表(通过 ggplot2)可导出到Word的漂亮 APA 表格(通过 flextable)轻松运行统计测试或检查假设,以及自动执行各种其他任务。具体见原文,https://cran.r-project.org/web/packages/rempsyc/vignettes.
利用dataframe表格,自动优化为.docx格式,会自动斜体和95%CI合并!很牛。下面小编演示手动版(复杂,自定义开发度高)、自动版broom (仅支持"t.test", "lm", "cor.test", and "wilcox.test")和report (仅支持"t.test", "lm", "cor.test" )过程。
nice_table还可以很好地与rempsyc包中的其他函数集成:nice_t_test、nice_mod、nice_slopes、nice_lm和nice_lm_slopes。小编感觉意义意义不是很大。
结合前面的图表,我们尝试采用自编cox分析输出的自定义汇总表格,再利用nice_table()进行三线表绘制,很方便!
四、R代码
#rempsyc包
#install.packages("rempsyc")
library(rempsyc)
nice_table(
mtcars[1:3, ],
title = c("Table 1", "Motor Trend Car Road Tests"),
note = c(
"The data was extracted from the 1974 Motor Trend US magazine.",
"* p < .05, ** p < .01, *** p < .001"
)
)
#----案例1 lm表格 -----
# 标化数据,list
mtcars.std <- lapply(mtcars, scale)
# linear model
model <- lm(mpg ~ cyl + wt * hp, mtcars.std)
# Gather summary statistics
stats.table <- as.data.frame(summary(model)$coefficients)
CI <- confint(model)
#合并
stats.table <- cbind(row.names(stats.table), stats.table, CI)
#重命名
names(stats.table) <- c("Term", "B", "SE", "t", "p", "CI_lower", "CI_upper")
#---nice_table绘制表格!
my_table <- rempsyc::nice_table(stats.table)
#直接word临时预览
print(my_table, preview = "docx")
#导出路径.docx
flextable::save_as_docx(my_table, path = "案例1nice_tablehere.docx")
##, stars = FALSE可删除*
#rempsyc::nice_table(stats.table, highlight = .05)
#----??案例2 表格 -----
library(broom)
model <- lm(mpg ~ cyl + wt * hp, mtcars.std)
#tidy函数抓取model信息
stats.table <- broom::tidy(model, conf.int = TRUE)
#直接console中输出
#(stats.table <- broom::tidy(model, conf.int = TRUE))
##----2.1 broom----
my_table2 <- rempsyc::nice_table(stats.table, broom = "lm")
#仅支持"t.test", "lm", "cor.test", and "wilcox.test"
my_table2
##----2.2 report::report_table-----
library(report)
model <- lm(mpg ~ cyl + wt * hp, mtcars.std)
stats.table <- report::report_table(model)
#全信息版; "t.test", "lm", and "cor.test"
nice_table(stats.table, report = "lm")
#简约版
nice_table(stats.table, report = "lm", short = TRUE)
#-----3 其他附带mice_ ----
##---- 3.1 nice_t_test-----
t_test <- nice_t_test(mtcars,
response = c("mpg", "hp", "wt"),
group = "vs") # 默认为 Welch t-test
t_test
my_table1 <- nice_table(t_test)
my_table1
##---- 3.2 nice_mod 调节效应-----
mod_test <- nice_mod(
data = mtcars,
response = "mpg",
predictor = "gear",
moderator = "wt"
)
mod_test
my_table2 <- nice_table(mod_test)
my_table2
“rempsyc::nice_table 一键转换三线表,推荐。
本公众号建立了学术交流群12345,仅供SCI学术交流,人数已满,需要后备注入群,请加小编好友且实名制:姓名-学校(单位)-专业