Joseph Coveney writes:
> But is there any way to improve accuracy in the following direction?
> Or is it that I'm not understanding these functions correctly?
>
> . display invbinomial(100, 50, Binomial(100, 50, 0.5))
> .50996594
>
> . display invbinomial(100, 50, Binomial(100, 50, 0.999))
> .85115891
Please note the definition for invbinomial() -- it is a bit unusual:
from -help invbinomial-:
invbinomial(n,k,P) returns the inverse binomial: for P<=0.5, probability p
such that the probability of observing k or more successes in n trials is P;
for P>0.5, probability p such that the probability of observing k or fewer
successes in n trials is 1-P.
Hope that helps.
--Jean Marie
jlinhart@stata.com
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/