Afsnit 4.11: Teste middelværdier ens når varianser er forskellige

Vi betragter nu Statistisk Model 4.9.1 med to grupper af normalfordelte observationer med hver sin middelværdi og hver sin varians, Vi ønsker stadig at teste de to middelværdier ens, det vil sige hypotesen
Hvis middelværdierne er ens, er de to fordelinger centreret omkring det samme punkt, men den ene fordeling spreder sig ud over et større område end den anden fordeling. Alternativet, hvor middelværdierne er forskellige, skal man passe på, hvordan man fortolker. Hvis for eksempel vil man nemt danne sig det mentale billede, at observationerne fra gruppe 1 vil ligge under observationerne fra gruppe 2, men hvis samtidigt vil der være et punkt således at for vil der være større sandsynlighed i gruppe 1 end i gruppe 2 for at få en værdi over Som et konkret eksempel kan vi sige og hvor der så gælder, at hvilket er 84.1%-fraktilen i begge fordelinger. Den følgende figur med normalfordelingstætheder illustrerer eksemplet, hvor den højre del er situationen med forskellig varians.
Ligesom i konstruktionen af de to -tests tidligere i dette kapitel starter vi med en standardisering:
Da og ikke kendes, kan den standardiserede størrelse ikke bruges direkte, og i stedet erstattes med den empiriske varians
I modsætning til tidligere er denne teststørrelse ikke -fordelt på grund af de forskellige varianser. Faktisk afhænger fordelingen af stadig af de ukendte varianser og Man har dog kunnet vise matematisk, at som en approksimation kan man bruge en -fordeling til beregning af -værdi, hvor antallet af frihedsgrader i -fordelingen afhænger af data gennem og
Resultat 4.11.1. (Welch's -test)
Til at teste hypotesen om ens middelværdier, mod alternativet i model hvor hver gruppe har sin egen varians, benyttes teststørrelsen fra (4.11.1) med tilhørende observeret værdi og en approksimativ -værdi beregnes som
Endvidere er et 95%-konfidensinterval for forskel i middelværdi, det vil sige for parameteren givet ved formlen
Man kan vise, at altid er større end eller lig med det mindste af de to frihedsgradsantal og og mindre end eller lig med Metoden givet ovenfor kaldes Welch's -test.
Spørgmålet, om der findes et eksakt test for hypotesen kendes under navnet Behrens–Fisher problem.
Eksempel 4.11.2. (Menneske-maskine interaktion)
Menneske-maskine interaktion (human-computer interaction, HCI) er et vigtigt område under informationteknologi. I eksemplet her vil jeg se på tidsforbruget til at skrive en SMS. For en del år siden lavede en australsk skoleklasse et eksperiment, hvor en række personer blev bedt om at skrive en SMS med teksten The quick brown fox jumps over the lazy dog. Teksten er et pangram, det vil sige, at alle alfabetets bogstaver optræder. Data var delt op på to aldersgrupper med en tydelig forskel i tidsforbrug mellem de to aldersgrupper. Desværre ser det ud til, at data ikke længere kan findes på nettet. Jeg vil i stedet bruge data inspireret af det australske eksperiment, som jeg selv har indsamlet ved en forelæsning. Jeg brugte et dansk pangram: "Dansk jomfru på Ærø kyler halvsexet quizbog ned i wc", og delte op efter, om der blev brugt en smartphone med fuldt tastatur eller en "gammeldags" mobil, hvor hver knap koder for tre bogstaver (tripelkodende mobil).
Data i dette eksempel vedrører udelukkende kvindelige deltagere i eksperimentet. Der var 27 kvinder, der brugte en smartphone, og 33 der brugte en tripelkodende mobil. Data er givet i sekunder.
Vi lader være tidsforbrug for den 'te kvinde i den 'te gruppe (gruppe sm: smartphone, gruppe tr: tripelkodende mobil). Vi benytter modellen idet data i qqplots snor sig nogenlunde om rette linjer (kør koden i kodevinduet nedenfor). Både qqplots og boxplots peger i retning af forskellige varianser i de to populationer.
Vi ønsker med eksperimentet at undersøge, om de to teknologier er lige gode med hensyn til at skrive en SMS-besked, formuleret som hypotesen at de to middelværdier er ens -teststørrelsen for denne hypotese, når varianserne er forskellige, bliver
baseret på følgende beregnede værdier
Den tilhørende -værdi fra en -fordeling er
Da denne er meget mindre 0.05, bliver konklusionen, at data strider mod samme middelværdi, og da tyder data altså på, at det er hurtigere at skrive en besked på en smartphone fremfor en tripelkodende mobil.
Forskellen i tidsforbruget ved brug af de to mobiltyper kan angives ved et 95%-konfidensinterval for forskellen i middelværdi. Hertil finder vi 97.5%-fraktilen i en -fordeling, og konfidensintervallet bliver
Middelværdien af tidsforbruget ved brug af smartphone ligger mellem 4.13 og 14.69 sekunder under middelværdien ved brug af en tripelkodende mobil med 95% sikkkerhed.

Beregninger i R

Nedenstående kodevindue laver qqplots og boxplots af de to datasæt, og laver de beregnede værdier benyttet ovenfor.

Se opstartskoden (til/fra)

Kommenter, ud fra den dannede figur, på forholdet mellem de to varianser og forholdet mellem de to middelværdier.

Svar: To datasæt

Der ser ud til at være forskellig hældning i de to qqplots, hvilket tyder på forskellig varians i de to grupper af observationer. Den samme tendens ses i boxplottene: det højre boxplot har større udstrækning. Boxplottene indikerer også, at der er forskel i middelværdi i de to grupper: kassen i det højre boxplot ligger højere end kassen i det venstre boxplot.

ForegåendeNæste