SASまたはSPSSの正確な二項テスト？

Question

私には400,000行があり、それぞれにID、試行数、成功の数があります。

私はPに対して（1つの）価値を持っています。これは、成功の仮定された確率です。

試行の成功とグローバル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 方式。