def estimate_with_fixed_beta(self, data, beta): mu = median(data) v = mean((data - mu) ** 2) alpha = sqrt(v * exp(gammaln(1. / beta) - gammaln(3. / beta))) return mu, alpha