Exact Binomialtest in SAS oder SPSS?

Question

Ich habe 400.000 Zeilen, die jeweils mit einer ID, die Anzahl der Versuche und die Anzahl der Erfolge.

Ich habe (eine) Wert für p, a hypothetisch angenommenen Erfolgswahrscheinlichkeit.

Ich mag einen genau berechnen, binomischen, einseitigen p-Wert für jede Zeile, seine Studien Erfolge mit, und die globalen p.

z. B.

1001 10 2

sollte mir geben

1001 10 2 2.639011e-01 für die Wahrscheinlichkeit von 2 Erfolgen oder mehr.

Ich würde eine Lösung in SAS bevorzugen, aber SPSS ist auch nützlich.

Danke!

Answer 1

Sie können die Binomialverteilung direkt verwenden:

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

Ich habe nicht geprüft, ob PROBBNML funktioniert mit -1 Antworten, so dass Sie für nsuccesses=0 zu Test benötigen könnten, bevor die p_higher Formel.