確信区間

確信区間を塗りつぶす

ベータ分布で、エラー率の確信区間を示すグラフを作る。

shape1 = 1; shape2 = 39;
d = tibble(x = rbeta(1e4, shape1, shape2), y = dbeta(x, shape1, shape2))
ggplot(data = d) + aes(x=x, y=y) + 
 geom_line() +
 geom_ribbon(data = filter(d, pbeta(x, shape1, shape2) < .90 ), aes(x = x,  ymax=y, ymin=0), fill = "blue", alpha = 0.5)