r - 按疾病分层的每列的 P 值

Question

我有以下品种的数据集；真实数据集中有> 40个并发症。

Disease <- c(1,1,1,1,0,0,0,0)
Complication1 <- c(1,1,1,1,0,0,0,0)
Complication2 <- c(0,0,0,0,1,1,1,1)
dataset <- data.frame(cbind(Complication1, Complication2, Disease))


 Complication1      Complcation2      Disease
    0                  1                 1
    0                  1                 1
    0                  1                 1
    0                  1                 1
    1                  0                 0
    1                  0                 0
    1                  0                 0
    1                  0                 0

我正在尝试根据疾病状态是否存在来计算每种并发症complication1的p 值。complication2

我使用了 TableOne，但这种方法存在一些问题，我想知道如何使用循环方法对每一列执行此操作。

运行 TableOne 包

table1 <- CreateTableOne(vars = xvars, strata = “Disease", data = dataset, test = TRUE)

xvars <- c(“Complciation1", “Complciation2”)

print(table1, smd = T)

对于这个例子，输出看起来像这样

                       Stratified by Disease
                            0           1             p    
  n                         4           4                          
  Complication1 (mean (SD)) 0.25 (0.00) 0.75 (0.00)   < 0.03     
  Complication2 (mean (SD)) 0.75 (0.00) 0.25 (0.00)   < 0.03   
  

标签： r