Exact test binomial SAS ou SPSS?

Question

Je 400.000 lignes, chacun avec un identifiant, nombre d'essais, et le nombre de succès.

Je (un) la valeur de p, la probabilité hypothétique de succès.

Je voudrais calculer exactement, binomiale, à sens unique valeur p pour chaque ligne, en utilisant ses succès des essais, et la p mondiale.

par exemple.

1001 10 2

devrait me donner

1001 10 2 2.639011e-01 pour la probabilité de succès de 2 ou plus.

Je préfère une solution SAS, SPSS, mais est également utile.

Merci!

Answer 1

Vous pouvez utiliser la distribution binomiale directement:

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

Je ne l'ai pas testé si PROBBNML travaillera avec -1 réponses, de sorte que vous pourriez avoir besoin de test pour nsuccesses=0 avant d'utiliser la formule p_higher.