Trabajar con salida de R (1)
El valor p que da la salida de lm es, por defecto, de dos colas. En el estudio de gemelos, puede tener más sentido seguir la hipótesis científica unilateral de que las puntuaciones de CI de los gemelos están asociadas de forma positiva. Como el valor p es la probabilidad de observar los datos o algo más extremo, el valor p de una prueba bilateral es el doble que el resultado unilateral. Es decir, para obtener un valor p unilateral a partir de la salida bilateral en R, divide el valor p entre dos.
En el script se proporciona el modelo de regresión lineal de Foster frente a Biological, model.
Este ejercicio forma parte del curso
Inferencia para la regresión lineal en R
Instrucciones del ejercicio
- Obtén el coeficiente de Biological a partir del modelo.
- Ordena el modelo con
tidy. - Filtra por el
term"Biological".
- Ordena el modelo con
- Añade una columna,
one_sided_p_value, con el valor p unilateral.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
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 = ___)