Я хочу показать числа риска на графике выживаемости или кумулятивной заболеваемости под осью X

Question

Под метками оси x я бы хотел отметить значения y, соответствующие этим значениям x.Лучше всего объяснить этот график: Пример

Answer 1

Если бы Фрэнк Харрелл был здесь соавтором, я бы подождал, пока он его опубликует, но поскольку его нет, вот небольшая модификация из первого примера в help(survplot) в его великолепном пакете "rms":

require(rms)
n <- 1000
set.seed(731)
age <- 50 + 12*rnorm(n)
label(age) <- "Age"
sex <- factor(sample(c('male','female'), n, TRUE))
cens <- 15*runif(n)
h <- .02*exp(.04*(age-50)+.8*(sex=='female'))
dt <- -log(runif(n))/h
label(dt) <- 'Follow-up Time'
e <- ifelse(dt <= cens,1,0)
dt <- pmin(dt, cens)
units(dt) <- "Year"
dd <- datadist(age, sex)
options(datadist='dd')
S <- Surv(dt,e)

# When age is in the model by itself and we predict at the mean age,
# approximate confidence intervals are ok

survplot(f, age=mean(age), conf.int=.95, 
               n.risk=TRUE, adj.n.risk=.75, y.n.risk=-0.15,
               xlab="")
title(main="Simulated Survival Plot Demonstrating Annotation of N-at-risk")
mtext("Follow-up Years", side=1, line=2)