Variational Inference in PRIMO