Create a helper function to simulate from the conditional normal distribution of new data given old data
sim_cond_normal(joint.mean, a, locs_new, locs_obs, kernel, ...)A function that takes in one argument n as the number of samples to draw from the condition normal distribution
of locs_new given locs_obs: either from rmvnorm for MVN or rnorm for univariate normal. The old and new data are assumed to follow a joint multivariate normal distribution.
The length n mean vector of the MVN distribution. By default mu1 is the first m elements of joint.mean
A vector of length n-m, the values of mu2 to condition on
A matrix containing the coordiantes of new locations
A matrix containing the coordinates of observed locations
A function (kernel function) that returns a matrix containing the similarity between the two arguments.
Hyperparameters to pass to the kernel function.
This serves as a helper function for spatialGEV_predict. The notations are consistent to the notations on the MVN wikipedia page