Model estimation was conducted for Main effect only (Model 1)
and Interaction effect (Model 2) in order to obtain model-fit indices.
Then, in order to select the best model, AIC and BIC indices were
compared because the likelihood ratio test tends to be unreliable
when using mixture analysis.