Я хочу проверить различные ограничения в анализе выживания. Вот мой код:
for (cutoff in seq(0.001,0.1,by=0.001)){
os$rna.nor.high= ifelse(os$rna.nor > cutoff, 1, 0)}
fit = survfit(Surv(os$OS,os$os.ce)~os$rna.nor.high)
ggsurvplot(fit, data= os)
Однако, он вернул только одну картинку и только протестировал cutoff=0.1
Я думаю, что мне нужно добавить один столбец (например, os$rna.nor.high
1-99) в мой фрейм данных os каждый раз, когда он проверяет новую обрезку, но я не знаю, как это сделать ....