2012-01-08 20 views
3

質問:プロット領域の境界を見つけるにはどうすればよいですか?プロットの境界を取得するR

私はcox回帰からスプラインをプロットするのにplotHR関数を使用しています。私はいくつかのアダプテーションを行ったhereと密度プロットは私を悩ますので、私はそれを2番目のプロットを使う必要なしに単純な正多角形に変えたいと思っています。問題は、2番目のプロットでは、媒体サイズが変わるとすぐに再配置する必要があることです。

今、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)) 
+2

'?par'を参照してください、特に' par( "usr") ' –

+0

の可能性があります[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-ples-axes-ranges-for-a-ggplot2 -object)は、プロットの境界を得るのが簡単ではないと仮定するのを嫌ってしまいます。 :) – Iterator

答えて

2

あなたはPKGを使用している場合:生存、(またはyopuが現在利用可能になっているplotHR機能がないように、ベースのグラフィックスに依存して任意のプロット関数をグレッグ・パッケージに)、これは動作するはずです:

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

「を左下の」XYポイントがpar("usr")[c(1, 3)]で、右上がpar("usr")[c(2, 4)]です。その機能をGregパッケージパッケージに公開していただきありがとうございます。

+0

ありがとうございます。プロットHRは私の質問でリンクされていますが、これはスプラインをプロットする小さくても非常にきれいな方法です。少なくとも私が今までに見つけた最高のものです –

関連する問題