[R] Peculiarity in non-central qchisq for ncp > 294.92 ...

Andrew Robinson andrewr at uidaho.edu
Mon Apr 14 13:49:32 CEST 2003


I'm sympathetic to the package author.  I'm running some pretty extreme
simulations, and am quite satisfied by artificially constraining the ncp
to be below 294.  This makes my conclusions slightly conservative, but
only very slightly, and the results are very unlikely to change much when
the proper value is used.

I respectfully submit that undocumented shortcomings are bugs.  This one
did me no damage because I was too impatient to wait for the result.
However, had my simulations been much longer term, I would have obtained
very inaccurate results!

I would be more than pleased to draft a note for the help file for this
and any other that you know of.  I'd view that as being higher priority
than fixing it (them).

Best wishes


On Mon, 14 Apr 2003, Peter Dalgaard BSA wrote:

> Martin Maechler <maechler at stat.math.ethz.ch> writes:
> > which is not really good either:
> >
> >       plot(function(x) pchisq(x, df=1, ncp=295), from=1,to=1e4, log='x')
> >
> > shows the pretty bad behavior (a clear bug).
> Well... shortcoming, I'd say. The author of those routines didn't
> expect people to use them with huge noncentralities. I've had this on
> my TODO list for a couple of years by now, ever since fixing the
> similar issue with dchisq--see Rnews 1/1. It looks fairly easy to do
> something along the same lines, but there is now a double recurrence
> relation and it requires a bit of concentration to get the tail bounds
> for the termination criterion.
> --
>    O__  ---- Peter Dalgaard             Blegdamsvej 3
>   c/ /'_ --- Dept. of Biostatistics     2200 Cph. N
>  (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
> ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907

Andrew Robinson			     Ph: 208 885 7115
Department of Forest Resources	     Fa: 208 885 6226
University of Idaho		     E : andrewr at uidaho.edu
PO Box 441133			     W : http://www.uidaho.edu/~andrewr
Moscow ID 83843			     Or: http://www.biometrics.uidaho.edu
No statement above necessarily represents my employer's opinion.

More information about the R-help mailing list