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