Test binomial exact dans SAS ou SPSS?

Question

J'ai 400 000 lignes, chacune avec un identifiant, le nombre d'essais et le nombre de succès.

J'ai une (valeur) pour p, une probabilité hypothétique de succès.

Je voudrais calculer une p-valeur exacte, binomiale, unilatérale pour chaque rangée, en utilisant ses succès d'essais, et le p global.

Par ex,

devrait me donner

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

Je préférerais une solution dans SAS, mais SPSS 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 les réponses -1, de sorte que vous pourriez avoir besoin de tester nsuccesses=0 avant d'utiliser la formule p_higher.