To compare any two or all three of the described methods it would be ideal to generate data sets with increasing deviations from each method’s assumptions while keeping sample and effect sizes constant (i.e. deviations from linear association, zero within-cluster covariance and multivariate normality, along with model misspecifications of LVMMs). The different methods can then be applied to the generated data, permitting a structured comparison.