R 출력 결과 다루기 (2)

과학적 연구 질문을 생각해 보면, IQ가 순전히 유전에 의해 결정된다면 두 쌍둥이 집합 간의 회귀선 기울기가 1일 것으로 예상할 수 있습니다. 기울기의 가설값인 1을 검정하려면, 관측된 기울기가 가설값 1에서 얼마나 벗어났는지를 평가하는 새로운 검정 통계량을 만들면 됩니다.

$$new_t = \frac{slope - 1}{SE}$$

기울기가 1이라는 가설이 참이라면, 새로운 검정 통계량은 t-분포를 따르며, 이를 이용해 p-값을 계산할 수 있습니다.

모델에서 생물학적 항은 biological_term으로 제공됩니다.

이 연습은 강의의 일부입니다

R을 활용한 선형 회귀 추론

강의 보기

연습 안내

twins 데이터셋의 자유도를 계산하세요.
참 기울기가 1과 다르다는 대립 가설에 대한 양측 p-값을 단계적으로 계산하세요.
- 기울기 estimate에서 1을 뺀 값을 표준 오차로 나눠 검정 통계량을 계산하세요.
- t-분포의 누적 분포 함수인 pt()를 사용해, 앞서 계산한 자유도와 test_statistic을 기반으로 단측 p-값을 계산하세요.
- 단측 p-값에 2를 곱하여 양측 p-값을 계산하세요.

실습형 인터랙티브 연습

이 예제를 이 샘플 코드를 완성하여 풀어보세요.

# Calculate the degrees of freedom of twins
degrees_of_freedom <- nrow(twins) - 2

biological_term %>%
  mutate(
    # Calculate the test statistic
    test_statistic = ___,
    # Calculate its one-sided p-value
    one_sided_p_value_of_test_statistic = ___,
    # ... and its two-sided p-value
    two_sided_p_value_of_test_statistic = ___
  )

코드 편집 및 실행