График выживания с годами вместо дней в качестве оси X - PullRequest
Новая
       9

График выживания с годами вместо дней в качестве оси X

0 голосов
/ 16 сентября 2018

Я пытаюсь построить график выживания с survplot, показывающим годы по оси x вместо дней (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 указана в днях, а не в годах. Как я могу это изменить?

...