4周期交叉生物等效性研究中gtsummary配对Wilcoxon检验报错解决

在4周期交叉设计中，每个受试者对应多个同处理的重复观测，而paired.wilcox.test要求每个配对组（Subject）仅包含一对（T/R）数据。直接使用原始多周期数据时，gtsummary无法正确识别配对结构，从而触发summary_fun未找到的错误。

1. 预处理数据：按受试者和处理分组汇总 #

先对每个受试者的每个Treatment计算TMAX的汇总值（均值、中位数均可，根据研究需求选择），确保每个Subject-Treatment组合仅保留1条记录，让数据符合配对检验的输入要求。

示例代码（以均值汇总为例）：

library(dplyr)
library(gtsummary)

# 按Subject和Treatment分组计算TMAX均值
tmax_summary <- tmax %>%
  group_by(Subject, Treatment) %>%
  summarise(TMAX = mean(TMAX), .groups = "drop")

2. 用预处理后的数据执行配对检验 #

使用处理后的数据集运行gtsummary代码，此时每个Subject对应T和R各一个值，paired.wilcox.test可正常执行：

tmax_summary %>% 
  tbl_summary(
    by = Treatment,
    include = -Subject,
    missing = 'no',
    type = list(TMAX ~ "continuous")
  ) %>% 
  add_p(
    test = list(all_continuous() ~ "paired.wilcox.test"),
    group = Subject, 
    pvalue_fun = ~ style_pvalue(.x, digits = 3)
  )

• 若需使用中位数汇总，只需将summarise中的mean(TMAX)替换为median(TMAX)。
• 该方法适用于所有多周期交叉设计的配对检验场景，核心是保证每个配对组仅包含一对对比数据，匹配配对检验的逻辑要求。

内容的提问来源于stack exchange，提问作者Roman Maksvytis