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)
