Distribuzione della convoluzione delle variabili quadrate normali e chi-quadrate?


14

il seguente problema è emerso di recente durante l'analisi dei dati. Se la variabile casuale X segue una distribuzione normale e Y segue una distribuzione (con n dof), come viene distribuito Z = X 2 + Y 2 ? Fino ad ora ho inventato il pdf di Y 2 : ψ 2 n ( x )χn2Z=X2+Y2Y2

ψn2(X)=F(X)X=(0Xtn/2-1e-t/22n/2Γ(n/2)dt)X'=12n/2Γ(n/2)(X)n/2-1e-X/2(X)X'=12n/2-1Γ(n/2)Xn/4-1e-X/2

nonché alcune semplificazioni per l'integrale di convoluzione ( ha il pdf χ 2 m con m dof):X2χm2

Kmn(t): =(χm2*ψn2)(t)=0tχm2(X)ψn2(t-X)dX=(2(n+m)2+1Γ(m2)Γ(n2))-10t(t-X)n4-1Xm2-1exp(-(t-X+X)/2)dX

Qualcuno vede un buon modo per calcolare questo integrale per qualsiasi t reale o deve essere calcolato numericamente? O mi sto perdendo una soluzione molto più semplice?


1
Se la non fosse quadrata, avrei qualche consiglio specifico. Non penso che questo sarà trattabile (né necessariamente particolarmente illuminante anche se dovesse dimostrarsi trattabile). Sarei tentato di esaminare approcci computazionali, come la convoluzione numerica o la simulazione, a seconda di cosa esattamente vuoi fare con il risultato. Y
Glen_b

È molto improbabile, secondo me, che l'integrale possa essere fatto.
Dave31415

1
@ Dave31415 Anche per e m , l'integrale può essere calcolato esplicitamente per valori integrali positivi di n e m . Equivarrà a una combinazione lineare di esponenziali e funzioni di errore con coefficienti che sono polinomi in nmnm . La valutazione può essere eseguita tramite la sostituzionex=t-u2. Ad esempio, conn=2,m=4otteniamo 1tx=tu2n=2,m=4. 14e18(2t+1)2(et2(2π(4t+3)(erfi(2t122)+erfi(122))+4e18)4et2+18(2t+1))
whuber

Bello. Per i numeri dispari, potresti probabilmente approssimarlo con la media del risultato per delimitare i numeri pari? O forse no.
Dave31415,

1
Grazie per le tue risposte! Per alcuni casi pari ho ottenuto un risultato simile con la funzione di Dawson, ma sembra che dovrò fare un po 'più di lavoro per una soluzione generale ...
Leo Szilard

Risposte:


2

Nel caso in cui aiuti, la variabile è una variabile gamma random generalizzata (vedi ad esempio Stacy 1962). La tua domanda è chiedere la distribuzione della somma di una variabile aleatoria chi-quadrato e una variabile aleatoria gamma generalizzata. A mia conoscenza, la densità della variabile risultante non ha espressione in forma chiusa. Quindi, la convoluzione che hai ottenuto è un integrale senza soluzione di forma chiusa. Penso che sarai bloccato con una soluzione numerica per questo.Y2


Stacy, EW (1962). Una generalizzazione della distribuzione gamma. Annali di statistiche matematiche 33 (3) , pagg. 1187-1192.


0

2


2
Potresti spiegare come questo risponde alla domanda? Non sembra direttamente correlato.
whuber

La convoluzione di Pearson di tipo III con se stessa può essere fatta. Per qualche ragione è più facile risolvere una cosa con se stessa, piuttosto che una cosa con un'altra. Ad esempio, ho risolto la convoluzione di Pearson di tipo III e ottenuto le convoluzioni di ND con GD, un problema correlato.
Carl,

Non sembra aver aiutato, eliminerà a breve.
Carl,
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.