This paper presents simple models developed to generate covariances between observable and latent variables. The methodology consists of using “variance penalty” terms as a measure of the contribution of the latent-variable uncertainties to the variance of a given calculated quantity z. This approach provides a useful understanding of how the observable and latent variables are related to each other and ensures the positive-definiteness of the covariance matrix. This work has been implemented in the nuclear data assimilation tool CONRAD. Performances of analytic and Monte Carlo models are illustrated with covariances calculated for neutron-induced capture reactions on stable xenon isotopes (124Xe, 126Xe, 128Xe, 129Xe, 130Xe, 132Xe, and 134Xe).