「統計学関連なんでもあり」の過去ログ--- 047

No.21680　グラフの塗りつぶし　　【コロン】　2015/05/24(Sun) 09:59

いつもお世話になっております。
ネットにありましたものをベースに，平均50，標準偏差10の正規分布を作成し，左側0〜20，右側80〜100を色づけしようとしているのですが，うまくいきません。ご指導いただけますでしょうか。paint_regionがネットにあったものです。

paint_region <- function( start, end, col="grey", length=50 ) {
# 多角形の頂点座標（直線部）
x <- seq( start, end, length=length )
y <- rep( 0, length )
# 多角形の頂点座標（曲線部）
revx <- rev( x )
revy <- dnorm( revx )
# 多角形のx座標，y座標（左回り）
xvals <- c( x, revx )
yvals <- c( y, revy )
# 多角形を追加
polygon( xvals, yvals, col=col )
}
curve(dnorm(x, 50, 10), from=0, to=100)
paint_region(0, 20, col="red", length=50)
paint_region(80, 100, col="blue", length=50)

No.21681　Re: グラフの塗りつぶし　　【青木繁伸】　2015/05/25(Mon) 05:22

冗長なところもありますが，
revy <- dnorm( revx, 50, 10)
とすれば所望のグラフは描けるでしょう。

No.21682　Re: グラフの塗りつぶし　　【コロン】　2015/05/25(Mon) 10:24

青木先生

早速のお返事ありがとうございます。希望のグラフが描けました！！