Given mean and covariance matrix of a multivariate normal distribution plus a part of observations, I really need to generate the conditional normal given these distributions. The most ideal package I have found is r-condMVNorm, which contains the method rcmvnorm (see for example https://rdrr.io/cran/condMVNorm/man/cmvnorm.html) doing exactly what I want.