Cómo comparar dos muestras de datos con la prueba t de r
Si desea utilizar R de t.test () función para comparar sus datos, primero hay que comprobar, entre otras cosas, si ambas muestras se distribuyen normalmente.
Si desea saber si la temperatura media se diferencia entre los períodos del castor es activa e inactiva, puede hacerlo con un simple comando:
gt; t.test (temp ~ activ, los datos = beaver2) Welch muestra de dos 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 del grupo 0 significa en el grupo 137.0968437.90306
Normalmente, sólo se puede llevar a cabo una prueba t de muestras para las que las varianzas son aproximadamente iguales. R utiliza variación de Welch en el t-test, que corrige para varianzas desiguales.
Se obtiene una gran cantidad 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 acuerdo p-valor. 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 valor p es menor que el límite de significación. En general, los científicos consideran que la hipótesis alternativa es cierto si el valor 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 que proporciona el medio de ambas muestras.
Uno lee la fórmula temp ~ Activ como “evaluar temperatura dentro de los grupos determinados por Activ.”Alternativamente, se pueden utilizar dos vectores separados para las muestras que desea comparar y pasar tanto a la función, como en el siguiente ejemplo:
gt; activetemp lt; - beaver2 $ temp [beaver2 $ activ == 1] gt; inactivetemp lt; - beaver2 $ temp [beaver2 $ activ == 0] gt; t.test (activetemp, inactivetemp)