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

Question

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

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

我想使用其试验成功率和全局 p 来计算每一行的精确、二项式、单侧 p 值。

例如，

1001 10 2

应该给我

对于 2 次或更多成功的概率，

1001 10 2 2.639011e-01。我更喜欢 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 。