Soit et deux iidrv où . Je voudrais connaître la distribution pour .
Le mieux que je puisse faire est de prendre la série Taylor des deux et de faire en sorte que la différence soit la somme de la différence entre deux VR normaux et deux RV chi carré en plus du reste de la différence entre les autres termes. Existe-t-il un moyen plus simple d'obtenir la distribution de la différence entre 2 VR log-normaux iid?
Réponses:
C'est un problème difficile. J'ai d'abord pensé à utiliser (une approximation de) la fonction de génération de moment de la distribution lognormale. Cela ne fonctionne pas, comme je vais l'expliquer. Mais d'abord une notation:
Soit la densité normale standard et Φ la fonction de distribution cumulative correspondante. Nous n'analyserons que la distribution log-normale des cas l n N ( 0 , 1 ) , qui a une fonction de densité f ( x ) = 1ϕ Φ lnN(0,1)
et fonction de distribution cumulative
F(x)=Φ(lnx)
Supposons queXetYsont des variables aléatoires indépendantes avec la distribution log-normale ci-dessus. Nous nous intéressons à la distribution deD=X-Y, qui est une distribution symétrique à zéro moyen. SoitM(t)=EetXest la fonction de génération de moment deX. Il est défini uniquement pourt
Cela signifie que nous aurons besoin d' une approche plus directe pour trouver des approximations pour la distribution de . Supposons que t ≥ 0 , calculez P ( D ≤ t )D t≥0
Cette expression peut être utilisée pour l'intégration numérique ou comme base de simulation. D'abord un test:
ce qui est clairement correct. Enveloppons cela dans une fonction:
qui donne:
On peut alors trouver la fonction de densité en se différenciant sous le signe intégral, en obtenant
que nous pouvons tester:
Et en traçant la densité que nous obtenons:
J'ai également essayé d'obtenir une approximation analytique, mais jusqu'à présent, je n'ai pas réussi, ce n'est pas un problème facile. Mais l'intégration numérique comme ci-dessus, programmée en R est très rapide sur du matériel moderne, c'est donc une bonne alternative qui devrait probablement être beaucoup plus utilisée.
la source
Cela ne répond pas strictement à votre question, mais ne serait-il pas plus facile de regarder le rapportX et Oui ? Vous arrivez alors simplement à
Depending on your application, this may serve your needs.
la source