Lucrul cu rezultatele din R (2)

Dacă IQ este determinat exclusiv de genetică, ne-am aștepta ca panta dreptei dintre cele două seturi de gemeni să fie 1. Testarea valorii ipotetice a pantei egale cu 1 se poate realiza prin construirea unei noi statistici de test, care măsoară cât de departe se află panta observată față de valoarea ipoteticată de 1.

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

Dacă ipoteza că panta este egală cu 1 este adevărată, noua statistică de test va urma o distribuție t, pe care o putem folosi pentru calculul p-valorii.

Termenul biologic din model este disponibil ca biological_term.

Acest exercițiu face parte din cursul

Inferență pentru regresia liniară în R

Vezi cursul

Instrucțiuni pentru exercițiu

Calculează gradele de libertate ale setului de date twins.
Calculează p-valoarea bilaterală pentru ipoteza alternativă că panta reală diferă de 1. Construiește calculul în etape.
- Calculează statistica de test ca estimate (panta) minus 1, totul împărțit la eroarea standard.
- Calculează p-valoarea unilaterală a statisticii de test folosind funcția de distribuție cumulativă a distribuției t, pt(), la valoarea test_statistic, cu gradele de libertate calculate anterior.
- Calculează p-valoarea bilaterală ca dublul p-valorii unilaterale.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# 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 = ___
  )

Editează și rulează codul