Approximating lognormal sum pdf (in R)

I have an application for which I need an approximation to the lognormal sum pdf for use as part of a likelihood function. The lognormal sum distribution has no closed form, and there are a bunch of papers in signal processing journals about different approximations. I have been using one of the simplest approximations (Fenton 1960), which involves replacing a sum of lognormals with a single lognormal with matching first and second moments. This is pretty straightforward to code, but judging by the literature on the subject that has been written in the last 50 years, this may not be the best approximation for all applications. I have no intuition for how to identify which approximations will lead to the best MLE estimates.

Does anyone know if
(A) There is different approximation I should be using for a maximum likelihood application?
(B) There is existing R code for any of the more computationally intensive approximations?

Update: For some background on the problem, see this review

Answer

To obtain a numerical version of distribution function for moderate N (say a dozen of r.vs or less), a simple approach is to compute the Discrete Fourier Transform (DFT) of each LN density, form the product and then use inverse DFT. The same grid must be used for all densities and it must be designed with some care. The computation can be done quite easily in a R function. However, do not expect to reach the remarkable precision of the classical distributions functions in R.

Attribution
Source : Link , Question Author : Ben Lauderdale , Answer Author : Yves

Leave a Comment