at some point, this should plot the predictive distribution rather than the 95% error interval. See http://nbviewer.ipython.org/5982623 for code and example of what I mean