[R] Round error?

[Bob Wheeler]

You might try pghyper() in the SuppDists package. It should give the 
same values for integer arguments.


Peter Dalgaard wrote:
> Aurora Torrente <aurora at ebi.ac.uk> writes:
> 
> 
>>Hi all,
>>
>>I have tried to compute a p-value for a hypergeometric distribution as:
>>
>>dhyper(x,k,l,n) + phyper(x,k,l,n,lower.tail=FALSE)
>>
>>and sometimes obtained negative values.  Do you know if it is because
>>a round error or am I doing something wrong?
> 
> 
> Could you give us specific x,k,l,n where this happens, please. It
> sounds like it shouldn't happen.
> 


-- 
Bob Wheeler --- http://www.bobwheeler.com/
         ECHIP, Inc. ---
Randomness comes in bunches.