[R] Ratio of huge products
Alberto Magni
alberto.magni86 at gmail.com
Thu Dec 15 19:35:11 CET 2011
Hello everybody,
I have to compute something in this form:
x = prod(a:b) / prod(c:d), where: a < c and b < d and obviously: a
< b and c < d
I cannot make assumptions on the relative position of c,b and a,d.
The problem is that a,b,c,d are large and the products are huge (R return Inf).
Their ratio is less than 1 but significantly higher than 0: it is a
non-tiny probability.
I need to find a way to simplify this ratio.
The only way to solve this that I see is to decompose into prime
factors all the
numbers in the numerator and the denominator and to remove the ones in common
Do you know a better way to do this ?
Thank you,
Alberto
More information about the R-help
mailing list