Trabajar con la salida de R (2)

Si pensamos en la pregunta de investigación científica, si el CI está determinado solo por la genética, esperaríamos que la pendiente de la recta entre los dos gemelos fuese 1. Comprobar la pendiente hipotética igual a 1 puede hacerse creando una nueva estadística de prueba que evalúe cuán lejos está la pendiente observada del valor hipotetizado de 1.

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

Si la hipótesis de que la pendiente es uno es cierta, entonces la nueva estadística de prueba seguirá una distribución t que podemos usar para calcular un p-valor.

El término biológico del modelo está disponible como biological_term.

Este ejercicio forma parte del curso

Inferencia para la regresión lineal en R

Ver curso

Instrucciones del ejercicio

Calcula los grados de libertad del conjunto de datos twins.
Calcula el p-valor bidireccional para la hipótesis alternativa de que la pendiente verdadera es distinta de 1. Ve construyendo el cálculo por etapas.
- Calcula la estadística de prueba como el estimate de la pendiente menos 1, todo dividido por el error estándar.
- Calcula el p-valor unidireccional de la estadística de prueba usando la función de distribución acumulada de la distribución t, pt(), en test_statistic, con los grados de libertad que acabas de calcular.
- Calcula el p-valor bidireccional como el doble del p-valor unidireccional.

Ejercicio interactivo práctico

Prueba este ejercicio y completa el código de muestra.

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

Editar y ejecutar código