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

Thomas Lumley tlumley at u.washington.edu
Mon Jan 22 20:37:01 CET 2007


On Mon, 22 Jan 2007, Lynette wrote:

> Dear all, especially to Thomas,
>
> I have figured out the problem. For the step function, something wrong with 
> my C codes. I should use the expression ((x>=0.25)&&(x<=0.75)) ? 2:1 instead 
> of ((x>=1/4)&&(x<=3/4)) ? 2:1 ). Have no idea why 0.25 makes difference from 
> 1/4 in C. But now I can go ahead with the correct integration in C. Thank you 
> all. And hope this helps to others.
>

1 and 4 are ints in C, so 1/4 is 0 (integer division).

 	-thomas



More information about the R-help mailing list