在R中获得地块边界

时间:2012-01-08 16:52:31

标签: r

我的问题:我如何找到情节区域的边界?

我正在使用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))

1 个答案:

答案 0 :(得分:2)

如果您正在使用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包中发布该功能。

相关问题