confidence interval for the force of infection

[choisy]

Would it be possible that the plot() method of the polynomial_model class shows the confidence interval of the force of infection as well? That would be so neat! (should be optional and specified by an argument of the plot() method).