问题描述
是否可以使用gtsummary R程序包创建一个包含2列的事前摘要表,以汇总2个不同时间点的多个变量?
我知道arsenal R package支持这一点,但我希望尽可能使用gtsummary,因为它支持tidyverse。
例如,是否可以使用gtsummary创建类似于此example中的表的事前摘要表?这是他们示例中的数据集的简单版本:
dat <- data.frame(
tp = paste0("Time Point ",c(1,2,1,2)),id = c(1,3,4,5,6),Cat = c("A","A","B",NA,"B"),Fac = factor(c("A","C","A")),Num = c(1,NA),stringsAsFactors = FALSE)
请注意,数据集为“长格式”:tp
是2个发布前的时间点,id
是2个重复测量的主题ID。为了制作表格,Cat
和Fac
是分类变量,将在每个时间点汇总为count(%),并使用McNemar检验比较它们是否随时间变化。 Num
是一个数字变量,将在每个时间点汇总为均值(标准差),并使用配对t检验评估随时间的变化。
解决方法
简短的回答是,目前没有简单的方法。您可以使用add_stat()
函数添加所需的任何统计信息,但是您需要编写一些自定义代码
我们一直在考虑实现显示差异(配对和非配对数据)的最佳方法。有两个功能确实会增加均值差异和比率/百分比差异。尽管它们不在gtsummary软件包中,但它们位于类似gtsummary的表中,将来它们的功能可能会迁移到gtsummary中。 http://www.danieldsjoberg.com/bstfun/reference/tbl_ancova.html http://www.danieldsjoberg.com/hotfun/reference/tbl_propdiff.html
最后,我想为此添加一个功能,例如tbl_summary(...) %>% add_difference()
。但是需要制定出最佳的API。