The parameters of the recovery model were estimated by backward stepwise likelihood ratio test in proportional hazards model. Censored patient was considered as not recovered. Twosided P values and 95% confidence interval of the parameters were calculated. The cumulative resolution rate was analyzed using Kaplan-Meier analysis and log-rank test and compared between subgroups of patients. The subgroups were based on gender, time from onset to hospital admission, general symptoms or respiratory symptoms. Geo-statistical maps were used to characterized the spatial and temporal distribution of weekly confirmed cases during May 7 th and August 12 th . Longitude and latitude coordinates of each province were identified from a national GIS database. The procedure was fulfilled in a GIS (ArcInfo 9.2, ESRI, Redlands, CA) by overlaying a national vector map.