Another bug fixed

This commit is contained in:
Ricardo 2014-10-03 12:02:58 +01:00
parent 45f8751769
commit c1d998e272

View file

@ -320,7 +320,6 @@ class HalfT(Prior):
return "hT(" + str(np.round(self.A)) + ', ' + str(np.round(self.nu)) + ')' return "hT(" + str(np.round(self.A)) + ', ' + str(np.round(self.nu)) + ')'
def lnpdf(self,theta): def lnpdf(self,theta):
stop
return (theta>0) * ( self.constant -.5*(self.nu+1) * np.log( 1.+ (1./self.nu) * (theta/self.A)**2 ) ) return (theta>0) * ( self.constant -.5*(self.nu+1) * np.log( 1.+ (1./self.nu) * (theta/self.A)**2 ) )
#theta = theta if isinstance(theta,np.ndarray) else np.array([theta]) #theta = theta if isinstance(theta,np.ndarray) else np.array([theta])