Miscellaneous

  • How to cite mixOmics?

It depends… See here to cite the paper corresponding to your need.

  • How to extract the names of variables selected when performing sPLS?

We will update a new version to extract theses genes in a most explicit manner. For now it is, for example for component 1 to extract the genes selected from data set X:

comp = 1
which(result$loadings$X[, comp] != 0)

You can do the same for component 2 and 3 and also for the data set Y.
Beware that the absolute value of these loadings indicate the importance of these genes so you can actually sort your list:

loadingsX1 = abs(result$loadings$X[, comp])
names(sort(loadingsX1, decreasing = T)[1:50])

will give you the names of the genes sorted by importance and selected by sPLS.

  • Is there a way the information on genes and metabolites connections can be incorporated in the mixOmics analysis?

mixOmics approaches are entirely data-driven, and we do not include any priori information. However, networks can be inferred by using similarity matrices that are a direct product of the methods RCCA or sPLS. So you can have a look at the network or cim functions (see more details on our website). The web
interface will also be able to give you such outputs.

  • Why does plsda return an error message even though I am using the same code as indicated in the tutorial or the help file?

This is highly likely due to a clash with another package that also provides a plsda function (caret or muma for instance; the spls package provides a splsda function) that was loaded before mixOmics. You can either unload these packages using the function detach() or use the command mixOmics::plsda(…) to be force the plsda function from mixOmics to be run.