Detailed derivations of the reproduction numbers are also given in [Additional file 1]. Model fit using 2005-2006 Taiwan winter seasonal influenza data and agespecific vaccination data as it was implemented during that flu season, as well as simulation studies of hypothetical scenarios, will be carried out.