From: Well Howell <whowell@superlink.net>
To: gsl-discuss@sourceware.org
Subject: Re: accuracy of gsl_cdf_binomial_P
Date: Sun, 05 Jun 2011 13:04:00 -0000 [thread overview]
Message-ID: <4DEB7EA2.7090500@superlink.net> (raw)
In-Reply-To: <459650.92078.qm@web110504.mail.gq1.yahoo.com>
An interesting (but "homework-like" ~;) question - and fun to answer too.
Anyway, I'd probably compare GSL results with those from other sources.
I had easy access to gsl_cdf_binomial_P (v 1.14), R pbinom(k,n,p),
binomCDF
(Excel 2007) and dcdflib (Fortran - Brown, Lovato & Russel; U. Texas;
November, 1997).
For a sample size of n=1000, a trial probability of p=0.01 and number of
successes of
s=1 thru 40, the CDF values from dcdclib and the R 2.13.0 stats package
pbinom()
function (http://cran.r-project.org/) show no difference.
Mean absolute deviations for these 40 tests, comparing pbinom with
gsl_cdf_binomial_P and with binomCDF, show MAD of 2.319E-15 and 3.296E-15
respectively.
My "commend"? Looks as if we all have to decide when to STOP
accumulating small terms, and some stop earlier than others. While I always
test functions in Excel against other sources before release in a report,
anything showing a MAD below 4E-15 sure beats using my slide rule
(which didn't have an incomplete beta function anyway ~;).
Well Howell
On 6/2/2011 12:49 AM, Z F wrote:
> Hello everybody,
>
> I was wondering if someone could comment on the accuracy of gsl_cdf_binomial_P() function gsl implementation for large n (n is about a few thousand).
> for different values of p and when the result of cdf is in the tails ( small less then 0.05 and large -- above 0.95)
>
> Thank you very much
>
> ZF
>
>
next prev parent reply other threads:[~2011-06-05 13:04 UTC|newest]
Thread overview: 4+ messages / expand[flat|nested] mbox.gz Atom feed top
2011-06-02 4:50 Z F
2011-06-05 13:04 ` Well Howell [this message]
2011-06-06 2:20 ` Z F
2011-06-06 10:48 ` Well Howell
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=4DEB7EA2.7090500@superlink.net \
--to=whowell@superlink.net \
--cc=gsl-discuss@sourceware.org \
--cc=well@wheatstone-analytics.com \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).