auto_scale

This commit is contained in:
Ricardo Andrade 2013-04-12 18:45:14 +01:00
parent a42d84274d
commit f6b98160a7

View file

@ -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):