Adequately handling unbalanced groups remains one of the major challenges for the analysis of multivariate data collected from multifactorial experimental designs. While partial least squares-based methods, such as analysis of variance multiblock orthogonal partial least squares (AMOPLS), can offer better discrimination between factor levels, they can be more heavily affected by this issue, and unbalanced designs of experiments may lead to a substantial confusion of the effects. Even state-of-the-art analysis of variance (ANOVA) decomposition methodologies using general linear models (GLM) lack the ability to efficiently disentangle these sources of variation when combined with AMOPLS.