def log_prob(self, x): mu = self.mu alpha = self.alpha beta = self.beta return log(beta / (2.0 * alpha)) - gammaln(1. / beta) - power(fabs(x - mu) / alpha, beta)