Работа с выводом R (2)

Если обратиться к научному вопросу исследования: если IQ определяется исключительно генетикой, то ожидаемый наклон прямой между двумя группами близнецов должен равняться 1. Проверить гипотезу о том, что наклон равен 1, можно с помощью новой тестовой статистики, которая оценивает, насколько наблюдаемый наклон отличается от гипотетического значения 1.

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

Если гипотеза о том, что наклон равен единице, верна, то новая тестовая статистика будет иметь t-распределение, которое можно использовать для вычисления p-значения.

Биологический член из модели доступен как biological_term.

Это упражнение является частью курса

Вывод для линейной регрессии в R

Посмотреть курс

Инструкции к упражнению

Вычислите степени свободы для набора данных twins.
Вычислите двустороннее p-значение для альтернативной гипотезы о том, что истинный наклон отличается от 1. Выполняйте вычисление поэтапно.
- Вычислите тестовую статистику как estimate минус 1, делённое на стандартную ошибку.
- Вычислите одностороннее p-значение тестовой статистики с помощью функции кумулятивного распределения t-распределения pt() при значении test_statistic и только что вычисленных степенях свободы.
- Вычислите двустороннее p-значение как удвоенное одностороннее 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 = ___
  )

Редактировать и запускать код