Added code to calculate expected number of observed cases with change in gamma