如何在 tbl_summary() 中包含 Pairwise wilcox 测试？

Question

我真的在网上找不到任何有用的东西来解决这个问题。 我有一个数据集，其中一个变量在 4 个不同的场合（CRP1、CRP4、CRP7、CRP10）进行了测量，并且我运行了成对的 Wilcoxon 测试来比较两组之间的 CRP 下降趋势。 我还制作了 ggplot 来显示显着的 p 值。

现在我想通过 tbl_summary() function显示我的p值，但似乎我无法做到这一点。 我为获取p值而运行的 Wilcoxon 测试是在我的数据集的长格式上执行的。

library(readxl)
Statistika <- read_excel("C:\\Users\\vladi\\OneDrive\\Desktop\\Statistica_copy.xlsx")
head(Statistika)

LongCRP <- Statistika %>%
  select(LEK, CRP1, CRP4, CRP7, CRP10) %>%
  filter(complete.cases(.)) %>%
  gather("Vreme", "CRP", CRP1, CRP4, CRP7, CRP10) %>%
  mutate(Vreme = factor(Vreme, levels = c("CRP1", "CRP4", "CRP7", "CRP10"))) #Long format

pairwise.wilcox.CRP <- LongCRP %>%
  group_by(LEK) %>%
  pairwise_wilcox_test(
    CRP ~ Vreme
  )

pairwise.wilcox.CRP <- pairwise.wilcox.CRP %>% add_xy_position(x = "LEK")

ggplot(LongCRP, aes(x = LEK, y = CRP, color = Vreme)) +
  geom_boxplot() +
  stat_pvalue_manual(pairwise.wilcox.CRP, label = "p.adj.signif", 
                     step.increase = 0.03,
                     bracket.nudge.y = 20) +
  theme_minimal() +
  ggtitle("Trend opadanja vrednosti CRP-a kod lečenih i nelečenih bolesnika") +
  labs(x = "Terapija") +
  scale_color_discrete(name = "Vreme\nmerenja CRP-a")

现在我知道我可以从我的 ggplot 图中提取p值，但我真的很想为此使用 tbl_summary() 。

感谢您的帮助！