Question: Cufflinks output - how to produce heatmap of diff exp genes as z-scores
0
gravatar for tahneesaunders
22 months 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.1k views
ADD COMMENTlink modified 21 months ago by domenico.somma20 • written 22 months ago by tahneesaunders0
0
gravatar for Kevin Blighe
22 months ago by
Kevin Blighe49k
Kevin Blighe49k 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 22 months ago by Kevin Blighe49k
0
gravatar for domenico.somma
21 months 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)

MyExp.heatmap(z_score=True)

Best Regards Dom

ADD COMMENTlink written 21 months ago by domenico.somma20
Please log in to add an answer.

Help
Access

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