個々の応答に対する予測区間
応答の期待値に対する区間推定に加えて、実際の個々の応答に対する区間推定が求められることも多いです。予測の枠組み自体は同じですが、予測点は回帰直線のまわりでより散らばるため、標準誤差は大きく計算されます。
期待される平均値の区間と同様に、個々の予測値の区間も、回帰直線の推定が中心付近でより安定するため、中央では小さく端へ行くほど大きくなります。平均応答の区間は、個々の応答の区間よりもかなり小さい点に注意してください。
モデルから係数レベルの情報を取り出す tidy()、観測レベルの情報を得る augment() をすでに見てきました。glance() はこの三位一体を完成させ、モデルレベルの情報を返します。
線形回帰モデルは model として与えられ、前の演習での予測は predictions に入っています。
この演習はコースの一部です
R による線形回帰の推測
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
twins_sigma <- model %>%
# Get model-level information
___ %>%
# Pull out sigma
___
predictions %>%
# Calculate the std err of the predictions
mutate(std_err_of_predictions = ___)