シミュレーション結果(散布図)の正解線重ね描き〜エラーバーを添えて〜

R

よく使うんだけど、ちゃんとメモっていなかったのでメモ。

#エラーバー描画用の関数(参考リンク参照)
add.error.bars <- function(x, y, error){
  y.upper <- y + error
  y.lower <- y - error
  w <- diff(range(x))/length(x)/5
  arrows(x, y.lower, x, y.upper, code=3, angle=90, length=w)
}
#適当なデータ
x <- 1:10
y <- (x-10)^2
y.with.error <- y + 5*rnorm(10)
#作ったデータの描画
plot(x, y.with.error)
lines(x, y)
legend(8, 60, legend=c("method A", "Correct"), lty=c(0,1), pch=c(1,-1))     
add.error.bars(x, y.with.error, 10*runif(10));

結果は以下の図のようなかんじで。

参考