我有一个coxph
函数,该函数迭代data.frame并输出data.frame中每一行的生存期。我想在coxph中加入分位数功能,以便coxph根据分位数定义组。我知道如何使用Surv
函数而不是coxph来做到这一点。
coxph <- capture.output(for(i in colnames(df)){
print(summary(coxph(as.formula(paste0("Surv(Time, Status)~", i,"+Age", "+ Sex")), data=as.data.frame(df))))
})
在上述类似于以下功能的函数中将分位数添加到i
:
survfit(Surv(Time, Status) ~ paste(cut(hsa_miR_320a, quantile(hsa_miR_320a, (0:2)/2), include.lowest=T)), data = as.data.frame((df)))