Сюжет: как изменить перерывы?

Question

Я строю кривую выживания с библиотекой survival, и я не нашел способа изменить диапазон разрывов.

Например:

library(survival)
temps <- c(5,15,25,30,18,16,38,40,40,40,40,40,40,40,40,40)
deces <- c(1,0,1,1,1,0,1,1,0,0,0,0,0,0,0,0)
plot(survfit(Surv(temps,deces)~1))    

Дает используя этот график:

Если я хочу, я могу изменить диапазон всей оси с помощью xlim=c(), шкала абсцисс с помощью xscale.

Но диапазоны х разрывов всегда будут оставаться на каждые 5 или 10 единиц (10,20,30,40). Например, невозможно изменить их на 12, как я хочу (12, 24, 36, ...), потому что они должны представлять месяцы.

Можно ли их изменить? Или мне нужно использовать библиотеку survminer?

Answer 1

Чтобы настроить оси, сохраните объект survfit и нанесите его с графическим параметром xaxt = "n", что означает, нет x ось . Затем нарисуйте ось x с помощью axis().

library(survival)

temps <- c(5,15,25,30,18,16,38,40,40,40,40,40,40,40,40,40)
deces <- c(1,0,1,1,1,0,1,1,0,0,0,0,0,0,0,0)

sv <- survfit(Surv(temps, deces) ~ 1)
plot(sv, xaxt = "n")
axis(1, at = seq(0, max(temps) + 12, by = 12))