» Graphics

Toulouse tutorial 6-7 Oct 2014: Details

sdejean — Wed, 17 Sep 2014 06:55:50 +0000

Training aims

The objective of this tutorial is to introduce the fundamental concepts behind projection-based approaches and illustrate their application on some exemplar studies using the R package mixOmics. In this tutorial, we will focus on the application of these approaches to medium and high throughput biological data (transcriptomics, metabolomics, proteomics data) using PCA, CCA, PLS, PLS-DA and the variants that the mixOmics team and collaborators have developed.

Training program^*

Key methodologies in mixOmics and their variants
1. Principal Component Analysis: PCA, sparse PCA, NIPALS
2. Canonical Correlation Analysis: CCA, regularized CCA
3. Partial Least Squares regression: PLS, sparse PLS
4. Partial Least Squares Discriminant Analysis: PLS-DA, sparse PLS-DA
5. Recent developments in mixOmics, including multilevel PLS for repeated measurements and introduction to the integration of multiple data sets
Review on the graphical outputs implemented in mixOmics
1. Sample plots
2. Correlation circles
3. Integrating two data sets: relevance networks and clustered image map
Case studies and applications
1. Example with PCA: Nutrimouse
2. Example with CCA: Multidrug
3. Example with PLS: Liver toxicity
4. Example with PLS-DA: SRBCT study

* The methodologies will be presented throughout the two days training session. Each methodology will be illustrated on a case study (we will alternate theory and application).

Course material will be available in a hard printed copy and online.

Other information

Details for registration are provided here.

The session will take place in the formation room of the INRA center of Toulouse-Auzeville.

Presentation about new developments in mixOmics

sdejean — Tue, 01 Oct 2013 07:17:10 +0000

Another presentation of mixOmics dating Sept 2013, which presents the latest developments including time course data, multi-group and multi-block analysis.

Go here.

Article published explaining correlation circle plots, relevance networks and CIM

Kim-Anh Lê Cao — Mon, 25 Feb 2013 05:24:28 +0000

Our manuscript ‘Insightful graphicalt outputs to explore relationships between two “omics” data sets has been published and explains how to interpret Correlation Circle plots, how relevance networks and CIM are generated from rCCA and sPLS.

Check this very colourful manuscripthere!

General presentation about mixOmics

Kim-Anh Lê Cao — Tue, 22 May 2012 02:14:31 +0000

A new general presentation about mixOmics is available (and should be updated for major update of the package) in the Presentation Section.

