Authors: Andrea Rau, Cathy Maugis-Rabusseau, and Antoine Godichon-Baggioni
Co-expression analysis for expression profiles arising from high-throughput sequencing data. Feature (e.g., gene) profiles are clustered using adapted transformations and mixture models or a K-means algorithm, and model selection criteria (to choose an appropriate number of clusters) are provided.
A typical call to coseq to apply the K-means algorithm to logCLR-transformed normalized RNA-seq profiles takes the following form:
library(coseq)
## The following two lines are equivalent:
run_kmeans <- coseq(counts, K=2:10)
run_kmeans <- coseq(counts, K=2:10, model="kmeans", transformation="logclr")
where counts represents a (n×q) matrix or data frame of read counts for n genes in q samples
and K=2:10 provides the desired range of numbers of clusters (here, 2 to 10).
A typical call to coseq to fit a Gaussian mixture model on arcsin- or logit-transformed normalized RNA-seq profiles takes the following form:
run_arcsin <- coseq(counts, K=2:10, model="Normal", transformation="arcsin")
run_logit <- coseq(counts, K=2:10, model="Normal", transformation="logit")
We note that this function directly calls the Rmixmod R package to fit Gaussian mixture models.
The output of the coseq function is an
S4 object of class coseqResults on which standard plot and summary functions can be directly applied; the former
uses functionalities from the ggplot2 package. The option of parallelization
via the BiocParallel Bioconductor package is also provided.
Rau, A. and Maugis-Rabusseau, C. (2018) Transformation and model choice for co-expression analysis of RNA-seq data. Briefings in Bioinformatics, 19(3)-425-436.
Godichon-Baggioni, A., Maugis-Rabusseau, C. and Rau, A. (2018) Clustering transformed compositional data using K-means, with applications in gene expression and bicycle sharing system data. Journal of Applied Statistics, doi:10.1080/02664763.2018.1454894.
The coseq package is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License, version 3, as published by the Free Software Foundation.
This program is distributed in the hope that it will be useful, but without any warranty; without even the implied warranty of merchantability or fitness for a particular purpose. See the GNU General Public License for more details.
A copy of the GNU General Public License, version 3, is available at http://www.r-project.org/Licenses/GPL-3.