> On Jan 26, 2017, at 12:51 PM, George Trojan - NOAA Federal > <george.tro...@noaa.gov> wrote: > > Why this happens? > >> library("R2Cuba") >> f <- function(x) 0 >> cuhre(ndim=1, ncomp=1, f) > Iteration 1: 11 integrand evaluations so far > [1] 0 +- 0 chisq 0 (0 df) > Iteration 2: 33 integrand evaluations so far > [1] 0 +- 0 chisq 0 (1 df) > integral: 0 (+-0) > nregions: 2; number of evaluations: 33; probability: 0 >> vegas(ndim=1, ncomp=1, f) > Iteration 1: 1000 integrand evaluations so far > [1] 0 +- 7.02518e-18 chisq 0 (0 df) > . . . > Iteration 13: 52000 integrand evaluations so far > [1] 0 +- 2.65413e-18 chisq 0 (12 df) > integral: 0 (+-2.7e-18) > number of evaluations: 52000; probability: 0 > failed with message ‘Dimension out of range’ >> suave(ndim=1, ncomp=1, f) > . . . > [1] 0 +- 7.02462e-17 chisq 0 (411 df) > integral: 0 (+-7e-17) > nregions: 50; number of evaluations: 50000; probability: 0 > failed with message ‘Dimension out of range’ > > Is this a bug?
If you look at the code you will see that the delivered message does not actually match what the documentation says should be delivered when 'ifail' is (positive) 1. Appears the package maintainer has not put in code to warn of failure to reach "the accuracy goal was not met within the allowed maximum number of integrand evaluations." That message is only supposed to be delivered when 'ifail' is -1. > maintainer("R2Cuba") [1] "Annie Bouvier <annie.bouv...@jouy.inra.fr>" > > George > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. David Winsemius Alameda, CA, USA ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.