Compute the variance of the data or the latent field
ngme_cov_matrix(ngme_object, model_name = "data", replicate = 1)
if model_name = "data", then return the covariance matrix of the data (without measurement noise) if the model_name is the name or index of the latent, then return the covariance matrix of the latent field
which replicate (default = 1)
ngme_model
a data.frame of posterior samples (mesh_size * n_post_samples)