Rのsegments関数でグラフを塗りつぶす

-
やりたいのは↓こういうことだったりします。


51. 低水準作図関数 ｜ R-Source より

やり方も上記リンク先に書いてあるので、ザッツオールです。

でも、領域を複数の台形に分けてpolygon関数で塗りつぶすという上記のやり方って、処理を頭の中で想像するのがけっこう難しい。

もっと単純に、縦棒をびっしり敷き詰めて塗りつぶすというやり方なら、頭の中で想像しやすいんじゃないかと思って、やってみました。




























































plot(dnorm, -4, 4)
xvals <- seq(2, 4, length=1000)
dvals <- dnorm(xvals)
segments( xvals,0, xvals,dvals )







最初の3行は、前述のリンク先の流用でして、塗りつぶしのところだけ、線分を描画するsegments関数を使いました。segments(x1, y1,  x2, y2) とやると、(x1,y1) から (x2,y2) まで線分を引いてくれるという関数です。引数がベクトルなら複数の線分を引いてくれます。



length=20 と少なくしてしまうと、この（あまり賢くない）やり方のボロが出てしまいます↓

























































なので、塗りつぶされている感じが出る程度まで、lengthを増やす必要があります。たくさんのグラフを描画するような場合は、このやり方は重くなってしまいそうですね。でも、グラフ1枚程度なら、1000本でも10000本でも大したことにはならないです。



塗りつぶしの方法としては理解しやすいでしょ。でも、やり方としてはpolygonで塗りつぶす方が賢いと思います。（実も蓋もない結論ですが・・・）




































コメント











コメントを投稿



















このブログの人気の投稿









Rのグラフで軸の目盛りの刻み幅を変更する方法






-














イメージ







Rでplotなどを使ってグラフを描くとき、x軸やy軸の目盛りは勝手に調整してくれて、大抵の場合はそれで問題ないのですが、たまにちょっと変えたい時があります。そのたびに必死で検索して調べているような気がするので、ここに書き留めておきます。   （最初の方の例はcurve関数を使ってますが、plot関数でも同様の方法でいけます）   もちろん、特に指定しなくても、目盛りと目盛りのラベル（例えば、x軸の -4, -2, 0, 2, 4）は書いてくれます↓    指定しなくても目盛りは表示される   # 目盛りの幅に関しては指定しない  curve(dnorm, xlim=c(-4,4))    上記では、2刻みになっていますが、例えば、1刻みにしたいなあという場合↓    目盛りが1刻みになるように指定   # -4から4までを8区間に分ける（1刻み）  curve(dnorm, xlim=c(-4, 4), xaxp=c(-4, 4, 8))    さらに刻みを細かくして0.5にしてみましたが、目盛りのラベルは適宜間引かれてしまうようです↓（作図領域を横に伸ばせばラベルが表示されるようになります）    目盛りが0.5刻みになるように指定   # 16区間に分けると0.5刻みになるが  # すべての目盛りに数値ラベルが書かれるわけではない  curve(dnorm, xlim=c(-4, 4), xaxp=c(-4, 4, 16))    いくつに分割するかちゃんと考えないと、↓こんなことになることも・・・    割り切れないところに目盛りがきてしまう例   # 8の幅を9で割るとか・・・  curve(dnorm, xlim=c(-4, 4), xaxp=c(-4, 4, 9))   次は違うやり方。   最初は軸の目盛りを書かないようにしておいて、あとから、axis関数をつかって目盛りを書き足すという方法です↓    axis関数で目盛りを後から追記したもの（見た目おんなじだけど・・・）   curve(dnorm, xlim=c(-4, 4), xaxt="n")  # 最初は目盛りを書かない  axis(side=1, at=-4:4)  # atに目盛りをベクトルで指定する   side=1というのがx軸の意味で、y軸に追記したい場合はsi








続きを読む










Rで繰り返しを含む数列の生成（rep関数、seq関数）






-















Rで、「123123123...」とか「111222333...」とか「13579...」（奇数）とか「02468...」（偶数）とか、そんな感じの繰り返しを含むような数字の列を作りたい場合、1行でさらっと書けてしまいます。   【繰り返し】  > rep(7, 10)   [1] 7 7 7 7 7 7 7 7 7 7   「○を○個だけrepeat」って感じで。    【1ずつ増加】  > 1:7   [1]  1  2  3  4  5  6  7   コロンでつなげば「○から○までインクリメント」    【○ずつ増加】  > seq(1, 7, 2)  [1]  1  3  5  7   sequence関数で、「○から○まで○ずつ増加」    【1ずつ増加する数列を繰り返す】    > rep(1:7, 3)   [1] 1 2 3 4 5 6 7 1 2 3 4 5 6 7 1 2 3 4 5 6 7   repeat関数の第1引数に、「1ずつ増加する数列」を指定すればOK。    【○ずつ増加する数列を繰り返す】  > rep(seq(1,7,2), 3)   [1] 1 3 5 7 1 3 5 7 1 3 5 7   repeat関数の第1引数に、「○ずつ増加する数列」を指定すればOK。    【それぞれ○個】    > rep(1:7, each=3)   [1] 1 1 1 2 2 2 3 3 3 4 4 4 5 5 5 6 6 6 7 7 7    【○ずつ増加、それぞれ○個】  > rep(seq(1,7,2), each=3)   [1] 1 1 1 3 3 3 5 5 5 7 7 7   コロンで指定していたところを、repにするだけなので簡単ですね。    【○ずつ増加、○回繰り返す、それぞれ○個】  > rep(seq(1,7,2), 2, each=3)   [1] 1 1 1 3 3 3 5 5 5 7 7 7 1 1 1 3 3 3 5 5 5 7 7 7   rep関数の第2引数は繰り返し回数で、「 rep(seq(1,7,2), times=2, each=3) 」と書いても同じ意味です。こう書いた方が、勘違いは少なそうですね。








続きを読む










reorderを使ってggplotの棒グラフの並び順を降順にする方法






-














イメージ







# サンプルデータの作成 fruits <- c("apple","durian","orange") count  <- c(2, 1, 3) df <- data.frame(fruits, count) df # 中身の確認   fruits count 1  apple     2 2 durian     1 3 orange     3   このデータを使って、ggplotで棒グラフを描いてみると、   library(ggplot2) ggplot(df, aes(x=fruits, y=count))    + geom_bar(stat="identity")     並び順はデータのまま    当然ながら、棒の順番はデータ通りに、２、１、３と並びますね。   これをソートしたい場合は、reorderを使います↓   ggplot(df, aes( x=reorder(fruits, count) , y=count)) +    geom_bar(stat="identity")     reorderすると昇順に並ぶ    x軸はfruitsなんだけど、reorder関数の第二引数であるcountの値で並べ替えてから使ってね、という指定です。   で、次の課題です。   reorderを使うと、昇順で１、２、３という並びになりましたが、これを降順の３、２，１という並びにしたいときは、どうやればいいか。   データフレームをソートするときなんかに使うorder関数でいうところの「decreasing=TRUE」みたいな指定があればいいのですが、reorder関数にはそのようなオプションはなさそうです。   で、実はごく簡単で、reorderの第二引数（count）の値の正負が逆になれば、順番も逆になるじゃんという理屈です。   ggplot(df, aes(x=reorder(fruits, -count ), y=count)) +    geom_bar(stat="identity")     reorder関数の第二引数にマイナスをつければ降順にできる    無事、降順に並んでくれました。








続きを読む