Cómo comparar dos muestras de datos con la prueba t de r
Si desea utilizar R t.test () función para comparar los datos, primero tiene que comprobar, entre otras cosas, si ambas muestras se distribuyen normalmente.
Si quieres saber si la temperatura promedio varía entre los períodos del castor está activo e inactivo, puede hacerlo con un simple comando:
> T.test (temp ~ activ, data = beaver2) Welch Two-Sample t-testdata: temp por activt = -18.5479, df = 80.852, p-valor lt; -2.2e 16alternative hipótesis: la verdadera diferencia en las medias no es igual al intervalo de confianza del 095 por ciento: -0.8927106 estimaciones -0.7197342sample: media en el grupo de 0 significa en el grupo 137.0968437.90306
Normalmente, sólo se puede llevar a cabo una prueba t de muestras para que las varianzas son aproximadamente iguales. R utiliza la variación de Welch en la prueba de t, que corrige la desigualdad de las diferencias.
Usted recibe un montón de información aquí:
La segunda línea le da la estadística de prueba (t para esta prueba), los grados de libertad (df), Y el valor de p de acuerdo. El muy pequeño p-valor indica que los medios de ambas muestras difieren significativamente.
La hipótesis alternativa le dice lo que se puede concluir si el p-valor es inferior al límite de significación. En general, los científicos consideran que la hipótesis alternativa es cierto si el valor de p es menor que 0.05.
El intervalo de confianza del 95 por ciento es el intervalo que contiene la diferencia entre las medias con 95 por ciento de probabilidad, por lo que en este caso la diferencia entre las medias se encuentra probablemente entre 0,72 y 0,89.
La última línea le da el medio de ambas muestras.
Usted lee la fórmula temp ~ activ como " evaluar temp dentro de los grupos determinados por activ." Alternativamente, puede utilizar dos vectores separados para las muestras que desea comparar y pasar tanto a la función, como en el siguiente ejemplo:
> Activetemp lt; - beaver2 $ temp [beaver2 $ activ == 1]> inactivetemp lt; - beaver2 $ temp [beaver2 $ activ == 0]> t.test (activetemp, inactivetemp)