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)