I condizionali completi possono determinare la distribuzione articolare?


9

Ho sentito che tutti i condizionali completi (usati nel campionamento di Gibbs) possono determinare la distribuzione congiunta. Ma non capisco perché e come. O ho sentito male? Grazie!

Risposte:


11

Questa domanda apparentemente semplice è più profonda di quanto sembri, portandoci fino al teorema di Hammersley-Clifford. Il fatto che possiamo recuperare la distribuzione congiunta dai condizionali completi è ciò che rende possibile il campionatore Gibbs. Può essere visto come un risultato sorprendente, se ricordiamo che i marginali non determinano la distribuzione congiunta.

Vediamo cosa succede se calcoliamo formalmente con le definizioni ben note delle densità articolari, condizionali e marginali. Da

fX,Y(X,y)=fX|Y(X|y)fY(y)=fY|X(y|X)fX(X),
noi abbiamo
fY|X(y|X)fX|Y(X|y)dy=fY(y)fX(X)dy=1fX(X),
e possiamo recuperare formalmente la densità articolare dalla creazione di tutti i condizionali
fX,Y(X,y)=fY|X(y|X)fY|X(y|X)/fX|Y(X|y)dy.(*)

Il problema con questo calcolo formale è che suppone l'esistenza di tutti gli oggetti coinvolti.

Ad esempio, considera cosa succede se ci viene dato

X|Y=y~Exp(y)eY|X=X~Exp(X).
Ne consegue che fY|X(y|X)/fX|Y(X|y)=X/ye l'integrale nel denominatore di (*) diverge.

Per garantire che possiamo recuperare la densità articolare dai condizionali completi utilizzando (*) abbiamo bisogno delle condizioni di compatibilità discusse in questo documento:

"Distribuzioni condizionate compatibili", Barry C. Arnold e S. James Press, Journal of American Statistical Association, Vol. 84, n. 405 (1989), pagg. 152-156.

Infine, leggi la discussione sul teorema di Hammersley-Clifford nel libro di Robert e Casella


1
Potresti chiarire cosa si intende per "l'integrale .... esiste"? Sembra che ci siano due diversi problemi qui, vale a dire. (i) fa l'integrale
fY|X(y|X)fX|Y(X|y)dy
esiste o no? e (ii) se l'integrale esiste, è il suo valore1fX(X)? O lo dici ogni voltaX e Y hanno densità condizionali tali che
fY|X(y|X)fX|Y(X|y)dy
esiste, quindi deve essere che il valore dell'integrale sia 1fX(X)?
Dilip Sarwate

Grazie @Zen! fY e fX|Y può determinare fX,Y, e fY|X e fX|Y può anche determinare fX,Y. (1) Quale fornisce ulteriori informazioni,fY o fY|X? (2) Quale fornisce meno informazioni ridondanti / sovrapposte fX|Y, fY o fY|X? (3) FuorifY e fY|X, uno di loro fornisce già le informazioni dell'altro (di cui dubito, perché ciò implicherebbe che l'uno conduca all'altro)? Immagino sia l'intersezione tra le informazioni di fY e di fY|X, che insieme a fX|Y determina fX,Y.
Tim

Ciao @Tim. fY ti rappresenta incertezza Y, mentre fY|X rappresenta la tua incertezza Y, dato che conosci il valore di X. "Quale contiene più informazioni?" non è una domanda facile. SefX|Y e fY|X sono compatibili (nel senso di Arnold e Press), quindi determinano fX,Y attraverso (*).
Zen

Attualmente sto lottando con lo stesso problema. Sono un po 'confuso dalla necessità di distribuzioni condizionate compatibili, poiché queste non sono mai menzionate in nessuna (almeno quelle che ho letto) introduzioni al Gibbs Sampling. O vale la necessità di distribuzioni condizionali compatibili solo se si cerca di recuperare formalmente le distribuzioni congiunte, ad esempio da (*). -> non si avvicina alla distribuzione congiunta di Gibbs Sampling?
sklingel,

In una normale impostazione di campionamento di Gibbs applicata a un problema statistico, si presume che esista la distribuzione della probabilità articolare (posteriore), quindi che i condizionali completi derivati ​​da questa distribuzione articolare sono compatibili. Al di fuori di questo caso, il campionamento di Gibbs non ha senso.
Xi'an,
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.