信頼区間（CI）を構築する

リサンプリングによって p-hat がどのように変動するかを1例見ましたが、ばらつきを適切に見積もるには、これを何度も繰り返す必要があります。ここでは、信頼区間の構築に用いる標準誤差（SE）を推定するために、完全なブートストラップ分布を計算します。さらに、infer の calculate() 動詞を使って、多くのデータセットから多くの統計量を計算する処理を簡潔にします。

まずは calculate の出力を確認しましょう。この関数は、データフレームを「stat」列と、それに対応する「replicate」列の2列だけに要約します。

ブートストラップ分布をプロットすると、釣鐘型の形状になっているはずです。この形状のおかげで、SE を2つ分、足し引きすることで95%区間を得られます。