def iX(i=i, X=X, p=p, pop=rate_model.regional_population(key)): return i * X * (1-p) * pop
def iX(i=i, X=X, p=p, pop=rate_model.regional_population(key)): birth_yld = np.zeros_like(p) birth_yld[0] = p[0] * pop[0] return i * X * (1 - p) * pop + birth_yld
def iX(i=i, X=X, p=p, pop=neg_binom_model.regional_population(key)): birth_yld = pl.zeros_like(p) birth_yld[0] = p[0] * pop[0] return i * X * (1-p) * pop + birth_yld