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

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

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

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

E.g.,

1001 10 2

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

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

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

Спасибо!


sas spss
person Gerry    schedule 21.06.2010    source источник

Ответы (1)


arrow_upward
3
arrow_downward

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

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.

person Aniko    schedule 23.06.2010