This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to make cluster from gene read count?

Hi

I have several questions and I appreciate for answer :

1- I wanted to know how can I add cluster analysis like below image based on read count. 2- Should I run the code only for the genes with minimum number of reads (e.g rowSum(gene) > 10) or I should run it only for DEGs ?

Tnx.

enter image description here

rna-seq clustering r deseq2 edger
• 3,273 views • link
updated by Chirag Parsania Chirag Parsania • written by ahmad mousavi ahmad mousavi

2 answers

Try out CLUST https://baselabujamous.com/computational-methods/clust/ Read counts, FPKM, TPM etc are accepted in that.

Regarding genes, you might want to leave out genes which show zero or low variance across samples over and above the minimum cutoff for genes.

• 1 views • link
written by pixie@bioinfo pixie@bioinfo

Thanks,

Do you have any idea of creating such plot with R?

• 1 views • link
written by ahmad mousavi ahmad mousavi

you need to calculate the z-scores ((x-mean)/sigma) for the counts. From this, follow any tutorial on plotting in R. You should be able to get the above plots. With CLUST, you can get similar plots. https://github.com/BaselAbujamous/clust

• 1 views • link
written by pixie@bioinfo pixie@bioinfo

Hi, First you can find the top genes which are highly variable across the conditions. One way to do this is calculate the standard deviation (SD) across the sample for each gene and then arrange from high to low using SD values. You can choose top 10% or 20% of the genes and later parform k-means on them could give you nice patterns across the samples. The plot I have shown in the another post which is exactly done by this approach.

• 1 views • link
written by Chirag Parsania Chirag Parsania

Log in to answer this question.

More posts like this