Re: [Rd] minor flaw in integrate()

2007-07-03 Thread Gabor Grothendieck
If integrate is changed it would be nice at the same time to make it into an S3 generic. deriv already is an S3 generic but strangely integrate is not. Ryacas provides a deriv method but for integrate Ryacas inconsistently provides Integrate since integrate is not generic. On 6/28/07, Peter Ruck

Re: [Rd] minor flaw in integrate()

2007-07-03 Thread Prof Brian Ripley
I think throwing an error is a better solution: this is rather unlikely to be deliberate and returning NaN might postpone the detection too long. On Tue, 3 Jul 2007, Duncan Murdoch wrote: > On 7/3/2007 11:55 AM, Martin Maechler wrote: >>> "PetRd" == Peter Ruckdeschel <[EMAIL PROTECTED]> >>>

Re: [Rd] minor flaw in integrate()

2007-07-03 Thread Duncan Murdoch
On 7/3/2007 11:55 AM, Martin Maechler wrote: >> "PetRd" == Peter Ruckdeschel <[EMAIL PROTECTED]> >> on Tue, 03 Jul 2007 17:26:43 +0200 writes: > > PetRd> Thanks Martin and Duncan for your > PetRd> comments, > > PetRd> Martin Maechler wrote: > >>> "DM" == Duncan Mur

Re: [Rd] minor flaw in integrate()

2007-07-03 Thread Martin Maechler
> "PetRd" == Peter Ruckdeschel <[EMAIL PROTECTED]> > on Tue, 03 Jul 2007 17:26:43 +0200 writes: PetRd> Thanks Martin and Duncan for your PetRd> comments, PetRd> Martin Maechler wrote: >>> "DM" == Duncan Murdoch <[EMAIL PROTECTED]> >>> on Mon, 02 Jul 2007 21

Re: [Rd] minor flaw in integrate()

2007-07-03 Thread Peter Ruckdeschel
Thanks Martin and Duncan for your comments, Martin Maechler wrote: >> "DM" == Duncan Murdoch <[EMAIL PROTECTED]> >> on Mon, 02 Jul 2007 21:56:23 -0400 writes: > > DM> On 28/06/2007 5:05 PM, Peter Ruckdeschel wrote: > >> Hi, > >> > >> I noticed a minor flaw in integrat

Re: [Rd] minor flaw in integrate()

2007-07-03 Thread Martin Maechler
> "DM" == Duncan Murdoch <[EMAIL PROTECTED]> > on Mon, 02 Jul 2007 21:56:23 -0400 writes: DM> On 28/06/2007 5:05 PM, Peter Ruckdeschel wrote: >> Hi, >> >> I noticed a minor flaw in integrate() from package stats: >> >> Taking up arguments lower and upper from

Re: [Rd] minor flaw in integrate()

2007-07-02 Thread Duncan Murdoch
On 28/06/2007 5:05 PM, Peter Ruckdeschel wrote: > Hi, > > I noticed a minor flaw in integrate() from package stats: > > Taking up arguments lower and upper from integrate(), > >if (lower == Inf) && (upper == Inf) > >or > >if (lower == -Inf) && (upper == -Inf) > > integrate()

[Rd] minor flaw in integrate()

2007-06-28 Thread Peter Ruckdeschel
Hi, I noticed a minor flaw in integrate() from package stats: Taking up arguments lower and upper from integrate(), if (lower == Inf) && (upper == Inf) or if (lower == -Inf) && (upper == -Inf) integrate() calculates the value for (lower==-Inf) && (upper==Inf). Rather, it shoul