gtsummary R程序包:带有两个2-sample测试配对的pre-post摘要表?

时间:2020-08-12 20:54:53

标签: r gtsummary

是否可以使用gtsummary R程序包创建一个包含2列的事前摘要表,以汇总2个不同时间点的多个变量?

我知道arsenal R package支持这一点,但我希望尽可能使用gtsummary,因为它支持tidyverse。

例如,是否可以使用gtsummary创建类似于此example中的表的事前摘要表?这是他们示例中的数据集的简单版本:

 dat <- data.frame(
  tp = paste0("Time Point ", c(1, 2, 1, 2, 1, 2, 1, 2, 1, 2)),
  id = c(1, 1, 2, 2, 3, 3, 4, 4, 5, 6),
  Cat = c("A", "A", "A", "B", "B", "B", "B", "A", NA, "B"),
  Fac = factor(c("A", "B", "C", "A", "B", "C", "A", "B", "C", "A")),
  Num = c(1, 2, 3, 4, 4, 3, 3, 4, 0, NA),
  stringsAsFactors = FALSE)

请注意,数据集为“长格式”:tp是2个发布前的时间点,id是2个重复测量的主题ID。为了制作表格,CatFac是分类变量,将在每个时间点汇总为count(%),并使用McNemar检验比较它们是否随时间变化。 Num是一个数字变量,将在每个时间点汇总为均值(标准差),并使用配对t检验评估随时间的变化。

1 个答案:

答案 0 :(得分:0)

简短的回答是,目前没有简单的方法。您可以使用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。

但是,我不了解此示例表中报告的差异,因此我没有计划使用类似的表。 enter image description here