Gibbsovo vzorkování

Gibbsovo vzorkování (Gibbsův výběrový plán, Šablona:Vjazyce2) je algoritmus, který umožňuje generovat náhodný výběr z mnohorozměrného rozdělení pomocí jednorozměrných plně podmíněných rozdělení. Tento algoritmus patří do třídy MCMC algoritmů, které náhodný výběr generují jako realizaci Markovova řetězce.

Implementace

Generování z rozdělení mnohorozměrné náhodné veličiny $𝐗=(X_1,\dots ,X_d)$ probíhá následujícím způsobem:

1. Zvolíme počáteční hodnotu ${𝐗}^{(0)}=(x_1^{(0)},\dots ,x_d^{(0)})$ a položíme $t=0$.
2. Začneme generováním prvního prvku $x_1^{(t+1)}$ z podmíněného rozdělení $X_1|x_2^{(t)},\dots ,x_d^{(t)}$. Pro druhý prvek podmiňujeme i pomocí nově generované hodnoty, tedy prvek $x_2^{(t+1)}$ je generován z podmíněného rozdělení $X_2|x_1^{(t+1)},x_3^{(t)},\dots ,x_d^{(t)}$. Postupně generujeme ostatní prvky až do $x_d^{(t+1)}$ z rozdělení $X_d|x_1^{(t+1)},\dots ,x_{d-1}^{(t+1)}$. Nové hodnoty můžeme složit do vektoru ${𝐗}^{(t+1)}=(x_1^{(t+1)},\dots ,x_d^{(t+1)})$.
3. Dokud není splněna předem stanovená podmínka (například počet iterací), položíme $t=t+1$ a vrátíme se ke druhému kroku.

Využití

Gibbsovo vzorkování je užitečné především v hierarchických bayesovských modelech, ve kterých jsou tvary plně podmíněných rozdělení přímočaré, ale úplné rozdělení parametrů může být komplikované. Počáteční hodnotu pro Gibbsovo vzorkování lze získat například pomocí EM algoritmu.

Modifikace

Gibbsovo vzorkování může být upraveno následovně:

• Místo pouze jednorozměrných plně podmíněných rozdělení můžeme používat i kombinaci vícerozměrných podmíněných rozdělení (např. $X_1|x_2,x_3$ a $X_2,X_3|x_1$).
• Místo generování od prvního prvku k poslednímu můžeme generovat v opačném pořadí; takový postup lze použít v důkazu existence limitního rozdělení.
• Index prvku, který aktualizujeme, může být vybrán náhodně (každý index má pravděpodobnost výběru $\frac{1}{d}$). V takovém případě mluvíme o náhodném procházení.

Literatura

• Šablona:Citation
• Šablona:Citation