Lavorare con l'output di R (1)
Per impostazione predefinita, il p-value fornito dall'output di lm è a due code. Nello studio sui gemelli, può essere più ragionevole adottare l'ipotesi scientifica unilaterale che i punteggi di QI dei gemelli siano associati in modo positivo. Poiché il p-value è la probabilità dei dati osservati o più estremi, il p-value del test a due code è il doppio rispetto a quello a una coda. Quindi, per ottenere un p-value a una coda dall'output a due code in R, dividi il p-value per due.
Nel codice è già fornito il modello di regressione lineare di Foster rispetto a Biological, model.
Questo esercizio fa parte del corso
Inferenza per la regressione lineare in R
Istruzioni dell'esercizio
- Estrai il coefficiente di Biological dal modello.
- Riordina (tidy) il modello.
- Filtra per il
term"Biological".
- Aggiungi una colonna,
one_sided_p_value, con il p-value a una coda.
esercizio interattivo pratico
Prova questo esercizio completando questo codice di esempio.
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 = ___)