Skip to content

Instantly share code, notes, and snippets.

@jcbozonier
Last active May 5, 2017 12:54
Show Gist options
  • Save jcbozonier/edc74d313eb882709c12e7b9295b9cb6 to your computer and use it in GitHub Desktop.
Save jcbozonier/edc74d313eb882709c12e7b9295b9cb6 to your computer and use it in GitHub Desktop.
Setting up our global priors
import pymc3 as pm
with pm.Model() as model:
global_m = pm.Normal('global_m', mu=0, sd=100**2)
global_m_sd = pm.Uniform('global_m_sd', lower=0, upper=1000)
global_b = pm.Normal('global_b', mu=0, sd=100**2)
global_b_sd = pm.Uniform('global_b_sd', lower=0, upper=1000)
degree_m = pm.Normal('degree_m', mu=global_m, sd=global_m_sd, shape=degree_count)
degree_m_sd = pm.Uniform('degree_m_sd', lower=0, upper=1000, shape=degree_count)
degree_b = pm.Normal('degree_b', mu=global_b, sd=global_b_sd, shape=degree_count)
degree_b_sd = pm.Uniform('degree_b_sd', lower=0, upper=1000, shape=degree_count)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment