lmdme: Linear Model on Designed Multivariate Experiments in R.

FRESNO, C.; BALZARINI, M.; FERNANDEZ, E.

JOURNAL OF STATISTICAL SOFTWARE

JOURNAL STATISTICAL SOFTWARE

Año: 2014 vol. 56 p. 1 - 16

1548-7660

The lmdme package implements analysis of variance (ANOVA) decomposition through linear models on designed multivariate experiments in R (R Development Core Team 2012), allowing ANOVA-principal component analysis (APCA) and ANOVA-simultaneous component analysis (ASCA). It also extends both methods with the application of partial least squares (PLS) through the specication of a desired output matrix. The package is freely available on the Bioconductor website (Gentleman et al. 2004), licenced under GNU general public license. ANOVA decomposition methods for multivariate designed experiments are becoming popular in omics" experiments (transcriptomics, metabolomics, etc.) where measure-ments are performed according to a predened experimental design (Smilde et al. 2005), with several experimental factors or including subject specic clinical covariates, such as those present in current clinical genomic studies. ANOVA-PCA and ASCA are well-suited methods to study interaction patterns on multidimensional datasets. However, current R implementation of APCA is only available for Spectra data (ChemoSpec), meanwhile ASCA (Nueda et al. 2007) is based on average calculations over the indexes of up to three design matrices. Thus, no statistical inference over estimated eects is provided. Moreover, ASCA is not available in R package format. Here, R implementation on ANOVA decomposition with PCA/PLS analysis is provided. It allows a flexible formula interface the specication on almost any linear model with appropriate inference over the estimated effects and display functions for both PCA and PLS. We will present the model, implementation and two high-throughput microarray examples: one applied on interaction pattern analysis and the other for quality assessment.