mixOmics 6.23.4

- Preamble
- 1 Introduction {#01}
- 2 Let’s get started {#02}
- 3 PCA on the
`multidrug`

study {#03}- 3.1 Load the data {#03:load-data}
- 3.2 Example: PCA {#03:pca}
- 3.2.1 Choose the number of components {#03:pca-ncomp}
- 3.2.2 PCA with fewer components {#03:pca-final}
- 3.2.3 Identify the informative variables {#03:pca-vars}
- 3.2.4 Sample plots {#03:pca-sample-plot}
- 3.2.5 Variable plot: correlation circle plot {#03:pca-variable-plot}
- 3.2.6 Biplot: samples and variables {#03:pca-biplot}

- 3.3 Example: sparse PCA {#03:spca}

- 4 PLS on the liver toxicity study {#04}
- 5 PLS-DA on the SRBCT case study {#05}
- 6 N-Integration {#06}
- 6.1 Block sPLS-DA on the TCGA case study {#06:diablo}
- 6.2 Load the data {#06:diablo-load-data}
- 6.3 Parameter choice {#06:diablo-params}
- 6.4 Final model {#06:diablo-final}
- 6.5 Sample plots {#06:diablo-sample-plots}
- 6.6 Variable plots {#06:diablo-variable-plots}
- 6.7 Model performance and prediction {#06:diablo-perf}

- 7 P-Integration {#07}
- 8 Session Information {#08}
- References

If you are following our online course, the following vignette will be useful as a complementary learning tool. This vignette also covers the essential use cases of various methods in this package for the general `mixOmcis`

user. The below methods will be covered:

- (s)PCA,
- PLS1 and PLS2,
- (s)PLS-DA,
- N-integration (multi-block sPLS-DA, aka. “DIABLO”), and
- P-integration (multi-group sPLS-DA, aka “MINT”).

As outlined in 1.3, this is not an exhaustive list of all the methods found within `mixOmics`

. More information can be found at our website and you can ask questions via our discourse forum.

`mixOmics`

is an R toolkit dedicated to the exploration and integration of biological data sets with a specific focus on variable selection. The package currently includes more than twenty multivariate methodologies, mostly developed by the `mixOmics`

team (see some of our references in 1.2.3). Originally, all methods were designed for omics data, however, their application is not limited to biological data only. Other applications where integration is required can be considered, but mostly for the case where the predictor variables are continuous (see also 1.1).

In `mixOmics`

, a strong focus is given to graphical representation to better translate and understand the relationships between the different data types and visualize the correlation structure at both sample and variable levels.

Note the data pre-processing requirements before analysing data with `mixOmics`

:

**Types of data**. Different types of biological data can be explored and integrated with`mixOmics`

. Our methods can handle molecular features measured on a continuous scale (e.g. microarray, mass spectrometry-based proteomics and metabolomics) or sequenced-based count data (RNA-seq, 16S, shotgun metagenomics) that become `continuous’ data after pre-processing and normalisation.**Normalisation**. The package does not handle normalisation as it is platform-specific and we cover a too wide variety of data! Prior to the analysis, we assume the data sets have been normalised using appropriate normalisation methods and pre-processed when applicable.**Prefiltering**. While`mixOmics`

methods can handle large data sets (several tens of thousands of predictors), we recommend pre-filtering the data to less than 10K predictor variables per data set, for example by using Median Absolute Deviation (Teng et al. 2016) for RNA-seq data, by removing consistently low counts in microbiome data sets (Lê Cao et al. 2016) or by removing near-zero variance predictors. Such step aims to lessen the computational time during the parameter tuning process.**Data format**. Our methods use matrix decomposition techniques. Therefore, the numeric data matrix or data frames have \(n\) observations or samples in rows and \(p\) predictors or variables (e.g. genes, proteins, OTUs) in columns.**Covariates**. In the current version of`mixOmics`

, covariates that may confound the analysis are not included in the methods. We recommend correcting for those covariates beforehand using appropriate univariate or multivariate methods for batch effect removal. Contact us for more details as we are currently working on this aspect.

We list here the main methodological or theoretical concepts you need to know to be able to efficiently apply `mixOmics`

:

**Individuals, observations or samples**: the experimental units on which information are collected, e.g. patients, cell lines, cells, faecal samples etc.**Variables, predictors**: read-out measured on each sample, e.g. gene (expression), protein or OTU (abundance), weight etc.**Variance**: measures the spread of one variable. In our methods, we estimate the variance of components rather that variable read-outs. A high variance indicates that the data points are very spread out from the mean, and from one another (scattered).**Covariance**: measures the strength of the relationship between two variables, i.e. whether they co-vary. A high covariance value indicates a strong relationship, e.g. weight and height in individuals frequently vary roughly in the same way; roughly, the heaviest are the tallest. A covariance value has no lower or upper bound.**Correlation**: a standardized version of the covariance that is bounded by -1 and 1.**Linear combination**: variables are combined by multiplying each of them by a coefficient and adding the results. A linear combination of height and weight could be \(2 * weight - 1.5 * height\) with the coefficients \(2\) and \(-1.5\) assigned with weight and height respectively.**Component**: an artificial variable built from a linear combination of the observed variables in a given data set. Variable coefficients are optimally defined based on some statistical criterion. For example in Principal Component Analysis, the coefficients of a (principal) component are defined so as to maximise the variance of the component.**Loadings**: variable coefficients used to define a component.**Sample plot**: representation of the samples projected in a small space spanned (defined) by the components. Samples coordinates are determined by their components values or scores.**Correlation circle plot**: representation of the variables in a space spanned by the components. Each variable coordinate is defined as the correlation between the original variable value and each component. A correlation circle plot enables to visualise the correlation between variables - negative or positive correlation, defined by the cosine angle between the centre of the circle and each variable point) and the contribution of each variable to each component - defined by the absolute value of the coordinate on each component. For this interpretation, data need to be centred and scaled (by default in most of our methods except PCA). For more details on this insightful graphic, see Figure 1 in (González et al. 2012).**Unsupervised analysis**: the method does not take into account any known sample groups and the analysis is exploratory. Examples of unsupervised methods covered in this vignette are Principal Component Analysis (PCA, Chapter 3), Projection to Latent Structures (PLS, Chapter 4), and also Canonical Correlation Analysis (CCA, not covered here but see the website page).**Supervised analysis**: the method includes a vector indicating the class membership of each sample. The aim is to discriminate sample groups and perform sample class prediction. Examples of supervised methods covered in this vignette are PLS Discriminant Analysis (PLS-DA, Chapter 5), DIABLO (Chapter 6) and also MINT (Chapter 7).

If the above descriptions were not comprehensive enough and you have some more questions, feel free to explore the glossary on our website.

Here is an overview of the most widely used methods in `mixOmics`

that will be further detailed in this vignette, with the exception of rCCA. We depict them along with the type of data set they can handle.