modelSelection_Q function

Selects the number of groups with ICL