New knol on the estimate of H1N1 cases in Mexico at the early stage of the pandemic
We recently published a knol in PLoS Currents Influenza about the estimate of H1N1 cases in Mexico at the early stage of the pandemic conducted with GLEaM:
Estimate of Novel Influenza A/H1N1 cases in Mexico at the early stage of the pandemic with a spatially structured epidemic model
V Colizza, A Vespignani, N Perra, C Poletto, B Gonçalves, H Hu, D Balcan, D Paolotti, W Van den Broeck, M Tizzoni, P Bajardi, JJ Ramasco. PLoS Currents: Influenza. 2009 Nov 11:RRN1129.
Determining the number of cases in an influenza epidemic represents a great challenge. Reliable figures for the actual number of cases are key to properly assess several parameters, such as mortality, morbidity or hospitalization rates. This is particularly relevant during the early phase of an outbreak, in order to better inform decision making process. Surveillance of cases inevitably suffers of several biases overall leading to an underascertainment of influenza cases. For example, many cases showing mild symptoms might not seek for medical attention, and would therefore not be included in the counting of cases. Moreover, the monitoring of cases is expected to change with time. After the initial alert at the beginning of an outbreak, an enhanced surveillance system is expected to have a large monitoring capacity. However, this situation changes with time, due to the large increase in the number of cases and the majority of resources being dedicated to the severe patients. This leads to an ascertainment mainly focused on the most severe cases, with the number of confirmed cases being a gross underestimation of the actual number of people infected by influenza. A significant example is provided by the time evolution of the monitored cases in Mexico during the first months of the epidemic: two studies   have assessed the size of the epidemic in the country and both of them found a significant underascertainment of the confirmed cases as reported by the mexican authorities. 
In order to address this point we use GLEaM to simulate the epidemic spreading and compute the number of cases in Mexico at the end of April and the beginning of May. We calibrate our model with a maximum likelihood estimate of the infection parameters, fitting the empirical arrival dates of the first infected in the countries seeded from Mexico (a detailed description of the estimate procedure is reported here). This allows us to provide an ab-initio calculation of the number of cases in Mexico independently by the estimation precedure. In the Table below the number of infectious individuals obtained in our simulations is reported, comparing it with the same number obtained in Refs  and with the confirmed cases reported by mexican authorities.
|Number of symptomatic cases in Mexico (Apr. the 30th)|
|Simulation Results||[121,000 - 1,394,000]|
|Lower bound range of Ref. ||113,000-375,000|
|Estimate of Ref. *||2,000 – 280,000|
|Mexican official report  (confirmed cases)||3,350|
Despite the different approximations used here and in Refs. , the three approaches are providing support to the possibility of a reporting ratio of infected cases in Mexico as low as 1 in 100. This finding is important when evaluating the massive amount of data which are now being collected in a large number of countries around the world. We can easily imagine that the reporting rate as well as the estimate of the cumulative attack rate in most of the countries could be easily underestimated in similar cases.