2012-01-08 5 views
3

Мой вопрос: Как найти границы участка?Получение границ графика в R

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

Теперь я знаю, что могу использовать plot_boundaries.y и plot_boundaries.x вместе с ограничениями и добавлять информацию о xaxs/yaxs, но это похоже на не самый интуитивный способ сделать это. Есть ли простой способ получить мои участки внизу слева x, y и top right 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)) 
+2

см. '? Par', особенно' parr ("usr") ' –

+0

возможный дубликат значений [bounding() рамки plot [] (http://stackoverflow.com/questions/7322301/get-plot-bounding- box-values) – joran

+0

@joran В интересах других, которые происходят в будущем, можете ли вы пояснить, почему вопрос о границах для 'plotHR()' является дубликатом вопроса на границах plot()? У меня нет понимания «plotHR», но [мой опыт работы с 'ggplot2'] (http://stackoverflow.com/questions/7705345/how-can-i-extract-plot-axes-ranges-for-a-ggplot2 -объект) заставляет меня неохотно предполагать, что границы сюжета тривиальны. :) – Iterator

ответ

2

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

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

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

+0

Спасибо. СюжетHR связан в моем вопросе, это небольшой, но очень аккуратный способ построения сплайнов. По крайней мере, это лучшее, что я нашел до сих пор –

Смежные вопросы