R çıktısıyla çalışmak (2)
Bilimsel araştırma sorusunu düşünürken, zekâ yalnızca genetikten etkileniyorsa, ikiz çiftleri arasındaki doğrunun eğiminin 1 olmasını bekleriz. Önerilen eğim değeri 1'i test etmek, gözlenen eğimin 1'den ne kadar uzak olduğunu değerlendiren yeni bir test istatistiği oluşturarak yapılabilir.
$$new_t = \frac{slope - 1}{SE}$$
Eğim bir'e eşittir hipotezi doğruysa, yeni test istatistiği p-değeri hesaplamak için kullanabileceğimiz bir t-dağılımına sahip olacaktır.
Modeldeki biyolojik terim biological_term olarak mevcut.
Bu egzersiz, kursun bir parçasıdır
R'de Doğrusal Regresyon için Çıkarım
Egzersiz talimatları
twinsveri kümesinin serbestlik derecesini hesapla.- Gerçek eğimin 1'den farklı olduğu alternatif hipotezi için çift yönlü p-değerini, adım adım hesapla.
- Test istatistiğini, eğim
estimatedeğerinden1çıkarıp tamamını standart hataya bölerek hesapla. - Tek yönlü p-değerini, t-dağılımının kümülatif dağılım fonksiyonu
pt()ile,test_statisticnoktasında ve az önce hesapladığın serbestlik derecesiyle hesapla. - Çift yönlü p-değerini, tek yönlü p-değerinin iki katı olarak hesapla.
- Test istatistiğini, eğim
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# 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 = ___
)