Valore atteso e varianza del registro (a)


20

Ho una variabile casuale X(a)=log(a) dove a è distribuito normalmente ( μ , σ 2 )N(μ,σ2) . Cosa posso dire di E(X) e Var(X) ? Anche un'approssimazione sarebbe utile.


4
Penso che la domanda riguardasse l '"inverso" del log-normale, cioè dove un normale rv A porta al log-normale X = exp (A), l'interrogatore stava chiedendo della distribuzione di X = log (A), che non è definito (a causa della necessità talvolta di registrare un numero negativo). Potrebbero esserci dei risultati per un normale troncato, ma è probabile che siano disordinati.
Martin O'Leary,

2
rocksportrocker, come sottolinea @Martin O'Leary, non è matematicamente possibile avere una tale variabile X , perché log(a) non è definito per valori negativi. Come minimo è necessario troncare a a valore non negativo. Ci può dire perché si crede a potrebbe essere normale?
whuber

Risposte:


23

Se consideriamo "approssimazione" in senso abbastanza generale possiamo arrivare da qualche parte.

Non dobbiamo presumere che abbiamo una distribuzione normale effettiva, ma qualcosa che è approssimativamente normale, tranne che la densità non può essere diversa da zero in un quartiere di 0.

Quindi diciamo che è "approssimativamente normale" (e concentrati vicino al * media) in un senso che possiamo handwave via le preoccupazioni per un avvicinarsi 0 (e il suo conseguente impatto sui momenti di log ( una ) , perché un doesn 't' get down near 0 '), ma con gli stessi momenti di ordine basso della distribuzione normale specificata, potremmo usare le serie di Taylor per approssimare i momenti della variabile casuale trasformata .aalog(a)a

Per alcune trasformazioni , ciò comporta l'espansione di g ( μ X + X - μ X ) come una serie di Taylor (pensa g ( x + h ) in cui μ X assume il ruolo di ' x ' e X - μ X assume il ruolo di " h ") e quindi prendere le aspettative e quindi calcolare la varianza o l'attesa del quadrato dell'espansione (da cui è possibile ottenere la varianza).g(X)g(μX+XμX)g(X+h)μXXX-μXh

Le aspettative e la varianza approssimative risultanti sono:

E[g(X)]g(μX)+g"(μX)2σX2

Var[g(X)](g'(μX))2σX2

e così (se non ho commesso errori), quando :g()=log()

E[log(un')]log(μun')-σun'22μun'2

Var[log(un')]σun'2/μun'2

* Perché questo sia una buona approssimazione in genere si desidera che la deviazione standard di essere piuttosto piccola rispetto alla media (basso coefficiente di variazione).un'


2
Poiché la serie Taylor per il registro ha un raggio di convergenza relativamente piccolo, si consiglia cautela nell'applicare queste approssimazioni.
whuber

@whuber per un'espansione attorno alla media, penso che ciò corrisponderebbe al consiglio che la "deviazione standard di dovrebbe essere piuttosto piccola rispetto alla media" con cui termina la mia risposta - se mi manca qualche ulteriore problema che tale consiglio non copre dovrei risolvere la mia risposta. un'
Glen_b -Restate Monica

3
L'approssimazione per la media funziona abbastanza bene per e quella per la varianza funziona abbastanza bene per μ / σ > 2,5 o giù di lì. μ/σ>1.5μ/σ>2.5
whuber

In ogni caso, vale sicuramente la pena di essere chiari sul fatto che stiamo indirettamente facendo affidamento sulla convergenza di (poiché ln ( μ + y - μ ) = ln [ μ { 1 + ( y - μ ) / μ } ] = ln ( μ ) + ln [ 1 + ( y - μ ) / μ ]ln(1+X)ln(μ+y-μ)=ln[μ{1+(y-μ)/μ}]=ln(μ)+ln[1+(y-μ)/μ]). Grazie anche per i valori espliciti suggeriti; semmai forse sono leggermente esagerato quando lo uso. Due commenti preziosi.
Glen_b -Restate Monica
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.