Biplot from PLSDA

Hi @enzo,

Thanks for the code. We have in fact implemented a customisable biplot function which is still in development :slight_smile:.

Below is an example with installation instructions:

BiocManager::install('ajabadi/mixOmics@devel', ask=FALSE, update = FALSE)
library(mixOmics)
data(breast.tumors)
X <- breast.tumors$gene.exp
colnames(X) <- paste0('GENE_', colnames(X))
rownames(X) <- paste0('SAMPLE_', rownames(X))
Y <- breast.tumors$sample$treatment

plsda.breast <- plsda(X, Y, ncomp = 2)
biplot(plsda.breast, cutoff = 0.72)

## remove arrows
biplot(plsda.breast, cutoff = 0.72, var.arrow.col = NULL, var.names.size = 4)

Created on 2020-10-27 by the reprex package (v0.3.0)

You can see more examples a ?mixOmics::biplot. Iā€™d appreciate your feedback on what I could further improve :slight_smile: .

Hope it helps,

Al

1 Like