Проверьте, соответствуют ли две выборки биномиальных распределений одному и тому же р

9

Предположим, я сделал:

  • p 1 k 1n1 независимых испытаний с неизвестной частотой успеха и наблюдаемыми успехами .p1k1
  • p 2 k 2n2 независимых испытаний с неизвестной частотой успеха и наблюдаемыми успехами .p2k2

Если теперь но все еще неизвестно, вероятность наблюдать для данного (или наоборот) пропорциональна \ 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} , поэтому, если я хочу проверить p_1 \ neq p_2 , мне нужно только посмотреть, в каком квантиле соответствующего распределения находятся мои наблюдения.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

Пока что изобретать велосипед. Теперь моя проблема в том, что я не могу найти это в литературе, и поэтому я хочу знать: каков технический термин для этого теста или что-то подобное?

Wrzlprmft
источник
2
Почему бы не использовать z-тест с двумя пропорциями ( en.wikipedia.org/wiki/Statistical_hypothesis_testing ) (если я правильно понимаю вашу проблему).
Верена Хауншмид
@ExpectoPatronum: На первый взгляд самая большая проблема состоит в том, что для этого теста требуется не менее 5 успехов и неудач для каждого наблюдения, что может быть не указано в моем приложении, а также указывает на то, что (ненужные) приближения сделаны.
Wrzlprmft
хорошо, это проблема, но большинство тестов имеют схожие требования.
Верена Хауншмид
@ExpectoPatronum: В любом случае, в поисках точной альтернативы двухпропорциональному z-критерию, я нашел точный критерий Фишера, который на первый взгляд очень похож (но я еще не разобрался в нем подробно).
Wrzlprmft
1
@ExpectoPatronum: Деление не имеет значения, так как большой член пропорционален только и - это в точности константа нормализации. Во всяком случае, теперь я подтвердил, что это точный тест Фишера, который я нашел благодаря вам. ( n 1 + n 2 + 1 )p(k2)(n1+n2+1)
Wrzlprmft

Ответы:

6

Статистика испытаний является то , что точного критерия Фишера .p(k2)

Так как нормализация может быть получена умножением на и, таким образом: n1+n2+1p(k2)= ( n1

k2n21n1+n2+1(n1k1)(n2k2)(n1+n2k1+k2)1=1n1+n2+1,
n1+n2+1
p(k2)=(n1k1)(n2k2)(n1+n2k1+k2)1.
Wrzlprmft
источник