Better tests for final size by susceptibility groups

epiverse-trace / finalsize

R package to calculate the final size of an SIR epidemic in populations with heterogeneity in social contacts and disease susceptibility

https://epiverse-trace.github.io/finalsize/

Other

11 stars 7 forks source link

Better tests for final size by susceptibility groups #43

Closed pratikunterwegs closed 2 years ago

pratikunterwegs commented 2 years ago

This PR tries to determine why multiple tests fail for final_size_grps_cpp. This is done by adding yet more tests to try and determine where failures occur. The best guess currently is numerical instability in solve_final_size_newton, which leads to NaNs and 0s returned for final sizes.