est_foi function

Estimate force of infection