Addnig текст ограниченных кубических сплайнов участок

Question

I имеют следующие гипотетические данные для построения графика ограничена кубических сплайнов

library(rms) 
library(survival) 
n <- 1000 
set.seed(731) 
age <- 50 + 12*rnorm(n) 
label(age) <- "Age" 
sex <- factor(sample(c('Male','Female'), n, rep=TRUE, prob=c(.6, .4))) 
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) 

f <- cph(S ~ rcs(age,4) + sex, x=TRUE, y=TRUE) 
cox.zph(f, "rank")    # tests of PH 
anova(f) 
plot(Predict(f, age, sex)) # plot age effect, 2 curves for 2 sexes 

Wit вышеприведенные данные, я получаю этот прекрасный сюжет:

Я хотел бы добавить следующий текст к нему text(60,-1,"p < 0.001"), но он не работает и дает мне следующее сообщение об ошибке: plot.new has not been called yet. Какие мысли могут быть проблемой? Любой другой способ добавить текст (или легенду) к нему?

Answer 1

Пакет rms переместился на платформу построения решетки. Вам нужно будет добавить текст в панельный вызов или одну из дополнительных функций решетки. См. ?llines или используйте trellis.focus.

[![ap <- function(...) { ltext(60, 1, "p < 0.001") } 

plot(Predict(f, age, sex), addpanel=ap)][1]][1]