This question is inspired by can't use emmeans inside map, and related to Map `joint_tests` to a list after fitting a `gls` model and `group_by` and keep grouping levels
