imd 统计教程 · 数据清理与基础统计
← 教程首页 返回本模块

一键自动化SCI论文格式三线表,rempsyc包

聂志强实战医学统计2024-08-25 22:03:02广东

今天群内小伙伴们介绍rempsyc包可以快速制作word文档。小编测试后发现非常不错,rempsyc包提供SCI论文类似统计表格,一般格式为APA的第7版表格(Times New Roman字体,size12,三线表水平线,双倍行距,正确的小数位数等)。

当然rempsyc包还可轻松定制各种统计图表(通过 ggplot2)可导出到Word的漂亮 APA 表格(通过 flextable)轻松运行统计测试或检查假设,以及自动执行各种其他任务。具体见原文,https://cran.r-project.org/web/packages/rempsyc/vignettes.

一、nice_table

利用dataframe表格,自动优化为.docx格式,会自动斜体和95%CI合并!很牛。下面小编演示手动版(复杂,自定义开发度高)、自动版broom (仅支持"t.test", "lm", "cor.test", and "wilcox.test")和report (仅支持"t.test", "lm", "cor.test" )过程。

二、其他nice模块

nice_table还可以很好地与rempsyc包中的其他函数集成:nice_t_test、nice_mod、nice_slopes、nice_lm和nice_lm_slopes。小编感觉意义意义不是很大。

三、任意表格nice_table

结合前面的图表,我们尝试采用自编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学术交流,人数已满,需要后备注入群,请加小编好友且实名制姓名-学校(单位)-专业