Use the function predict to compute 95% prediction intervals for the 30 observed values using
- the model where we only regress on hardness
- the model where we regress on hardness and tension. Compare the prediction intervals. It can be useful to plot the prediction intervals versus the fitted values.