r - 在 R 中获取绘图边界

Question

我的问题:如何找到地块区域的边界？

我正在使用 plotHR用于从 cox 回归绘制样条曲线的函数。我做了一些改编发现 here密度图让我很烦恼，所以我希望能够将其更改为简单的正多边形，而无需使用第二个图。问题是，一旦中等尺寸发生变化，第二个图就需要重新定位。

现在我知道我可以将 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:survival，(或任何依赖于基础图形的绘图函数，就像 yopu 现在在 greg 包中提供的 plotHR 函数一样)，那么这应该可以工作:

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 包中发布该函数。