From da5de6b814cfca1b8de609ac3819d7fa9b240555 Mon Sep 17 00:00:00 2001 From: Alan Saul Date: Tue, 18 Aug 2015 16:29:44 +0100 Subject: [PATCH] Added fallback for non-identity link function for Gaussian variational expectations --- GPy/likelihoods/gaussian.py | 3 +++ 1 file changed, 3 insertions(+) diff --git a/GPy/likelihoods/gaussian.py b/GPy/likelihoods/gaussian.py index 9abb8cde..86aad330 100644 --- a/GPy/likelihoods/gaussian.py +++ b/GPy/likelihoods/gaussian.py @@ -315,6 +315,9 @@ class Gaussian(Likelihood): return -0.5*np.log(2*np.pi) -0.5*np.log(v) - 0.5*np.square(y_test - mu_star)/v def variational_expectations(self, Y, m, v, gh_points=None, Y_metadata=None): + if not isinstance(self.gp_link, link_functions.Identity): + return super(Gaussian, self).variational_expectations(Y=Y, m=m, v=v, gh_points=gh_points, Y_metadata=Y_metadata) + lik_var = float(self.variance) F = -0.5*np.log(2*np.pi) -0.5*np.log(lik_var) - 0.5*(np.square(Y) + np.square(m) + v - 2*m*Y)/lik_var dF_dmu = (Y - m)/lik_var