R0 function

Computes reproduction numbers from fitted models