where μ i is the age-dependent base mortality probability of infected i, ρ i is her status of antiviral therapy (0 or 1), and τ is the antiviral efficacy measured as the relative decrease in the base probability [44] . We assume that a recovered case develops full immunity but continues circulating in the region.