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

Я строю кривую выживания с помощью библиотеки 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.

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

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


r range survival break
person Kolia51    schedule 20.01.2020    source источник

Ответы (1)


arrow_upward
3
arrow_downward

Чтобы настроить оси, сохраните объект 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))

введите здесь описание изображения

person Rui Barradas    schedule 20.01.2020
comment
Большое спасибо! Я никогда не привык так делать, теперь я научился! - person Kolia51; 20.01.2020