パッケージを使っているのですが、関数geoR
を使って得られたグラフを1つだけ表示したいです。plot
パッケージで提供されるデータセットを使用する:
library(geoR); data(elevation)
plot(elevation)
これにより、以下のように 2 x 2 グリッド上に 4 つのプロットが得られます。右下のプロットを単独で使用したいのですが、このプロットを単独で取得する方法がわかりません。
だから私は最初からそれをプロットしようとしました:
axExFact <-1.1 # to set fitting y-axis limits
ymax <- max(c(0, signif(max(elevation[[2]])*axExFact, digits=1)))
ymin <- min(c(0, signif(min(elevation[[2]])*axExFact, digits=1)))
with(elevation, hist(data, main='Plot', xlab = 'Value ',cex.lab=1.25, font.lab=2,ylim=c(ymin, ymax)))
データの範囲を完全にカバーするように y 軸の範囲を調整することはできますが、ヒストグラムに沿って密度推定値を追加することはできません。でできると思ったlines(density(as.numeric[[elevation]]))
のですが、うまくいきません。
したがって、関数で取得したグラフを取得する方がはるかに簡単plot
です。その場合、唯一の問題はy軸を調整することです。どんな提案も歓迎します。