问题描述
我将计算矩阵法向密度以在对数似然公式中使用它们。
在我的计算中，我需要计算大数（数千万）的指数。我意识到，当 i >=710 时，R 为 exp(i) 返回无穷大。
无论如何我可以手动强制 R 不显示无穷大，或者它的日志相应地不为 0？
感谢您纠正我的错字
解决方法
您应该知道 R 中的最大双精度数是 .Machine$double.xmax。如果你有比这更好的数字，它肯定会给你Inf
> log(.Machine$double.xmax)
[1] 709.7827
还有一个贡献包 Brobdingnag 可以处理比 exp(710) 大得多的数字。
library(Brobdingnag)
exp(as.brob(710))
#[1] +exp(710)
exp(as.brob(10^seq(100,310,by = 5)))
# [1] +exp(1e+100) +exp(1e+105) +exp(1e+110) +exp(1e+115) +exp(1e+120)
# [6] +exp(1e+125) +exp(1e+130) +exp(1e+135) +exp(1e+140) +exp(1e+145)
#[11] +exp(1e+150) +exp(1e+155) +exp(1e+160) +exp(1e+165) +exp(1e+170)
#[16] +exp(1e+175) +exp(1e+180) +exp(1e+185) +exp(1e+190) +exp(1e+195)
#[21] +exp(1e+200) +exp(1e+205) +exp(1e+210) +exp(1e+215) +exp(1e+220)
#[26] +exp(1e+225) +exp(1e+230) +exp(1e+235) +exp(1e+240) +exp(1e+245)
#[31] +exp(1e+250) +exp(1e+255) +exp(1e+260) +exp(1e+265) +exp(1e+270)
#[36] +exp(1e+275) +exp(1e+280) +exp(1e+285) +exp(1e+290) +exp(1e+295)
#[41] +exp(1e+300) +exp(1e+305) +exp(Inf)
而且 10^305 比 710 大得多。请参阅包 vignette。