设置使用 R 中的 coxph 函数计算的危险比的参考组

Question

当我coxph(Surv(Time, Status)~Class, data = df)在如下数据集上运行时，它总是按字母顺序设置参考组，在这种情况下，MutantA 将是参考组。有没有办法告诉它让 WT 成为参考组？

df

Time    Status  Class
3   1   WT
4   1   WT
5   1   WT
7   1   WT
7   1   WT
7   1   WT
7   1   WT
2   1   WT
2   1   WT
2   1   WT
5   1   WT
6   1   WT
7   1   WT
8   1   MutantA
9   1   MutantA
2   1   MutantA
12  1   MutantA
3   1   MutantA
4   1   MutantA
5   1   MutantA
7   1   MutantA
7   1   MutantA
7   1   MutantA
7   0   MutantA
2   1   MutantA
2   1   MutantA
2   1   MutantA
5   1   MutantB
6   0   MutantB
7   1   MutantB
8   1   MutantB
9   1   MutantB
2   1   MutantB
12  1   MutantB
3   1   MutantB
4   0   MutantB
5   1   MutantB
7   1   MutantB
7   0   MutantB
7   1   MutantB
7   0   MutantB
2   1   MutantB
2   1   MutantB
2   1   MutantB
5   1   MutantB
6   1   MutantC
7   1   MutantC
8   1   MutantC
9   1   MutantC
2   1   MutantC
12  1   MutantC
11  1   MutantC
4   1   MutantC
3   1   MutantC
6   1   MutantC
7   1   MutantC
11  1   MutantC
3   1   MutantC
2   1   MutantC
6   1   MutantC

Answer 1

是的，您只需要Class通过创建WT第一个级别来更改因子列中级别的顺序  ：

df$Class <- factor(df$Class, levels = c("WT","MutantA","MutantB","MutantC"))
# or using Gregor's comment:
df$Class = relevel(df$Class, ref = "WT")
coxph(Surv(Time, Status)~Class, data  = df)

Call:
coxph(formula = Surv(Time, Status) ~ Class, data = df)

               coef exp(coef) se(coef)     z     p
ClassMutantA -0.557     0.573    0.406 -1.37 0.171
ClassMutantB -0.688     0.503    0.395 -1.74 0.082
ClassMutantC -0.689     0.502    0.401 -1.72 0.085