Question: Cufflinks output - how to produce heatmap of diff exp genes as z-scores
gravatar for tahneesaunders
2.6 years ago by
tahneesaunders0 wrote:

Hi all,

I have performed followed the Tophat alignment and Cufflinks differential analysis pipelines for my RNAseq data. I am not playing around in R with Cummerbund - I can produce pretty nice volcano plots and extract lists of significant genes.

I was wondering what the best way to produce a heatmap of z-scores using the cufflinks output is?

Appreciate all your suggestions.

Thanks. Tahnee.

(I have searched for about an hour to find an 'easy' solution, however cannot find one to my skill level, without being an R expert...).

zscore rna-seq cufflinks heatmap R • 1.5k views
ADD COMMENTlink modified 2.6 years ago by domenico.somma20 • written 2.6 years ago by tahneesaunders0
gravatar for Kevin Blighe
2.6 years ago by
Kevin Blighe63k
Kevin Blighe63k wrote:

Hi, you may want to take a look at my answer here, where I have done specifically this: Heatmap based with FPKM values

You will have to get your data in this format:

      Sample1 Sample2 Sample3 Sample4 ...
Gene1 45      453       2       4
Gene2 45       43     444      53
Gene3 40        3      23     232
Gene4 50       56      43     224
ADD COMMENTlink written 2.6 years ago by Kevin Blighe63k
gravatar for domenico.somma
2.6 years ago by
Glasgow University
domenico.somma20 wrote:

Hello, I wrote a python module named Papillon for this purpose. Python normally is easier to learn than R. You can find Papillon here. There is a tutorial, you can install it with pip or conda, and make a heatmap is very easy:

import papillon as pp

MyExp=pp.read_db("My Experiment")                # "My Experiment" is the folder name of your data

MyExp.get_gene()             # you can select gene/isoform o specific gene name (see the tutorial)


Best Regards Dom

ADD COMMENTlink written 2.6 years ago by domenico.somma20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 674 users visited in the last hour