我的问题:我如何找到情节区域的边界?
我正在使用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))
答案 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包中发布该功能。