Finally, as an illustration in Fig 3 the prediction of the raw rate (cases × 100000/pop) during the 2010–11 flu epidemic season is provided for two counties (Vigo and Santiago) as a result of reversing the log transform of the response in the preceding models. In both counties, the peak is achieved at week 2011–5 (first week of February). The two considered horizons (t + 1 and t + 2) are shown by rows. In each case, the raw rate is compared with the prediction obtained one or two weeks before with the models LM, Rate(w); GLS–AR(1), Rate(w); GLS–AR(1), Temp(t) and GLS–AR(p), Temp(t). Focusing on t + 1, the comparison among the two dependence structures (AR(1) and AR(p), lines green and blue, respectively) associated with Temp(t) shows a big difference for Vigo but no for Santiago. This suggests that for Santiago an AR(1) is enough whereas for Vigo it seems more adequate a general AR(p) specification. Respect to the models including the Rate(t) (lines red and gray), the model using GLS reacts faster than the LM model providing better predictions of the peak. Predictions for medium or low intensities (below 125) are quite similar. For t + 2, no clear patterns are shown, although the specification GLS–AR(p) seems to do slightly better.