Wie lassen sich zwei Datenproben mit R-T-Test vergleichen
Wenn Sie möchten, R zu verwenden t.test () Funktion Ihre Daten zu vergleichen, müssen Sie zunächst überprüfen, unter anderem, ob beide Proben normal verteilt sind.
Wenn Sie, wenn die durchschnittliche Temperatur zwischen den Perioden unterscheidet wissen wollen, dass der Biber aktiv und inaktiv ist, können Sie dies mit einem einfachen Befehl:
> T.test (Temp ~ activ, data = beaver2) Welch bei zwei Stichproben t-Testdaten: temp von activt = -18,5479, df = 80,852, p-Wert lt; 2.2E-16alternative Hypothese: wahre Unterschied in den Mitteln nicht gleich 095 Prozent Konfidenzintervall: -,8927106 -0.7197342sample Schätzungen: Mittelwert in Gruppe bedeuten 0 in Gruppe 137.0968437.90306
Normalerweise können Sie nur ein T-Test auf Proben, für die Durchführung der Varianzen etwa gleich groß sind. R verwendet Welch Veränderung auf dem t-Test, der für ungleiche Varianzen korrigiert.
Sie erhalten eine ganze Menge Informationen hier:
Die zweite Zeile gibt Ihnen die Teststatistik (t Für diesen Test), die Freiheitsgrade (df) Und die nach p-Wert. Der sehr kleine p-Wert zeigt an, daß die Mittel beider Proben signifikant unterscheiden.
Die alternative Hypothese sagt Ihnen, was Sie daraus schließen kann, wenn der p-Wert niedriger als der Grenzwert für Bedeutung ist. Im Allgemeinen betrachten die Wissenschaftler die alternative Hypothese wahr zu sein, wenn der p-Wert niedriger als 0,05 ist.
Die 95 Prozent Vertrauensintervall ist das Intervall, das den Unterschied zwischen dem Mittel mit 95 Prozent Wahrscheinlichkeit enthält, so in diesem Fall die Differenz zwischen den Mitteln wahrscheinlich zwischen 0,72 und 0,89 liegt.
Die letzte Zeile gibt Ihnen die Mittel beider Proben.
Sie lesen die Formel Temp ~ activ wie # 147 zu bewerten Temp innerhalb von Gruppen bestimmt durch activ.# 148- Alternativ können Sie zwei getrennte Vektoren für die Proben verwenden Sie die Funktion beide miteinander vergleichen wollen und übergeben, wie im folgenden Beispiel:
> activetemp lt; - beaver2 $ temp [beaver2 $ activ == 1]> inactivetemp lt; - beaver2 $ temp [beaver2 $ activ == 0]> t.test (activetemp, inactivetemp)