h_incidence_rate function

Helper functions for incidence rate