gtsummary
包提供了一种优雅而灵活的方式来创建可供发表的分析表和汇总表。gtsummary
包可以总结数据集、回归模型结果等。总结数据集即对数据进行描述统计,与我们之前介绍的table1
包相似。详见:R实战 | 文章第一表:三线表的绘制。今天我们主要介绍回归模型结果表达绘制。
![efa2857a7b4a042c7d0dd8cd292cde61.png](https://img-blog.csdnimg.cn/img_convert/efa2857a7b4a042c7d0dd8cd292cde61.png)
绘制
# 安装并加载包
install.packages("gtsummary")
library(gtsummary)
# 示例数据
head(trial)
> head(trial)
# A tibble: 6 × 8
trt age marker stage grade response death ttdeath
<chr> <dbl> <dbl> <fct> <fct> <int> <int> <dbl>
1 Drug A 23 0.16 T1 II 0 0 24
2 Drug B 9 1.11 T2 I 1 0 24
3 Drug A 31 0.277 T1 II 0 0 24
4 Drug A NA 2.07 T3 III 1 1 17.6
5 Drug A 51 2.77 T4 III 1 1 16.4
6 Drug B 39 0.613 T4 I 0 1 15.6
# 建模
mod1 <- glm(response ~ trt + age + grade, trial, family = binomial)
# tbl_regression()绘制表格
t1 <- tbl_regression(mod1, exponentiate = TRUE)
t1
![78b7b8543965bff93f1e94a2ca0c313f.png](https://img-blog.csdnimg.cn/img_convert/78b7b8543965bff93f1e94a2ca0c313f.png)
# 合并多种回归结果表
library(survival)
# 生存模型
t2 <-
coxph(Surv(ttdeath, death) ~ trt + grade + age, trial) %>%
tbl_regression(exponentiate = TRUE) # 这里用了管道符简化
t2
# merge tables
tbl_merge_ex1 <-
tbl_merge(
tbls = list(t1, t2),
tab_spanner = c("**Tumor Response**", "**Time to Death**")
)
tbl_merge_ex1
![af5706f0a8cd49a288ab5c72fa506395.png](https://img-blog.csdnimg.cn/img_convert/af5706f0a8cd49a288ab5c72fa506395.png)
导出结果
# 保存为.html .tex .ltx .rtf
tbl_merge_ex1 %>%
as_gt() %>%
gt::gtsave(filename = "tbl_merge_ex1.html") # use extensions .html .tex .ltx .rtf
# 保存为word
install.packages('gdtools')
install.packages('flextable')
tf <- tempfile(fileext = ".docx")
tbl_merge_ex1 %>%
as_flex_table() %>%
flextable::save_as_docx(path = tf)
更多详细参数设置:Tutorial: tbl_summary • gtsummary (danieldsjoberg.com)(https://www.danieldsjoberg.com/gtsummary/articles/tbl_summary.html)
示例数据和代码领取
点赞
、在看
本文,分享至朋友圈集赞20个
并保留30分钟
,截图发至微信mzbj0002
领取。
木舟笔记2022年度VIP可免费领取。
木舟笔记2022年度VIP企划
权益:
2022年度木舟笔记所有推文示例数据及代码(在VIP群里实时更新)。
资源合集 木舟笔记科研交流群。
半价购买
跟着Cell学作图系列合集
(免费教程+代码领取)|跟着Cell学作图系列合集。
收费:
99¥/人。可添加微信:mzbj0002
转账,或直接在文末打赏。
![ab04d45f6f4976f3d26600f6ffec9b80.png](https://img-blog.csdnimg.cn/img_convert/ab04d45f6f4976f3d26600f6ffec9b80.png)
参考
Presentation-Ready Data Summary and Analytic Result Tables • gtsummary (danieldsjoberg.com)
往期内容
![97ac4ab0bfd62bc2c1b907590f440fdd.png](https://img-blog.csdnimg.cn/img_convert/97ac4ab0bfd62bc2c1b907590f440fdd.png)