2, I need to calculate e.g. exp(a)/(exp(b)+c), where both a and b
are very
large numbers (>>1000, e.g a=1000, b=1007, and c=5). R gives me NaN
when I
use the following command:
exp(1000)/(exp(1007)+5)
[1] NaN
I am pretty sure this should be close to zero.
No. It should be close to 1. Try 1000/(1000+5) for a simpler example.
Shouldn't it rather be 1 / exp(7), since exp(1007) = exp(1000) * exp(7)?
You can get the correct solution if you rewrite the fraction
algebraically first, dividing all terms by exp(1000).
1 / (exp(7) + 5 * exp(-1000))
[1] 0.000911882
Of course, you have to know which of the numbers are large in order to
do this right ...
Best regards,
Stefan Evert
[ stefan.ev...@uos.de | http://purl.org/stefan.evert ]
______________________________________________
R-help@r-project.org 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.
