From 8d98652e8ba401cf20ce8dd9e159830a43c45125 Mon Sep 17 00:00:00 2001 From: Ricardo Date: Mon, 17 Mar 2014 10:29:25 +0000 Subject: [PATCH] Re-definition of the week --- GPy/likelihoods/likelihood.py | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/GPy/likelihoods/likelihood.py b/GPy/likelihoods/likelihood.py index 3eafedb1..f8fc17d2 100644 --- a/GPy/likelihoods/likelihood.py +++ b/GPy/likelihoods/likelihood.py @@ -387,13 +387,14 @@ class Likelihood(Parameterized): return pred_mean, pred_var - def predictive_quantiles(self, mu, var, quantiles, Y_metadata): + def predictive_quantiles(self, mu, var, quantiles, Y_metadata=None): #compute the quantiles by sampling!!! N_samp = 1000 s = np.random.randn(mu.shape[0], N_samp)*np.sqrt(var) + mu - ss_f = s.flatten() - ss_y = self.samples(ss_f) - ss_y = ss_y.reshape(mu.shape[0], N_samp) + #ss_f = s.flatten() + #ss_y = self.samples(ss_f, Y_metadata) + ss_y = self.samples(s, Y_metadata) + #ss_y = ss_y.reshape(mu.shape[0], N_samp) return [np.percentile(ss_y ,q, axis=1)[:,None] for q in quantiles]