regional_mean works for 2-D and 3-D with `group` method