1

我需要在标准正态分布中计算 50 的密度,实际上,结果会返回 0 pdf('norm',50,0,1),然后我尝试显示它的对数,它不适用于log(pdf('norm',50,0,1))。在 R 中,我们可以dnorm(50,log = TRUE)在 MATLAB 中做类似的函数吗?我没有找到它。谢谢

4

1 回答 1

2

如果您计算评估为 50 的 pdf,由于精度有限(具体而言,结果小于realmin),您会得到 0。您应该直接计算对数。根据定义,标准正态 pdf 在 50 处评估的对数是

>> format long %// to show more decimals
>> -50^2/2 - .5*log(2*pi)
ans =
   -1.250918938533205e+003
于 2014-06-14T09:47:20.107 回答