I'd suggest that the original sin here is calling some particular
numerical integration routine 'integrate', which gives the user an
illusory sense of power Functions have to be well-behaved in
various ways for quadrature to work well, and you've got to expect
things like
Hi,
I was recently given some interesting tips on a similar issue, see
R-help puzzle with integrate over infinite range
http://www.r-help.com/list/85/713882.html
Maybe fails can be a bit misleading here (fails to produce the
actual result vs. returning an error message). As a result of this
Prof Brian Ripley rip...@stats.ox.ac.uk
on Tue, 7 Dec 2010 07:41:16 + (GMT) writes:
On Mon, 6 Dec 2010, Spencer Graves wrote:
Hello:
The example integrate(dnorm,0,2) says it fails on many systems.
I just got 0 from it, when I should have gotten
Ripley rip...@stats.ox.ac.uk
From: Martin Maechler
Sent by: r-devel-boun...@r-project.org
Date: 12/07/2010 03:29AM
Subject: Re: [Rd] 0.5 != integrate(dnorm,0,2) = 0
Prof Brian Ripley rip...@stats.ox.ac.uk
on Tue, 7 Dec 2010 07:41:16 + (GMT) writes:
On Mon, 6 Dec 2010, Spencer
-project.org
Date: 12/07/2010 03:29AM
Subject: Re: [Rd] 0.5 != integrate(dnorm,0,2) = 0
Prof Brian Ripleyrip...@stats.ox.ac.uk
on Tue, 7 Dec 2010 07:41:16 + (GMT) writes:
On Mon, 6 Dec 2010, Spencer Graves wrote:
Hello:
The example integrate(dnorm
and
pick an appropriate region of integration.
John Nolan, American U.
-r-devel-boun...@r-project.org wrote: -
To: r-devel@r-project.org
From: Pierre Chausse
Sent by: r-devel-boun...@r-project.org
Date: 12/07/2010 09:46AM
Subject: Re: [Rd] 0.5 != integrate(dnorm,0,2) = 0
The warning
On Mon, 6 Dec 2010, Spencer Graves wrote:
Hello:
The example integrate(dnorm,0,2) says it fails on many systems.
I just got 0 from it, when I should have gotten either an error or something
close to 0.5. I got this with R 2.12.0 under both Windows Vista_x64 and
Linux (Fedora 13);