Find centralized, trusted content and collaborate around the technologies you use most.
Teams
Q&A for work
Connect and share knowledge within a single location that is structured and easy to search.
400,000 行あり、それぞれに ID、試行回数、成功回数が含まれています。
仮定された成功確率である p の (1 つの) 値があります。
試行の成功とグローバルpを使用して、各行の正確な二項片側p値を計算したいと思います。
例えば、
1001 10 2
私に与えるべきです
1001 10 2 2.639011e-01 2 回以上の成功の確率。
SAS でのソリューションが望ましいですが、SPSS も便利です。
ありがとう!
二項分布を直接使用できます。
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
nsuccesses=0
p_higher