SAS或SPSS中的准确二项式检验？

Question

我有400,000行，每行都有一个id，试验次数和成功次数。

我有（一）p值，一个假设的成功概率。

我想计算一个确切的，二项式，单边p值为每行，使用它的试验成功，和全球p。

例如，

应该给我

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的反应，所以你可能需要使用p_higher公式，用于测试nsuccesses=0。