Werken met R-output (1)
De p-waarde die lm geeft, is standaard tweezijdig. In de tweelingenstudie is het misschien logischer om de eenzijdige wetenschappelijke hypothese te volgen dat de IQ-scores van de tweelingen positief samenhangen. Omdat de p-waarde de kans is op de geobserveerde data of extremer, is de p-waarde van de tweezijdige test twee keer zo groot als die van de eenzijdige test. Met andere woorden: om in R een eenzijdige p-waarde te krijgen uit een tweezijdige uitvoer, deel je de p-waarde door twee.
Het lineaire regressiemodel van Foster tegen Biological, model, is in het script beschikbaar.
Deze oefening maakt deel uit van de cursus
Inferentie voor lineaire regressie in R
Oefeninstructies
- Haal de coëfficiënt voor Biological uit het model.
- Maak het model netjes met tidy.
- Filter op de
term"Biological".
- Voeg een kolom toe,
one_sided_p_value, met de eenzijdige p-waarde.
Interactieve oefening met praktijkervaring
Probeer deze oefening door deze voorbeeldcode aan te vullen.
model <- lm(Foster ~ Biological, data = twins)
# Get the Biological model coefficient
biological_term <- model %>%
# Tidy the model
___ %>%
# Filter for term equal to "Biological"
___
biological_term %>%
# Add a column of one-sided p-values
___(one_sided_p_value = ___)