Concurrent ordination: Simultaneous unconstrained and constrained latent variable modelling

In community ecology, unconstrained ordination can be used to indirectly explore drivers of community composition, while constrained ordination can be used to directly relate predictors to an ecological community. However, existing constrained ordination methods do not explicitly account for community composition that cannot be explained by the predictors, so that they have the potential to misrepresent community composition if not all predictors are available in the data. We propose and develop a set of new methods for ordination and joint species distribution modelling (JSDM) as part of the generalized linear latent variable model (GLLVM) framework, that incorporate predictors directly into an ordination. This includes a new ordination method that we refer to as concurrent ordination, as it simultaneously constructs unconstrained and constrained latent variables. Both unmeasured residual covariation and predictors are incorporated into the ordination by simultaneously imposing reduced rank structures on the residual covariance matrix and on fixed-effects. We evaluate the method with a simulation study, and show that the proposed developments outperform canonical correspondence analysis (CCA) for Poisson and Bernoulli responses, and perform similar to redundancy analysis (RDA) for normally distributed responses, the two most popular methods for constrained ordination in community ecology. Two examples with real data further demonstrate the benefits of concurrent ordination, and the need to account for residual covariation in the analysis of multivariate data. This article contextualizes the role of constrained ordination in the GLLVM and JSDM frameworks, while developing a new ordination method that incorporates the best of unconstrained and constrained ordination, and which overcomes some of the deficiencies of existing classical ordination methods.