- How to cite mixOmics?
It depends… See http://perso.math.univ-toulouse.fr/mixomics/presentations/publications/ 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.