У меня есть функция 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)))