[R] Make check fails on d-p-q-r-tests.R...
Prof Brian Ripley
ripley at stats.ox.ac.uk
Tue Sep 27 23:38:24 CEST 2005
I am afraid this does look like a real problem, if a minor one. We have
for the first problem
x <- 10^(ex <- c(1,2,5*(1:5),50,100,200,300,Inf))
ex <- -c(rev(1/x), ex)
qcauchy(ex, log=TRUE)
The first entry of the result should be Inf and the last -Inf. From the
output you have shown us I would guess the last is NaN.
If this is what is going on, it is a problem but a somewhat esoteric one
that probably reflects a lack of IEC60559 (aka IEEE754) conformance. It
needs more hands-on debugging than we can provide on a help list.
On Tue, 27 Sep 2005, Jeff Ross wrote:
> On 1:03:39 am 09/27/05 Prof Brian Ripley <ripley at stats.ox.ac.uk> wrote:
>
>> You will have to show us the error! It will be shown in
>> d-p-q-r-tests.Rout.fail (unless this was a segfault or similar).
>>
>> It is not OK to skip the test, but note that this test is random and
>> does fail about 1 in 50 times, so you could just try rerunning it.
>
> I am so sorry. I didn't even think to look for a .fail file.
>
> Here is the last bit of it:
>
>>> ## for PR#7902:
>> ex <- -c(rev(1/x), ex)
>> All.eq(-x, qcauchy(pcauchy(-x)))
> [1] TRUE
>> All.eq(+x, qcauchy(pcauchy(+x, log=TRUE), log=TRUE))
> [1] TRUE
>> All.eq(1/x, pcauchy(qcauchy(1/x)))
> [1] TRUE
>> All.eq(ex, pcauchy(qcauchy(ex, log=TRUE), log=TRUE))
> [1] "`is.NA' value mismatches: 1 in current, 0 in target"
> Warning message:
> NaNs produced in: qcauchy(p, location, scale, lower.tail, log.p)
>> II <- c(-Inf,Inf)
>> stopifnot(pcauchy(II) == 0:1, ## qcauchy(0:1) == II,
> + pcauchy(II, log=TRUE) == c(-Inf,0),
> + qcauchy(c(-Inf,0), log=TRUE) == II)
> Error in if (!(is.logical(r <- eval(ll[[i]])) && all(r)))
> stop(paste(deparse(mc[[i + :
> missing value where TRUE/FALSE needed
> In addition: Warning message:
> NaNs produced in: qcauchy(p, location, scale, lower.tail, log.p)
> Execution halted
>
> This is from the latest R-patched source tar ball, and it is the identical
> error as R-2.1.1. I've consistently gotten the same error, even with
> configuration file tweaks and making clean between runs.
>
> Thanks!
>
> Jeff Ross
--
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
More information about the R-help
mailing list