I'm conducting LCGA to explore whether some classes in my sample are different on distal outcome.
In preliminary analysis I've identified one significant covariate likely to impact my distal outcome and I'd like to include it in comparisons tests between classes.
On the basis on the recommandation of Jung & Wickrama (2008), the first step is to choose the best-fitting model (based on fit indexes), by performing a series of models with an increasing number of classes.
I'm wondering if it's right to perform these series of models including the covariate. Or should I perform a series of models with an increasing number of classes without covariate, choose the best-fitting model with X classes, and after perform the chosen model with X classes and include at this moment the covariate in the model to explore difference between classes on distal outcome?
Does anybody have any suggestion? Or references using this type of data analysis with covariate?
Many thanks