Page 132 - Applied Statistics with R
P. 132

132     CHAPTER 8. INFERENCE FOR SIMPLE LINEAR REGRESSION

                                                              ̂
                                 We then repeat the process for    .
                                                              0
                                 mean(beta_0_hats) # empirical mean

                                 ## [1] 3.001147


                                 beta_0             # true mean

                                 ## [1] 3

                                 var(beta_0_hats)   # empirical variance


                                 ## [1] 0.04017924

                                 var_beta_0_hat     # true variance


                                 ## [1] 0.04

                                 hist(beta_0_hats, prob = TRUE, breaks = 25,
                                      xlab = expression(hat(beta)[0]), main = "", border = "dodgerblue")
                                 curve(dnorm(x, mean = beta_0, sd = sqrt(var_beta_0_hat)),
                                       col = "darkorange", add = TRUE, lwd = 3)







                                       2.0


                                       1.5

                                   Density  1.0



                                       0.5


                                       0.0

                                                    2.5              3.0               3.5
                                                                     ^
                                                                     β 0
   127   128   129   130   131   132   133   134   135   136   137