我正在将我的数据集 $x \in [60,80]$ 的归一化直方图拟合到 Nakagami 分布。首先,我通过以下 MLE 代码估计了 packagednaka
的比例和形状参数:VGAM
ll <- function(par) {
if(par[1]>0 & par[2]>0) {return(-sum(log(dnaka(x, scale = par[1], shape = par[2]) ) ) )} # m=shape, ohm or spread = scale
else return(Inf)
}
mle = optim(c(1000,1), ll)
然后,我通过以下代码根据估计的参数估计对数似然值:
lik = sum(log(dnaka(x, shape = mle$par[1], scale = mle$par[2]) ) )
但是对数似然值 lik 是-Inf
。我知道这个无限值是由于 Nakagami 分布的 PDF 方程中的 exp(.) 项。有没有办法估计我的数据集 $x \in [60,80]$ 的 Nakagami 分布的有限对数似然值?谢谢你。