This will be a breaking release that breaks up the contact_matrix() function and, in doing so, changes the workflow for extracting contact matrices. Data wrangling will be moved into separate functions that will be expected to be called before contact_matrix()