Afsnit 1.2: Binomialmodellen

Binomialmodellen beskriver situationen, hvor der laves uafhængige delforsøg, og i hvert af disse er der to mulige udfald. Til sidst tælles der op, hvor mange af de delforsøg der giver et bestemt af de to udfald. Hvis sandsynligheden for dette udfald i det enkelte delforsøg er og antallet af denne type udfald betegnes med den stokastiske variabel skriver vi Generelt kan man sige, at binomialmodellen beskriver situationen, hvor vi kaster en skæv mønt gange, hvor sandsynligheden for krone er og tæller op, hvor mange gange vi fik krone. Udover at er sandsynligheden i det enkelte kast, kan også tolkes som andelen af krone i uendelig mange kast af mønten, eller som andelen af krone i "populationen" af møntkast. I binomialfordelingen, kalder vi for antalsværdien og for sandsynlighedsparameteren. Et vigtigt element i Introducerende statistik og datananalyse med MATLAB er at vænne jer til at beskrive et eksperiment gennem en statistisk model. En statistisk model involverer en eller flere stokastiske variable, og det er vigtigt at I vænner jer til dette begreb. En stokastisk variabel repræsenterer udfaldet af et eksperiment før dette gennemføres, og vil give en ny værdi hver gang eksperimentet gentages.
Definition 1.2.1. (Fordelingsnotation og statistisk model)
Hvis er en stokastisk variabel, angiver vi fordelingen ved at skrive Stokastiske variable angives som hovedregel ved store bogstaver. Når der er indsamlet data, bruges et lille bogstav til at angive den målte værdi.
I en statistisk model til beskrivelse af et eksperiment angives de stokastiske variable, der måles, og man angiver deres fordelinger. Disse fordelinger vil typisk indeholde ukendte parametre, som man netop ønsker at sige noget om gennem eksperimentet. Variationsområdet for parametrene er en del af den statistiske model.
For at kunne referere til binomialmodellen kommer her en formel angivelse.
Statistisk Model 1.2.2. (Binomialmodellen)
Lad den stokastiske variabel angive antallet ud af delforsøg med et bestemt udfald, hvor der er sandynlighed for dette udfald. Binomialmodellen til beskrivelse af dette skrives som
For binomialmodellen kender I fra jeres calculuskursus følgende resultater:
Vi kan nu vende tilbage til Kipping og Popes krystaleksperiment fra afsnit 1.1 og beskrive situationen således, at vi har observationen fra en stokastisk variabel Hypotesen om en bestemt andel af D-krystaller er på formen hvor er den værdi, vi ønsker at undersøge. Når man tester en hypotese, bør det også overvejes, hvad er alternativet, hvis hypotesen ikke er sand. I de fleste tilfælde i disse noter vil alternativet blot være alle mulige andre værdier af parameteren, her I nogle situationer kan man også være interesseret i ensidede alternativer som for eksempel
Jeg formulerer nu generelt test i binomialfordelingen og tilhørende -værdi.
Resultat 1.2.3. (-værdi for test af andel)
Betragt den statistiske model og hypotesen mod alternativet Som teststørrelse bruger vi og -værdien for en observation med tilhørende teststørrelse er

1.2.4 Beregning i python af p-værdi for test af andel

I python beregnes sandsynligheder og fordelingsfunktionen i tilfældet med kommandoerne binom.pmf(x,n,p) og binom.cdf(x,n,p) (MATLAB: binopdf(x,n,p) og binocdf(x,n,p)). Til beregning af -værdien skal vi også bruge 1-binom.cdf(x,n,p) som i python kan beregnes med kommandoen binom.sf(x,n,p) (MATLAB: binocdf(x,n,p,'upper')).
I kodevinduet nedenfor beregnes -værdien for test af Kipping og Popes hypotese for data i det foregående afsnit. Koden kan bruges generelt, hvis man selv indskriver det observerede antal (xobs), antal forsøg () og værdien af sandsynlighedsparameteren (p0), der ønskes testet.

MATLAB-kode

Prøv at køre koden som den står. Prøv dernæst at fjerne linjerne med if-delen i udregningen af pval og kør koden igen. Kan du forklare, hvorfor det er nødvendigt at if-delen med i beregningen?

Svar: Binomial p-værdi

I eksemplet er det forventede antal 17, og -værdien er sandsynligheden for værdier og . Hvis imidlertid vi skriver binom.sf(21,34,0.5), får vi sandsynligheden og ikke Derfor må vi skrive et tal lidt under 21, som gøres i programmet ved at trække 0.5 fra. Dog skal vi kun trække noget fra, hvis den øvre grænse xupper er et heltal, hvorfor der i koden er brugt et "if"-statement.

ForegåendeNæste