Testen Sie, ob zwei Stichproben von Binomialverteilungen mit demselben p übereinstimmen


9

Angenommen, ich habe getan:

  • p 1 k 1n1 unabhängige Versuche mit einer unbekannten Erfolgsrate und beobachteten Erfolgen.p1k1
  • n2 unabhängige Versuche mit unbekannter Erfolgsrate und beobachteten Erfolgen.k 2p2k2

Wenn jetzt aber immer noch unbekannt, ist die Wahrscheinlichkeit , für ein gegebenes (oder umgekehrt zu beobachten, proportional zu \ int_0 ^ 1 B (n_1, p, k_1) B (n_2, p, k_2) \ text {d} p = \ frac {1} {n_1 + n_2 + 1} \ binom {n_1} {k_1} \ binom {n_2} {k_2} \ binom {n_1 + n_2} {k_1 + k_2 } ^ {- 1} Wenn ich also auf p_1 \ neq p_2 testen möchte , muss ich nur nachsehen, in welchem ​​Quantil der entsprechenden Verteilung meine Beobachtungen sind.p ( k 2 ) k 2 k 1p1=p2=:pp(k2)k2k101B(n1,p,k1)B(n2,p,k2)dp=1n1+n2+1(n1k1)(n2k2)(n1+n2k1+k2)1p1p2

Bisher zur Neuerfindung des Rades. Mein Problem ist nun, dass ich dies in der Literatur nicht finde und daher wissen möchte: Was ist der Fachbegriff für diesen Test oder ähnliches?


2
Warum nicht den Zwei-Proportionen-Z-Test verwenden ( en.wikipedia.org/wiki/Statistical_hypothesis_testing ) (Wenn ich Ihr Problem richtig verstehe).
Verena Haunschmid

@ExpectoPatronum: Auf den ersten Blick besteht das größte Problem darin, dass dieser Test mindestens 5 Erfolge und Misserfolge für jede Beobachtung erfordert, die in meiner Anwendung möglicherweise nicht angegeben sind und auch darauf hinweisen, dass (unnötige) Annäherungen vorgenommen werden.
Wrzlprmft

ok, das ist ein Problem, aber die meisten Tests haben ähnliche Anforderungen.
Verena Haunschmid

@ExpectoPatronum: Auf der Suche nach einer exakten Alternative zum Zwei-Proportionen-Z-Test habe ich den exakten Fisher-Test gefunden, der auf den ersten Blick sehr ähnlich aussieht (aber ich muss ihn noch im Detail untersuchen).
Wrzlprmft

1
@ExpectoPatronum: Die Division spielt keine Rolle, da der große Term nur proportional zu und genau die Normalisierungskonstante ist. Wie auch immer, ich habe jetzt bestätigt, dass dies Fisher's Exact Test ist, den ich dank Ihnen gefunden habe. p(k2)(n1+n2+1)
Wrzlprmft

Antworten:


6

Die Teststatistik ist die des Fisher's Exact Test .p(k2)

Da Normalisierung kann durch Multiplikation mit und somit:

k2n21n1+n2+1(n1k1)(n2k2)(n1+n2k1+k2)1=1n1+n2+1,
n1+n2+1
p(k2)=(n1k1)(n2k2)(n1+n2k1+k2)1.
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.