def logp (value, tree, nregime): locs = np.empty(nregime, dtype=object) for reg in range(nregime): locs[reg] = [i+1 for i,v in enumerate(value) if v==reg] nswitches = nshifts(tree, value) # WHAT SHOULD BE THE PRIOR ON NUMBER OF SWITCHES? return pymc.exponential_like(nswitches, beta=1)
def l(value=0.1, rate=1.0): """Rate stochastic of poisson distribution.""" return exponential_like(value, rate)
def l(value=.1, rate=1.): """Rate stochastic of poisson distribution.""" return exponential_like(value, rate)