Time Series Forecasting, Bayesian Machine Learning.
len(y)-1 mod, samples = analysis(Y=y, X=x, family='poisson', forecast_start=forecast_start, forecast_end=forecast_end, k=k, ntrend=1, nsamps=5000, seasPeriods=[12], seasHarmComponents=[[1,2]], prior_length=4, deltrend=0.94, delregn=0.90, delVar=0.98, delSeas=0.98, rho=.6, ) forecast =