R 我想在x轴下的生存率或累积发病率图中显示有风险的数字

在x轴标记下，我想注释对应于这些x值的y值。这张图最能说明问题：

---

如果Frank Harrell是这里的撰稿人，我会等他发帖，但因为他不是，这里有一个小小的修改，来自于他伟大的“rms”软件包中的

帮助（survplot）

中的第一个示例：

require（rms）
n请发布一个可复制的代码示例。这样做的指导是，图中的x轴点比下面的摘要多。因为没有明确的摘要到x轴值的映射，所以它甚至不能表示您所请求的内容。从更高的层面来看，这似乎是一种不好的数据呈现方式。
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)