p_group_by function

Stratifying your analysis