我有一个应用程序,我需要一个对数正态和 pdf 的近似值,以用作似然函数的一部分。对数正态和分布没有封闭形式,信号处理期刊上有一堆关于不同近似的论文。我一直在使用最简单的近似值之一(Fenton 1960),它涉及用匹配第一和第二时刻的单个对数正态替换对数正态的总和。这对代码来说非常简单,但从过去 50 年中有关该主题的文献来看,这可能不是所有应用程序的最佳近似值。对于如何确定哪些近似值将导致最佳 MLE 估计,我没有直觉。
有谁知道(A)我应该对最大似然应用程序使用不同的近似值吗?(B) 是否存在任何计算密集型近似值的 R 代码?
更新:有关该问题的一些背景信息,请参阅此评论