Methods for exploring treatment effect heterogeneity in subgroup