This question concerns Gelman’s Rat Tumor example in BDA 3 (bottom of page 109). I’m trying to replicate the inference made in this example. Here is some code to get started: import numpy as np import pandas as pd import pymc3 as pm import theano.tensor as T import matplotlib.pyplot as plt # rat…

Actually, the problem is because you didnt model a vector of theta: with pm.Model() as model: # Uninformative prior for alpha and beta u = pm.Uniform('u', 0, 1) v = pm.Uniform('v', 0, 1) alpha = pm.Deterministic('a', u / tt.pow(v, 2)) beta = pm.Deterministic('b', (1 - u) / …

Reparametizing Gelman's rat tumour example

Questions

junpenglao June 24, 2018, 8:39am 2

Hmm, if I understand correctly what BDA described is uniform on the domain of the parameter, not just on [0, 1]

Topic		Replies	Views
Gelman's (alpha + beta)^-5/2 parametrization for beta prior in AB test Questions	7	1092	April 3, 2021
Apparent pathological behavior Questions	6	875	August 10, 2018
Advice on Reparameterization Questions	2	816	December 10, 2018
Conversion rate model: Beta parameters from pm.Uniform Questions	5	1154	February 7, 2021
What to do with divergences in my case Questions	7	2495	September 21, 2017

Reparametizing Gelman's rat tumour example

Related topics