Variational learning for latent Gaussian model of discrete data