Точный биномиальный тест в SAS или SPSS?

Question

У меня 400 000 строк, каждая с идентификатором, количеством испытаний и количеством успехов.

У меня есть (одно) значение для p, предполагаемая вероятность успеха.

Я бы хотел вычислить точное, биномиальное одностороннее значение p для каждой строки, используя результаты испытаний и глобальное значение p.

Например,

1001 10 2

должен дать мне

1001 10 2 2.639011e-01 для вероятности 2 или более успехов.

Я бы предпочел решение в SAS, но SPSS также полезен.

Спасибо!

Answer 1

Вы можете использовать биномиальное распределение напрямую:

data yourdata;
   set yourdata;
   p_lower = PROBBNML(globalp, ntrials, nsuccesses);
   p_higher = 1 - PROBBNML(globalp, ntrials, nsuccesses-1);
run;

Я не проверял, будет ли PROBBNML работать с -1 ответами, поэтому вам может потребоваться проверить на nsuccesses=0 перед использованием формулы p_higher.