Question: How to add FPKM values in cells of heatmap?
0
gravatar for Mehmet
14 months ago by
Mehmet460
Japan
Mehmet460 wrote:

Dear All,

I want to add FPKM values into cells of heatmap using complex heatmap package. I made a heatmap, I tried to use:

cell_fun = function(j, i, x, y, width, height, fill) {grid.text(sprintf("%.1f", df[i, j]), x, y, gp = gpar(fontsize = 10))},

df= my data matrix having FPKM values of each gene in each sample.

Some values are seen in cell, while others are different from df file.

rna-seq R • 576 views
ADD COMMENTlink modified 14 months ago by EagleEye6.3k • written 14 months ago by Mehmet460

I already figured out.

ADD REPLYlink written 14 months ago by Mehmet460
1

For the solution:

the data frame that has FPKM values of each gene must have only FPKM values of each gene in each sample/condition. Previously, I kept gene name in the data frame, but order of FPKM values in cells of heatmap was different that the data frame.

ADD REPLYlink written 14 months ago by Mehmet460

Great, can you paste the solution for everyone else? I, only now, just saw your request in the other question (I was traveling today)

ADD REPLYlink written 14 months ago by Kevin Blighe42k

Hi Kevin,

I would like to have your opinion about differential gene expression analysis workflow;

I have ~50 genes and have RNA reads obtained from four tissues. I want to see expression differences of those genes in the four tissues. I am planning to perform the steps after quality control: 1. Alignment of each RNA reads to reference genome separately. This will produce 4 .bam files. 2. Read counting of four .bam files. in this step, should I combine all bam files or should I do read counting for each .bam files separately? 3. Normalization. Which method do you recommend? 4. DEG analysis; edgeR or DESeq.

is my workflow suitable?

ADD REPLYlink written 12 months ago by Mehmet460

Hello Mehmet, which read counting method do you plan to use? For normalisation, my preference is DESeq2.

Your sample size is just 4?, that is, 1 sample per tissue?

ADD REPLYlink written 12 months ago by Kevin Blighe42k
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: 1215 users visited in the last hour