Working with numbers smaller tha 10^-308

5 次查看（过去 30 天）

sadoune sadoune 2013-5-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75223-working-with-numbers-smaller-tha-10-308

Hi,

I would like to know if it is possible to use Matlab for very small numbers (smaller than 10^-308 ). It's very important for my application "statistical genetics". in this field, somes probabilities are smaller than 10^-500.

Thank you

4 个评论
显示 2更早的评论隐藏 2更早的评论

sadoune sadoune 2013-5-9

Hi Robert,

You can imagine a vector length 500, each case represent a probability, and the probability i would like to compute is the product of all this probabilities (independance).

John Doe 2013-5-9

Ok, then I see where you get the number from, but if you are ending up with a probability that low, it is literally impossible for it to happen. Isn't OK to just let it be equal to 0?

I don't mean to be rude or condescending, I just don't get how it would make sense to keep working with numbers that low.

Anyhow, I hope someone can provide an answer for you, and good luck =)

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Sean de Wolski 2013-5-9

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75223-working-with-numbers-smaller-tha-10-308#answer_84933

编辑：Sean de Wolski 2013-5-9

在 MATLAB Online 中打开

You can do this using vpa() in the Symbolic Math Toolbox:

v = vpa('3.14159^-5000')

And for more info:

doc vpa

You could also use exact symbolic arithmetic:

sym('3.1415^-34334')

doc sym

3 个评论
显示 1更早的评论隐藏 1更早的评论

Sean de Wolski 2013-5-10

在 MATLAB Online 中打开

That's because you're doing the computation as a double and then converting it to vpa(). Convert it to vpa before doing the computation:

for i=400:-1:1
  tab(i)=vpa(0.1);
end
v = (prod(tab))

sadoune sadoune 2013-5-10

Thank you Wolski. it's OK now

请先登录，再进行评论。

Peter Perkins 2013-5-11

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75223-working-with-numbers-smaller-tha-10-308#answer_85077

Quite often in statistics when one has very small probabilities, one works on the log scale. That's why one typically maximizes the LOG-likelihood, and not the likelihood.