@as_op(itypes=[tt.lscalar, tt.dscalar, tt.dscalar], otypes=[tt.dvector])
def rateFunc(switchpoint, early_mean, late_mean):
 out = empty(years)
 out[:switchpoint] = early_mean
 out[switchpoint:] = late_mean
 return out

@ as_op ( itypes = [ tt . lscalar , tt . dscalar , tt . dscalar ] , otypes = [ tt . dvector ] ) def rateFunc ( switchpoint , early_mean , late_mean ) : out = empty ( years ) out [ : switchpoint ] = early_mean out [ switchpoint : ] = late_mean return out
concatenate poisson means .