小编Nel*_*lly的帖子

当我尝试使用 add_p() 函数获取 by 变量（有 10 个级别）和有两个级别的分类变量（是/否）之间差异的 p 值时，出现以下错误。我不确定如何提供可重现的示例。根据试验数据，我想象我的变量将是具有 10 个级别的“T 阶段”变量，分类变量将是：(1)“化疗治疗”具有 2 个级别，(2)“化疗治疗 2”具有 4 个级别水平。但这是我运行的代码。

library(gtsummary)
library(tidyverse)
miro_def %>% 
  select(mheim, age_dx, time_t1d_yrs, gender, collard, fhist_pandz) %>% 
  tbl_summary(by = mheim, missing = "no",
              type = list(c(gender, collard, fhist_pandz, mheim) ~ "categorical"),
              label = list(gender ~ "Gender", 
                           fhist_pandz ~ "Family history of PD", 
                           age_dx ~ "Age at diagnosis", 
                           time_t1d_yrs ~ "Follow-up(years)")) %>% 
  add_p() %>% 
  # style the output with custom header 
  #modify_header(stat_by = "{level}") %>% 
  # convert to …