thanks for the followup - nice link collection!
Hi Guys
I am trying to read more about methods available/recommended for clustering gene expression data. Basically the input is a matrix of gene expression values from RNA-Seq differential expression analysis. In this case we have 10 different time points/ conditions.
We are trying to find out clusters which represent the states of the set of the genes across these different time points. Google/paper search did lead me to few hits but I was wondering if I can reach out to folks who have done this analysis in the past and learn from their experience.
Thanks! -Abhi
6 answers
Hey Guys
Just for completeness I found the following resources to be useful.
Concise review (somewhat old but useful) How does gene expression clustering work ? Patrik D'heaseleer : link
Data Clustering : A review by Jain, Murty And Flynn : link
Specific to gene expression clustering : Cluster Analysis for Gene Expression Data: A Survey :link
Judging the Quality of Gene Expression-Based Clustering Methods Using Gene Annotation : link
Tool: The C Clustering Library : link
Estimating the number of clusters in a data set via the gap statistic : link basically a method to determine an optimum #K(partitions) for partition based clustering
Problems in gene clustering based on gene expression data. Journal of Multivariate Analysis 90, 44–66.
Clustering tools available in R ( somewhat exhaustive list as far as the basic methods are concerned) link
A roadmap of clustering algorithms: finding a match for a biomedical applicationlink
-Abhi
I'm not sure if you are looking for something different than "normal" clustering? I mean, does the fact that this is timecourse data influence the clustering approach you want to take?
Anyway, I think Section 7 of the DESeq vignette could be helpful, as they first shoot their data through a variance stabilized transformation before doing exploratory analysis (heatmaps, PCA, etc). which (by the looks of what they present there) seem like a good thing to try.
Yes, the best bet for now is probably to use DESeq's variance-stabilizing transformation and then use tools designed for microarray analysis.
Sean Thomas developed a procedure called rank expectation that he used to cluster chromatin accessibility in areas of the Drosophila genome at five stages of development (see Fig. 4). The test procedure is generally applicable to situations where you want to discern groups or clusters of data by signal, like gene expression, chromatin accessibility, etc.
you need too set your data in 10 time point like this way 0-1 0-2 0-3 0-4 0-5 ... you can do it even with deseq2 and get output for other analysis
Log in to answer this question.