Question: Microarray Analysis Using Dchip-Problem After Log Transforming
0
gravatar for pixie@bioinfo
6.1 years ago by
pixie@bioinfo1.4k
pixie@bioinfo1.4k wrote:

I am using dChip for analyzing Affy Cel files. I want to find differentially expressed genes in my data. I load the cel files, normalize and log transform it. However, I am not able to get any differentially expressed genes in my samples (using 2 fold change), even though there are around 57k probe ids and the paper corresponding to the dataset reports around 5K DEGs with 5 fold change. If I do not log-transform it, I get a good number of DEGs. I am not able to understand what is happening ..kindly help!

microarray • 1.8k views
ADD COMMENTlink modified 6.1 years ago by Michael Dondrup46k • written 6.1 years ago by pixie@bioinfo1.4k
0
gravatar for Michael Dondrup
6.1 years ago by
Bergen, Norway
Michael Dondrup46k wrote:

The short answer: something is different between your methods and their methods. If you want to replicate the outcome of an analysis you need to do exactly the same as the paper you want to replicate. You have to find out, what the differences are and eliminate them.

The rest is too much guesswork for us, because there are many things that can go wrong and you neither told us what you did exactly nor which paper you are referring to. That said (and no offence), it is quite possible that you have made a simple mistake in your calculus with the logarithm.

When you want to relate fold-change cutoffs on logged values to linear fold changes you have to calculate and filter for log differences (-) not log fractions (/). For example if you want to look for genes that have two-fold change on linear scale you have to look for all genes with absolute log2-difference ≥ 1 or for an absolute log2 difference of ~2.321928, because

log(x/y) = log(x) - log(y) != log(x)/log(y)

That would explain the discrepancy.

ADD COMMENTlink modified 6.1 years ago • written 6.1 years ago by Michael Dondrup46k

yes you are right...I had taken log2 transformed it while i was applying a fold change for log10 base. Thanks!

ADD REPLYlink written 6.1 years ago by pixie@bioinfo1.4k
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: 908 users visited in the last hour