def predictive_quantiles(self, mu, var, quantiles, Y_metadata=None):
        return  [stats.norm.ppf(q/100.)*np.sqrt(var + self.variance) + mu for q in quantiles]
dataset/ETHPy150Open SheffieldML/GPy/GPy/likelihoods/gaussian.py/Gaussian.predictive_quantiles/original