Робота з 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 = ___
)