使用 gtsummary 为变量的子类别生成 p 值

Question

我正在使用 gtsummary 生成特定变量/子类别的分层表。但是，我遇到了生成 p 值的问题。如此处所示，种族有一个 p 值，但我更希望子类别的 p 值与种族相关，并且希望种族变量标题与子类别分开。

我试过代码：

data %>% tbl_summary(by = Gender) %>% add_overall %>% add_p()

但似乎出现错误。本质上，我需要修复此问题，以便将“种族”标题与种族的子类别分开（因此左侧的新列带有种族标题，然后右侧的列将具有子类别及其 p 值）。有人对如何解决此问题有任何建议吗？任何建议表示赞赏！我附上了我需要的输出图片。

我在这里附上了一些示例数据的快照：

我尝试先选择我需要的列（上图是部分数据的快照），将列重命名为我想要的格式，并尝试使用以下代码创建表：

new_data = select(data, O_Gender, O_Age, O_Race)

names(new_data)[1] = "Gender"
names(new_data)[1] = "Age"
names(new_data)[1] = "Race"

new_table <- new_data %>% tbl_summary(by = Gender) %>% add_overall() %>% add_p()

手头的问题是我得到了年龄的 p 值，但我需要种族中包含的值的 p 值。我不确定最佳路线，因为我知道种族是变量，其他是子类别，但我需要子类别的 p 值。