シミュレーションに基づく傾きの信頼区間

R による線形回帰の推測

Jo Hardin

Professor, Pomona College

置換 vs. ブートストラップのばらつき

置換データからの傾き

ブートストラップデータからの傾き

置換 vs. ブートストラップのコード

置換：

twins %>% 
  specify(Foster ~ Biological) %>% 
  hypothesize(null = "independence") %>% 
  generate(reps = 100, type = "permute") %>% 
  calculate(stat = "slope")

ブートストラップ：

twins %>% 
  specify(Foster ~ Biological) %>%  
  generate(reps = 100, type = "bootstrap") %>% 
  calculate(stat = "slope")

標本分布：無作為化 vs. ブートストラップ

置換データからの傾き

ブートストラップデータからの傾き

さあ、練習しましょう！

R による線形回帰の推測