Variational Gibbs inference for statistical model estimation from incomplete data