以年而不是天为x轴的生存图

Question

我尝试用survplot绘制一个生存图，在x轴上显示年而不是天（dod：死亡日期/结束日期； dob：出生日期； dead：身份；开始：开始观察；暴露：暴露类别）。

原始数据集中的日期为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轴以天为单位，而不是以年为单位。我该如何更改？