def predictive_quantiles(self, mu, var, quantiles, Y_metadata=None):
p = self.predictive_mean(mu, var)<EOL>return [np.asarray(p>(q/<NUM_LIT>), dtype=np.int32) for q in quantiles]<EOL>
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!