Práce s výstupem R (2)

Pokud je IQ podmíněno pouze genetikou, pak bychom očekávali, že sklon přímky mezi dvěma skupinami dvojčat bude roven 1. Testování hypotézy o hodnotě sklonu 1 provedeme tak, že vytvoříme novou testovou statistiku, která měří, jak daleko je pozorovaný sklon od hypotetické hodnoty 1.

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

Pokud hypotéza o sklonu rovném jedné platí, bude mít nová testová statistika t-rozdělení, které můžeme využít k výpočtu p-hodnoty.

Biologický člen z modelu je dostupný jako biological_term.

Toto cvičení je součástí kurzu

Inference for Linear Regression in R

Zobrazit kurz

Pokyny k cvičení

Vypočítej stupně volnosti datasetu twins.
Vypočítej oboustrannou p-hodnotu pro alternativní hypotézu, že skutečný sklon se liší od 1. Postup rozděl do jednotlivých kroků.
- Vypočítej testovou statistiku jako odhad sklonu estimate minus 1, celé vydělené standardní chybou.
- Vypočítej jednostrannou p-hodnotu testové statistiky pomocí funkce kumulativního rozdělení t-distribuce, pt(), v bodě test_statistic s právě vypočítanými stupni volnosti.
- Vypočítej oboustrannou p-hodnotu jako dvojnásobek jednostranné p-hodnoty.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

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

Upravit a spustit kód