Thank you for the quick responses! The individual noise terms are not required but I indeed have missing values in the new data for which I hoped to use the imputation posterior all_x_new to sample. As you posted in OOS predictions with missing input values I could fill the gaps of x_with_missing by sampling them from an appropriate distribution and putting the result together with x_with_missing. That helps a lot!