est_seroincidence function

Find the maximum likelihood estimate of the incidence rate parameter