| Rの棒グラフと折れ線の描画について |
|
| 投稿者：波音 2012/01/26(Thu) 13:21 No. 16277 |
以前にRで棒グラフと折れ線を描く方法について質問させて頂いてました。
http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc039/01537.html
棒グラフに追加でpoints()を使って折れ線を書く場合、x軸の値を適切に取得して指定しないと変なところにデータポイントが打たれてしまいます。過去に頂いた回答だと
「1〜カテゴリー数の倍数になっている」
ということで、例示したデータだと
points(1:9*3-1, x,ylim=c(0,140000), pch=20, cex=2, col= color3)
となっていました。この 1:9*3-1 の部分において、 *3-1 という部分（1:9はカテゴリ数）はどうやって分かるのでしょうか？
例えばカテゴリ数が2の場合、15の場合、21の場合などはどうなるのでしょうか。。。
何か倍数部分の数値を得る方法があるのでしょうか？（それとも試行錯誤で見つけるしかないのでしょうか）
|
| Re: Rの棒グラフと折れ線の描画について |
投稿者：青木繁伸 2012/01/26(Thu) 21:15 No. 16278  |
あのときのコメントは，バーグラフの個数が 2 という特定のときについてのものでしたね。
一般的にバーグラフの個数が k のときのマークの x 座標値は i * (k+1) - k/2 ではないかと思います。
特定の k のときの x 座標は，locator 関数を使って，マークを描画すべき x 座標を得て（当然小数点の付いた数値になりますが，それは 0.5 単位の数値に近い数値になっていると思われます），バーグラフの個数との一般式を推定してプログラムに用いれば宜しいかと思います。 |
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：波音 2012/01/28(Sat) 01:24 No. 16290 |
回答ありがとうございます。
今し方、小手調べで得た情報ですが（多少、力技ですが）棒グラフに重ねて折れ線を描く場合：
#### 水準数に応じたx軸の座標を返す関数
x.value <- function(x) 1.1984*x - 0.5149
でx軸の座標値を得られるみたいです（これは棒グラフの棒が1個のとき）。
#### 出力サンプル
# 例えば0:37（38カテゴリ）でも自動的にx軸を取得して折れ線を描ける
tbl <- table(sample(0:37, 300, replace=TRUE))
barplot(tbl)
x.jiku <- x.value(1:length(names(tbl)))
points(x.jiku, tbl, col="red", type="b")
|
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：青木繁伸 2012/01/29(Sun) 08:27 No. 16294 |
graphics:::barplot.default でプログラムを見ると，space のデフォルトは 0.2 で，棒の幅は1であることがわかり，i 番目（i=1,2,...）の棒の 左辺の座標は (i-1)*1.2+0.2，右辺の座標は i*1.2 であり，中点は i*1.2-0.5 であることが分かります。
x.value <- function(x) 1.2*x - 0.5
ですね。 |
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：波音 2012/01/29(Sun) 20:28 No. 16295 |
ありがとうございます。barplotのソースコードでspaceに代入されている部分をみれば正確なプロットの間隔を確認できたわけですね。
私はlocator関数でカテゴリ数が1〜20までの棒グラフを作成して1つ1つ座標軸を取得して、回帰して定数項と係数の部分を得ていました(^_^;) |
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：波音 2012/02/04(Sat) 03:20 No. 16336 |
最終的には添付図のようなグラフを描きたいのですが。。。
barplot()やplot()で描かれる棒グラフの横軸の値を取得できたは良いものの、また新たな問題が出てきてしまいました(泣)
例えば以下のようにして
y <- as.factor(rep(c(0, 1), c(40, 50)))
plot(y)
prb <- c(0.3, 0.5)
par(new=TRUE)
x.value <- function(x) 1.2*x - 0.5
plot(x.value(1:length(prb)), prb, ylim=c(0, 1), ann=FALSE, axes=FALSE, las=1)
棒グラフに比率の折れ線を重ね描きしたいのですが、par(new=TRUE)として2回目にplot()を使ったときに横軸の位置（棒グラフの真ん中の位置）が変わってしまいます。
locator()で確認してみると、par(new=TRUE)とした時点では
0.6872666 1.8931549
となっているのですが、2回目のplot()を使った後で同じくX軸座標の値を取得すると
0.9657818 1.6186847
となってしまっています。。。
何かうまい方法はないでしょうか？
|
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：青木繁伸 2012/02/04(Sat) 08:49 No. 16337 |
図を書き足すときに par(new=TRUE) を使わないことです。
plot, barplot などの高水準描画関数の後は，lines, points など低水準描画関数を使います。また，今回の例のような場合には，第2の縦軸を用意します。第1軸の目盛りにあうようにプロット値を変換して描画し，軸目盛りには，変換前の目盛りを振ります。y <- as.factor(rep(c(0, 1), c(40, 50)))
plot(y, col=4)
prb <- c(0.3, 0.5)
#par(new=TRUE)
x.value <- function(x) 1.2*x - 0.5
points(x.value(1:length(prb)), prb*50, pch=15, col=2)
lines(x.value(1:length(prb)), prb*50, col=2)
axis(4, at=0:5*25, labels=0:5*0.5) 
|
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：青木繁伸 2012/02/04(Sat) 17:17 No. 16342 |
barplot のヘルプを見ていたら，barplot は，棒の中点（x座標）を返すと書いてありました。
ヘルプはよく読むべきだなあと改めて思い知らされました。
Value
A numeric vector (or matrix, when beside = TRUE), say mp, giving the coordinates of all the bar midpoints drawn, useful for adding to the graph.
If beside is true, use colMeans(mp) for the midpoints of each group of bars, see example. |
| Re: Rの棒グラフと折れ線の描画について |
| 投稿者：波音 2012/02/05(Sun) 22:05 No. 16357 |
毎度の返信ありがとうございます。
ホ、ホントですね。。。
今やってみたら確かにbarplot()はx軸の座標（つまり棒の真ん中の値）を返してくれました。そういえば、ヒストグラムを描くhist()も $mid として、棒の中心座標をオブジェクトとして返していますから、barplot()も同じようなことがないのかと検討をつけてみるべきでした。
まずはヘルプをよく見るということを私も再認識いたします。 |
|
