incidence_rate function

Incidence rate estimation