اختبار ذي الحدين الدقيق في SAS أو SPSS؟

Question

لديّ 400000 صف ، ولكل منها هوية ، وعدد من التجارب ، وعدد النجاحات.

لدي قيمة (واحدة) لـ P ، احتمال النجاح المفترض.

أرغب في حساب قيمة P ذات الحدين ، من جانب واحد لكل صف ، باستخدام نجاحات التجارب ، و P.

على سبيل المثال ،

1001 10 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 ، لذلك قد تحتاج إلى اختبار nsuccesses=0 قبل استخدام p_higher معادلة.