code: def predict_quantiles(self, Xnew=None, quantiles=(2.5, 97.5), balance=None, **kw):
          
          if balance is None:
              p_balance = self.balance
          else:
              p_balance = balance
          
          
          mu, var = self._raw_predict(Xnew, p_balance=p_balance)
          #import pdb; pdb.set_trace()
          return  [stats.norm.ppf(q/100.)*np.sqrt(var + float(self.Gaussian_noise.variance)) + mu for q in quantiles]
docs: Inputs:
          ------------------
          
          balance: bool
          Whether to balance or not the model as a whole