Hello all –
I understand this has been asked, but would it make sense to run per-OTU linear models that include adjustment covariates (e.g., sex) and use the scaled residuals as input for SPLS(DA)? There are a few papers that have done something similar in the context of Spearman correlations, for instance.
Tanks very much.