График выживания с годами вместо дней в качестве оси X
Я пытаюсь построить сюжет выживания с survplot
показывая годы на оси х вместо дней (dod: дата смерти / конец f'up; dob: дата рождения; мертвый: статус; начало: начало наблюдения; экспозиция: категория экспозиции).
Даты в исходном наборе данных имеют формат даты R - преобразуются в дни as.numeric
,
library(rms)
n <- 1000
sex <- as.factor(round(runif(n,0,1),0))
dob <- as.Date(round(runif(n,-9000,-4000)), origin="1970-01-01")
dead <- as.factor(rnorm(n,0.5,1)>1)
start <- as.Date(round(runif(n,3000,4000)), origin="1970-01-01")
dod <- rep(as.Date(15000, origin="1970-01-01"),n)
df <- data.frame(cbind(sex,dob,start,dod,dead))
df$dod[which(df$dead==T)] <- as.Date(round(runif(length(which(df$dead==T)), 5000,15000)), origin="1970-01-01")
df$surv <- Surv(as.numeric(df$start), as.numeric(df$dod), origin=as.numeric(df$dob), df$dead)
ddist <- datadist(df); options(datadist="ddist")
f <- cph(surv ~ sex, data=df, x=T, y=T, surv=T)
survplot(f, sex)
Это хорошо работает, но ось X указана в днях, а не в годах. Как я могу это изменить?