def plot_beta(alpha, beta, linestyle, **params):
pl.plot(pi, [pl.exp(mc.beta_like(pi_i, alpha, beta)) for pi_i in pi],
label=r'$\alpha$ $=$ $%.1f$, $\beta$ $=$ $%.1f$'%(alpha, beta),
linewidth=2,
alpha=.9,
linestyle=linestyle,
color='black',
**params)
def plot_beta(alpha, beta, linestyle, **params):
pl.plot(pi, [pl.exp(mc.beta_like(pi_i, alpha, beta)) for pi_i in pi],
label=r'$\alpha$ $=$ $%.1f$, $\beta$ $=$ $%.1f$' % (alpha, beta),
linewidth=2,
alpha=.9,
linestyle=linestyle,
color='black',
**params)
def in_factor(a=alpha_in, b=beta_in, p=in_prob):
"""The 'observation' that the range fills the entire EO region"""
return pm.beta_like(p,a,b)
def out_factor(a=alpha_out, b=beta_out, p=out_prob):
"""The 'observation' that the range does not extend outside the EO region"""
return pm.beta_like(p,a,b)
