def f(sp_sub, a, b, n=n):
p = pm.invlogit(sp_sub)
h = pm.rbeta(a,b,size=len(sp_sub))
p_def = g6pd.p_fem_def(p,h)
return pm.rbinomial(n=n, p=p)
def f(sp_sub, a, b, n=n):
p = pm.invlogit(sp_sub)
h = pm.rbeta(a, b, size=len(sp_sub))
p_def = g6pd.p_fem_def(p, h)
return pm.rbinomial(n=n, p=p)
def f(sp_sub, x, a, b):
p = pm.invlogit(sp_sub(x))
h = pm.rbeta(a,b,size=len(p))
return g6pd.p_fem_def(p,h)
def f(sp_sub, x, a, b):
p = pm.invlogit(sp_sub(x))
h = pm.rbeta(a, b, size=len(p))
return g6pd.p_fem_def(p, h)