Lê Cao K.-A. Unravelling `omics’ data with the mixOmics R package, Illustration on several studies. General presentation on mixOmics (last updated 05/04/2012) [Presentation]

(s)IPCA

Kim-Anh Lê Cao — Fri, 06 Apr 2012 01:42:53 +0000

Independent Principal Component Analysis (IPCA)

In some case studies, we have identified some limitations when using PCA:

PCA assumes that gene expression follows a multivariate normal distribution and recent studies have demonstrated that microarray gene expression measurements follow instead a super-Gaussian distribution
PCA decomposes the data based on the maximization of its variance. In some cases, the biological question may not be related to the highest variance in the data

Instead, we propose to apply Independent Principal Component Analysis (IPCA) which combines the advantages of both PCA and Independent Component Analysis (ICA). It uses ICA as a denoising process of the loading vectors produced by PCA to better highlight the important biological entities and reveal insightful patterns in the data.

IPCA offers a better visualization of the data than ICA and with a smaller number of components than PCA.

How to choose the number of components:

The kurtosis measure is used to order the loading vectors to order the Independent Principal Components. We have shown that the kurtosis value is a good post hoc indicator of the number of components to choose, as a sudden drop in the values corresponds to irrelevant dimensions.

Sparse Independent Principal Component Analysis (sIPCA)

Similar to the sparse PCA version implemented in mixOmics, soft-thresholding is applied in the independent loading vectors in IPCA to perform internal variable selection.

How to choose the number of variables to select:

The number of variables to select is still an open issue. In our paper we proposed to use the Davies Bouldinmeasure which is an index of crisp cluster validity. This index compares the within-cluster scatter with the between-cluster separation.

More details about how to use the ipca.R function in the case study.

References

Yao F., Coquery J., Lê Cao K.-A. (2012) Independent Principal Component Analysis for biologically meaningful dimension reduction of large biological data sets, BMC Bioinformatics 13:24.
Comon P: Independent component analysis, a new concept? Signal Process 1994, 36:287-314.
Hyvärinen A, Oja E: Indepedent Component Analysis: Algorithms and Applications. Neural Networks 2000, 13(4-5):411-430

New methods: multilevel analyses

Kim-Anh Lê Cao — Thu, 05 Apr 2012 04:38:10 +0000

A multilevel approach has been added for cross-over design experiments (up to two cross factors), in collaboration with A/Prof B. Liquet (Universite de Bordeaux, France). This approach takes into account the complex structure of repeated measurements from different assays, where different treatments are applied on the same subjects to highlight the treatment effects within subject separately from the biological variation between subject.

Two different frameworks are proposed:

a discriminant analysis (method = ‘splsda’) enables the selection of features separating the different treatments
a integrative analysis (method = ‘spls’) enables the interaction of two matched data sets and the selection of subset of correlated variables (positively or negatively) across the samples. The approach is unsupervised: no prior knowledge about the samples groups is included.

The multilevel function first decomposes the variance in the data sets X (and Y) and applies either sPLS-DA or sPLS on the within-subject deviation. One or two-factor analyses are available for sPLS-DA.

Associated functions include: multilevel.R, tune.multilevel.R, pheatmap.multilevel.R (see examples in methods, graphics and case studies).

This is our first step towards repeated measurements designs.

The package has been updated to version 4.0-1 to implement these methodologies. It now requires the library ‘pheatmap’.

Web-interface

chua — Mon, 01 Aug 2011 00:00:27 +0000

R package and Methods: IPCA and sparse IPCA functions have been implemented (as well as their associated S3 functions). IPCA stands for Principal Component Analysis with Independent Loadings. It is a combination of the advantages of both PCA and Independent Component Analysis (ICA). PCA is a powerful exploratory tool if the biological question is related to the highest variance. ICA was recently proposed in the literature as an alternative to PCA as it optimizes an independence condition that can give more meaningful components. A preprint can be available upon request.
R package and Data: The Liver Toxicity study data has been updated to provide geneBank IDs and gene titles
R package and Data: Two other data sets have been added: Prostate Tumor study (gene expression) and Metabolomic study of Yeast (metabolomics).
Web interface: We are making good progress on our associated web-interface (now deployed on http://mixomics.qfab.org). Few illustrative examples are also available, and you can download the illustrative examples and run any type of analysis trough the interface. We are currently developing a ‘next level analysis’ to provide pathway enrichment analyses and give the functional annotation of the selected genes using the iHOP database. Do not hesitate to give us some feedback!
‘sletter: we now have a newsletter, to subscribe, send an email to mixomics[at]math.univ-toulouse.fr with no subject in the body.

New Graphics: network & cim

chua — Wed, 01 Jun 2011 00:00:33 +0000

New S3 method network and cim for results from PLS model
New code for the valid function to PLS-DA and SPLS-DA models validation
The S3 method plot.valid was modified to display graphical results from valid function for PLS-DA and SPLS-DA models
cim and network functions were modified to obtain the similarity matrix in return value
The S3 method plotVar was modified to obtain the coordinates for X and Y variables in return value
The predict function has been modified to simultaneously run either several or all prediction methods available to predict the classes of the test data from PLS-DA and SPLS-DA models

New Function: (s)PCA added

chua — Tue, 01 Mar 2011 00:00:52 +0000

New function pca and spca are now available to perform Principal Component Analysis (PCA) and sparse PCA for variable selection
The S3 methods plotVar, plot3dVar, plotIndiv, plot3dIndiv were modified to generate graphical results for pca and spca

New function: plot.valid

chua — Mon, 01 Nov 2010 00:00:12 +0000

New function plot.valid to display the results of the valid function
New code for imgCor function for a nicer representation of the correlation matrices
In predict function the argument 'method' were replaced by method = c("max.dist", "class.dist", "centroids.dist", "mahalanobis.dist")
The arguments dendrogram, ColSideColors and RowSideColors were added to the cim function
valid function can also been performed with missing values
Functions pls, plsda, spls and splsda were modified to identify zero- or near-zero variance predictors
The functions plotVar and plot3dVar were modified to represent only the X variables in the case of PLS-DA and SPLS-DA
The pca function has been improved so that the S3 methods plotIndiv, plot3dIndiv, plotVarand plot3dVar can be used with these new classe