Как вручную переместить легенду в riskRegression :: plotCalibration ()? - PullRequest
Новая
       44

Как вручную переместить легенду в riskRegression :: plotCalibration ()?

0 голосов
/ 08 мая 2020

Я произвел следующие riskRegression::plotCalibration. Легенда распечатывается автоматически и, к сожалению, портит графикi c. Я не могу понять, как его переместить.

Вопрос: есть идеи, как переместить легенду в нижний правый угол графика?

В настоящее время график выглядит следующим образом:

enter image description here

Написано с помощью 

library(rms)
i <- datadist(y)
options(datadist = "i")

gc <- cph(Surv(os.neck,mors)~score.group, data=y, y=TRUE, x=TRUE)

и 

library(riskRegression)
x <- Score(list(nomogram=gc), 
           Hist(os.neck,mors)~1,data=y, plots=c("cal"),
           times=60, metrics=c("auc","brier"),summary="ipa")

plotCalibration(x, cens.method = "jackknife", col=c("#2C77BF", "#E38072", "green", "yellow"))

Мои данные y

y <- structure(list(score.group = structure(c(2L, 2L, 2L, 4L, 2L, 
1L, 3L, 1L, 3L, 2L, 4L, 4L, 3L, 2L, 3L, 2L, 2L, 2L, 2L, 3L, 2L, 
2L, 2L, 2L, 2L, 2L, 3L, 3L, 2L, 2L, 4L, 1L, 2L, 3L, 2L, 2L, 2L, 
2L, 3L, 1L, 2L, 3L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 2L, 2L, 3L, 1L, 
3L, 3L, 2L, 1L, 2L, 2L, 2L, 2L, 4L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 
2L, 3L, 1L, 2L, 3L, 3L, 2L, 2L, 2L, 2L, 3L, 3L, 4L, 2L, 2L, 2L, 
4L, 3L, 3L, 4L, 1L, 3L, 5L, 2L, 3L, 2L, 2L, 2L, 2L, 2L, 3L), .Label = c("1", 
"2", "3", "4", "5"), class = "factor"), os.neck = c(9.63, 7.03, 
9.17, 10.48, 7.69, 15.18, 13.5, 16.33, 15.31, 12.09, 12.35, 22.28, 
15.77, 14.39, 10.02, 14.52, 8.44, 23.82, 5.95, 3.78, 19.32, 20.14, 
15.51, 19.78, 12.98, 32.92, 9.76, 5.65, 30.75, 2.79, 33.58, 27.53, 
27.63, 14.62, 29.17, 25.4, 18.43, 5.29, 30.75, 28.48, 14.69, 
13.14, 6.6, 26.81, 40.74, 11.63, 13.31, 10.41, 9.56, 17.51, 35.78, 
35.75, 37.62, 33.25, 36.96, 34.56, 40.05, 41.26, 24.34, 37.49, 
40.94, 24.11, 39.33, 11.24, 39.1, 19.75, 38.93, 39.36, 36.34, 
48, 29.17, 47.93, 3.68, 24.21, 46.36, 49.12, 50.96, 14.16, 54.01, 
19.88, 50.86, 1.87, 54.24, 13.93, 11.6, 10.05, 23.1, 62.78, 12.58, 
39, 59.83, 6.77, 60.39, 18.46, 61.77, 58.41, 49.45, 64.26, 2.4, 
26.51), mors = c(0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 
0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 
1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 
0L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 
0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 
0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 
1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L)), row.names = c(NA, 100L), class = "data.frame")
...