This is also a pretty comprehensive set of examples: Forecasting with Structural AR Timeseries — PyMC example gallery