From f6b98160a7c0ace6ca5f795aeb878d30b8aaf6a4 Mon Sep 17 00:00:00 2001 From: Ricardo Andrade Date: Fri, 12 Apr 2013 18:45:14 +0100 Subject: [PATCH] auto_scale --- GPy/models/sparse_GP.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/GPy/models/sparse_GP.py b/GPy/models/sparse_GP.py index cebcba0b..b816e684 100644 --- a/GPy/models/sparse_GP.py +++ b/GPy/models/sparse_GP.py @@ -185,7 +185,11 @@ class sparse_GP(GP): self.likelihood._set_params(p[self.Z.size+self.kern.Nparam:]) self._compute_kernel_matrices() if self.auto_scale_factor: - self.scale_factor = np.sqrt(self.psi2.sum(0).mean()*self.likelihood.precision) + if self.likelihood.is_heteroscedastic: + self.scale_factor = max(100.,(self.psi2_beta_scaled.sum(0).max())) + print self.scale_factor + else: + self.scale_factor = np.sqrt(self.psi2.sum(0).mean()*self.likelihood.precision) self._computations() def _get_params(self):