Question: Which method is good for getting logFC for pairwise comparsion of replicates?
gravatar for mike
5.2 years ago by
mike60 wrote:


I have 4 induced and 4 control samples (technical replicates). I ran cuffdiff to find differentially expressed genes. Now I want to create Heat Map, showing differentially expressed genes across each replicates (not condition) using Log fold change values. Is there any method to extract FPKM values for each replicate and calculate Log fold change? Or would it be appropriate to use EdgeR to perform pairwise comparison (induced1 vs control1 , induced1 vs control2 ....)  between each replicate to get the log fold changes for the genes?

Note: I tried treating each replicate as a sample in cuffdiff but this method reduced the power of detecting differentially expressed genes because of the low coverage for replicates.


rna-seq • 1.9k views
ADD COMMENTlink modified 2.2 years ago by O.rka190 • written 5.2 years ago by mike60

Any progress on this? I'm running the metrics on a bunch of distances measures right now to see the best one for my data set.

ADD REPLYlink written 2.2 years ago by O.rka190

This thread is 3 years old. If you are using FPKM values and want to plot them in a heatmap, then convert them to the Z scale using zFPKM package, and then use Euclidean distance. Using Euclidean distance on the un-transformed FPKM values would be inappropriate; however, you could use something like 1 minus the correlation distance.

Se here: A: Heatmap based with FPKM values

ADD REPLYlink written 2.2 years ago by Kevin Blighe63k
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: 1687 users visited in the last hour