Displaying 1 result from an estimated 1 matches for "r_q_p01_check".
2005 May 27
1
qcauchy accuracy (PR#7902)
...cauchy has been fixed, it is becoming clear that qcauchy suffers from
the same problems.
qcauchy(pcauchy(1e100,0,1,FALSE,TRUE),0,1,FALSE,TRUE)
should yield 1e100 back, but I get 1.633178e+16. The code below does much
better. Notes:
1. p need not be finite. -Inf is ok in the log_p case and R_Q_P01_check
already checks things.
2. No need to disallow scale=0 and infinite location.
3. The code below uses isnan and finite directly. It needs to be adapted to
the
R way of doing that.
double
qcauchy (double p, double location, double scale, int lower_tail, int log_p)
{
if (isnan(p) || isnan(l...