如何防止R显示exp和日志值inf和zero

Question

我将计算矩阵正态密度，以便在对数似然公式中使用它们
在我的计算中，我需要计算大数(数千万(的指数。我意识到，当I>710.
我是否可以手动强制R不显示无穷大，或者它的log相应地不为0？

感谢您更正我的拼写错误

Answer 1

您应该知道R中的最大二重数是.Machine$double.xmax。如果你有比这更好的数字，它肯定会给你Inf

> log(.Machine$double.xmax)  
[1] 709.7827

Answer 2

还有一个贡献包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大得多。请参阅包装小插曲。