prior = pymc3.Gamma('prior', alpha=1, beta=0.99)
likelihood = pymc3.Pareto('likelihood', alpha=prior, m=0.1, observed=play_data)