이상값의 영향

R을 활용한 선형 회귀 추론

Jo Hardin

Professor, Pomona College

다른 회귀선

다른 회귀 모델

starbucks_lowFib <- starbucks %>% filter(Fiber < 15)
lm(Protein ~ Fiber, data = starbucks) %>% tidy()

         term estimate std.error statistic      p.value
1 (Intercept) 7.526138 0.9924180  7.583637 1.101756e-11
2       Fiber 1.383684 0.2451395  5.644476 1.286752e-07

lm(Protein ~ Fiber, data = starbucks_lowFib) %>% tidy()

         term estimate std.error statistic      p.value
1 (Intercept) 6.537053 1.0633640  6.147521 1.292803e-08
2       Fiber 1.796844 0.2995901  5.997675 2.600224e-08

회귀 무작위화 검정 비교

전체 데이터셋

perm_slope %>% mutate(
    abs_perm_slope = abs(stat)) %>%
  summarize(
    p_value = mean(
      abs_perm_slope > abs(obs_slope)
    )
  )

 A tibble: 1 x 1
  p_value
    <dbl>
1       0

저섬유 데이터셋

perm_slope_lowFib %>% mutate(
  abs_perm_slope = abs(stat)) %>%
  summarize(
    p_value = mean(
      abs_perm_slope > abs(obs_slope_lowFib)
     )
   )

 A tibble: 1 x 1
  p_value
    <dbl>
1       0

연습해 봅시다!

R을 활용한 선형 회귀 추론