Gibbs Sampling with Missing Data

Gibbs sampling with missing data treats unobserved values as additional unknowns alongside model parameters and samples all of them jointly within a Markov chain Monte Carlo loop. The method alternates between drawing the missing values from their conditional distribution given the parameters and drawing the parameters from their conditional distribution given the completed data, producing a posterior over both simultaneously.