Получение границ участка в R

Question

Мой вопрос: Как найти границы для участка?

Я использую функцию plotHR для построения сплайнов из регрессии Кокса. Я сделал некоторые адаптации, найденные здесь , и график плотности меня раздражает, поэтому я хочу иметь возможность изменить его на простой правильный многоугольник без необходимости использования второго графика. Проблема в том, что второй график необходимо переместить, как только изменится средний размер.

Теперь я знаю, что могу использовать plot_boundaries.y и plot_boundaries.x вместе с ограничениями и добавлением информации xaxs / yaxs, но это, похоже, не самый интуитивный способ сделать это. Есть ли простой способ получить мои графики внизу слева x, y и вверху справа x, y?

Вот пример того, как использовать plotHR:

library(survival)

hmohiv<-read.table("http://www.ats.ucla.edu/stat/R/examples/asa/hmohiv.csv", sep=",", header = TRUE)

surv <- with(hmohiv, Surv(time, censor))
fit <- coxph(surv~ pspline(age), data=hmohiv)
par(xaxs="i", yaxs="i")
plotHR(fit, bty="l", ylim=c(.4, 5), y.ticks=c(.5, 1, 1.5, 2, 3, 4), xlim=c(25, 55))

Answer 1

Если вы используете pkg: выживание (или любую функцию построения графиков, которая зависит от базовой графики, как и функция plotHR, которую yopu теперь сделала доступной в пакете greg), то это должно работать:

plot(survfit(fit))
par("usr")
# [1]  0.0 62.4  0.0  1.0

«Левая нижняя» точка x-y равна par("usr")[c(1, 3)], а верхняя правая - par("usr")[c(2, 4)]. Спасибо за публикацию этой функции в пакете пакета Greg.