code: def predictive_quantiles(self, mu, var, quantiles, Y_metadata=None):
          
          p = self.predictive_mean(mu, var)
          return [np.asarray(p>(q/100.), dtype=np.int32) for q in quantiles]
docs: Get the "quantiles" of the binary labels (Bernoulli draws). all the
          quantiles must be either 0 or 1, since those are the only values the
          draw can take!