Comment comparer deux échantillons de données avec le test t de r

Si vous souhaitez utiliser les R t.test () fonction de comparer vos données, vous devez d'abord vérifier, entre autres choses, si les deux échantillons sont distribués normalement.

Si vous voulez savoir si la température moyenne diffère entre les périodes Le castor est actif et inactif, vous pouvez le faire avec une simple commande:

> T.test (temp ~ activ, data = beaver2) Welch deux échantillons t-testdata: temp par activt = -18,5479, df = 80,852, p-valeur lt; 2.2e-16alternative hypothèse: vraie différence dans les moyens est pas égal à un intervalle de confiance de 095 pour cent: -0.8927106 estimations -0.7197342sample: moyen dans le groupe 0 signifie dans le groupe 137.0968437.90306

Normalement, vous ne pouvez effectuer un test t sur des échantillons pour lesquels les écarts sont à peu près égale. R utilise la variation de Welch sur le t-test, qui corrige des variances inégales.




Vous obtenez tout un tas d'informations ici:

  • La deuxième ligne vous donne la statistique de test (t pour ce test), les degrés de liberté (df), Et la p-value en fonction. La très faible valeur p indique que le moyen de deux échantillons diffèrent de façon significative.

  • L'hypothèse alternative vous dit ce que vous pouvez conclure si la p-valeur est inférieure à la limite de signification. Généralement, les scientifiques considèrent l'hypothèse alternative pour être vrai si la p-valeur est inférieure à 0,05.

  • L'intervalle de confiance de 95 pour cent est l'intervalle qui contient la différence entre les moyens avec 95 pour cent de probabilité, alors dans ce cas, la différence entre les moyens se situe probablement entre 0,72 et 0,89.

  • La dernière ligne vous donne les moyens des deux échantillons.

Vous avez lu la formule Temp ~ activ comme “ évaluer Temp au sein des groupes déterminés par activ.” Alternativement, vous pouvez utiliser deux vecteurs séparés pour les échantillons que vous voulez comparer et passer à la fois à la fonction, comme dans l'exemple suivant:

> Activetemp lt; - beaver2 $ temp [beaver2 $ activ == 1]> inactivetemp lt; - beaver2 $ temp [beaver2 $ activ == 0]> t.test (activetemp, inactivetemp)

» » » » Comment comparer deux échantillons de données avec le test t de r