Hoe de t-test van welch uit te voeren in stata
De meest gebruikelijke manier om gemiddelden tussen twee onafhankelijke groepen te vergelijken, is door een t-test met twee steekproeven te gebruiken. Deze test gaat er echter van uit dat de varianties tussen de twee groepen gelijk zijn.
Als u denkt dat de variantie tussen de twee groepen niet gelijk is , kunt u de t-test van Welch gebruiken, wat het niet-parametrische equivalent is van de t-test met twee steekproeven.
In deze tutorial wordt uitgelegd hoe u de Welch’s t-test in Stata uitvoert.
Voorbeeld: Welch’s t-test in Stata
Voor dit voorbeeld gebruiken we de Fuel3- gegevensset, die de mpg bevat van 12 auto’s die een bepaalde brandstofbehandeling hebben ondergaan en van 12 auto’s die dat niet hebben gedaan.
Volg de volgende stappen om een Welch’t T-test uit te voeren om te bepalen of er een verschil is in de gemiddelde mpg tussen de twee groepen.
Stap 1: Gegevens laden en weergeven.
Laad eerst de gegevensset door de volgende opdracht in het vak Opdracht te typen:
gebruik https://www.stata-press.com/data/r13/fuel3
Bekijk de onbewerkte gegevens met behulp van de volgende opdracht:
lijst
Stap 2: Visualiseer de gegevens.
Voordat we de t-test van Welch uitvoeren, maken we eerst twee boxplots om de mpg-verdeling voor elke groep te visualiseren:
mpg grafische doos, over (verwerkt)
We kunnen zien dat de mpg van groep 1 (de groep die de brandstofbehandeling heeft ondergaan) doorgaans hoger is dan die van groep 0. We kunnen ook zien dat de variantie van groep 1 iets kleiner lijkt dan die van groep 0 (de breedte van de doos is kleiner).
Stap 3: Voer de t-test van Welch uit
Gebruik de volgende syntaxis om de Welch t-test uit te voeren:
ttest variabele_naar_maat, door (groepering_variabele) welch
Hier is de syntaxis voor ons specifieke voorbeeld:
ttest mpg, door (verwerkt) Welch
Zo interpreteert u het resultaat:
- De gemiddelde mpg voor Groep 0 was 21 . Het 95% betrouwbaarheidsinterval voor het werkelijke populatiegemiddelde was (19,26525, 22,73745) .
- De gemiddelde mpg voor Groep 1 was 22,75 . Het 95% betrouwbaarheidsinterval voor het werkelijke populatiegemiddelde was (20,68449, 24,81551) .
- Het gemiddelde verschil in mpg voor Groep 0 – Groep 1 was -1,75 . Het 95% betrouwbaarheidsinterval voor het werkelijke verschil tussen de populatiegemiddelden was (-4,28369, 0,7836902) .
- De teststatistiek, t , voor de t-test van Welch was -1,4280 .
- Omdat we geïnteresseerd zijn in de alternatieve hypothese dat de gemiddelde mpg eenvoudigweg verschillend was tussen de twee groepen, zullen we kijken naar de p-waarde geassocieerd met Ha:diff != 0, die 0,1666 blijkt te zijn. Aangezien deze waarde niet minder dan 0,05 bedraagt, hebben we niet voldoende bewijs om te zeggen dat de gemiddelde mpg tussen de twee groepen verschillend is.
Stap 4: Rapporteer de resultaten.
Ten slotte willen we graag de resultaten van onze Welch’s t-test rapporteren. Hier is een voorbeeld van hoe u dit kunt doen:
Er werd een Welch’s t-test uitgevoerd om te bepalen of er een statistisch significant verschil in mpg was tussen een groep auto’s die een brandstofbehandeling kregen en een groep die dat niet deed. De steekproefomvang voor beide groepen bedroeg 12 auto’s.
De t-test van Welch onthulde dat er geen statistisch significant verschil was in de gemiddelden (t = -1,4280, p = 0,1666) tussen de twee groepen.
Het 95% betrouwbaarheidsinterval voor het werkelijke gemiddelde verschil tussen Groep 0 (niet-behandelde groep) en Groep 1 (behandelde groep) bleek (-4,28369, 0,7836902) te zijn.