Re: [R] How to get correct integration in C for step function?

From: Thomas Lumley <tlumley_at_u.washington.edu>
Date: Mon 22 Jan 2007 - 15:48:14 GMT

On Sun, 21 Jan 2007, Lynette wrote:

> Dear all,
>
> I am using Rdqags in C to realize the integration. It seems for the
> continous C function I can get correct results. However, for step functions,
> the results are not correct. For example, the following one, when integrated
> from 0 to 1 gives 1 instead of the correct 1.5
>

Using integrate() in R for an R-defined step function gives the right answer (eg on the example in ?ecdf).

This suggests a problem in your C code, since integrate() just calls dqags.

         -thomas



R-help@stat.math.ethz.ch mailing list
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. Received on Tue Jan 23 02:51:29 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Mon 22 Jan 2007 - 17:30:32 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.