近似对数正态和 pdf(在 R 中)

机器算法验证 r 对数正态分布
2022-03-26 07:01:01

我有一个应用程序,我需要一个对数正态和 pdf 的近似值,以用作似然函数的一部分。对数正态和分布没有封闭形式,信号处理期刊上有一堆关于不同近似的论文。我一直在使用最简单的近似值之一(Fenton 1960),它涉及用匹配第一和第二时刻的单个对数正态替换对数正态的总和。这对代码来说非常简单,但从过去 50 年中有关该主题的文献来看,这可能不是所有应用程序的最佳近似值。对于如何确定哪些近似值将导致最佳 MLE 估计,我没有直觉。

有谁知道(A)我应该对最大似然应用程序使用不同的近似值吗?(B) 是否存在任何计算密集型近似值的 R 代码?

更新:有关该问题的一些背景信息,请参阅此评论

1个回答

获得中等分布函数的数值版本N(比如十几个 r.vs 或更少),一种简单的方法是计算每个 LN 密度的离散傅里叶变换 (DFT),形成乘积,然后使用逆 DFT。所有密度都必须使用相同的网格,并且必须小心设计。计算可以很容易地在 R 函数中完成。但是,不要期望达到 R 中经典分布函数的显着精度。