×
[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
plot(dnorm, -4, 4)
xvals <- seq(-4, 4, length=1000) # 領域をx軸方向に1000個の多角形(台形)に等分割
dvals <- dnorm(xvals) # 対応するグラフの高さ
polygon(c(xvals,rev(xvals)),
c(rep(0,1000),rev(dvals)),col="skyblue") # 塗りつぶす
参考というか自分用のコピペです。元は以下URL。感謝
http://cse.naro.affrc.go.jp/takezawa/r-tips/r/51.html
xvals <- seq(-4, 4, length=1000) # 領域をx軸方向に1000個の多角形(台形)に等分割
dvals <- dnorm(xvals) # 対応するグラフの高さ
polygon(c(xvals,rev(xvals)),
c(rep(0,1000),rev(dvals)),col="skyblue") # 塗りつぶす
参考というか自分用のコピペです。元は以下URL。感謝
http://cse.naro.affrc.go.jp/takezawa/r-tips/r/51.html
PR
Comment
Trackback
Trackback URL
Comment form
