Ungefähre logarithmische Summe pdf (in R)

8

Ich habe eine Anwendung, für die ich eine Annäherung an die lognormale Summe pdf zur Verwendung als Teil einer Wahrscheinlichkeitsfunktion benötige. Die logarithmische Normalverteilung hat keine geschlossene Form, und es gibt eine Reihe von Artikeln in Signalverarbeitungsjournalen über verschiedene Näherungen. Ich habe eine der einfachsten Näherungen verwendet (Fenton 1960), bei der eine Summe von Lognormalen durch ein einzelnes Lognormal ersetzt wird, wobei der erste und der zweite Moment übereinstimmen. Dies ist ziemlich einfach zu codieren, aber nach der Literatur zu dem Thema, die in den letzten 50 Jahren geschrieben wurde, ist dies möglicherweise nicht die beste Annäherung für alle Anwendungen. Ich habe keine Ahnung, wie ich feststellen kann, welche Annäherungen zu den besten MLE-Schätzungen führen.

Weiß jemand, ob (A) es eine andere Annäherung gibt, die ich für eine Anwendung mit maximaler Wahrscheinlichkeit verwenden sollte? (B) Gibt es einen R-Code für eine der rechenintensiveren Näherungen?

Update: Hintergrundinformationen zum Problem finden Sie in dieser Rezension

Ben Lauderdale
quelle
Können Sie nur eine Berührung klarstellen? Ist das, was Sie als "lognormale Summe pdf" bezeichnen, die Dichtefunktion von wobei mit den Parametern und lognormal ist ? X n μ σ 2Y=X1++XnXnμσ2
Kardinal
Ja, das PDF für die Summe der logarithmischen Normalvariablen variiert.
Ben Lauderdale
Wie groß ist in Ihrer Anwendung? n
Kardinal
Ich interessiere mich am meisten für Fälle, in denen N klein ist, <10 oder so. Es wäre jedoch sehr hilfreich, wenn ich mindestens N bis zu 100 oder so verwalten könnte.
Ben Lauderdale
1
Ein Moment, der einem Lognormal entspricht, klingt an der Oberfläche wie eine seltsame Idee. Dies liegt daran, dass das Lognormal nicht durch seine Momente gekennzeichnet ist. Ich werde hier nachsehen, aber vielleicht gibt es eine Möglichkeit, das Problem ein wenig umzukehren. Sei eine "Standard" , ). Definieren für . Dann ist ein PDF und und haben für jedes solche die gleichen Momente . μ = 0 σ = 1 b ( - 1 , 1 ) f b ( x ) = f 0 ( x ) ( 1 + b sin ( 2 π log x ) ) f b f 0 f b bf0(x)μ=0σ=1b(1,1)fb(x)=f0(x)(1+bsin(2πlogx))fbf0fbb
Kardinal

Antworten:

2

Um eine numerische Version der Verteilungsfunktion für moderates (z. B. ein Dutzend rvs oder weniger) zu erhalten, besteht ein einfacher Ansatz darin, die diskrete Fourier-Transformation (DFT) jeder LN-Dichte zu berechnen, das Produkt zu bilden und dann die inverse DFT zu verwenden. Für alle Dichten muss das gleiche Gitter verwendet werden, und es muss mit einiger Sorgfalt entworfen werden. Die Berechnung kann ganz einfach in einer R-Funktion durchgeführt werden. Erwarten Sie jedoch nicht, die bemerkenswerte Präzision der klassischen Verteilungsfunktionen in R zu erreichen.N

Yves
quelle