[BioC] EdgeR: Log FC does not changes when I swap the groups

saad khan [guest] guest at bioconductor.org
Mon Oct 10 21:41:30 CEST 2011


I am running edgeR on RNA-seq data with two groups Control(C1,C2,C3) and Treated(P1,P2,P3) having three replicates each. I need the log FC ratio with Treated(P1,P2,P3) value in the numerator and Control(C1,C2,C3) in the denominator but the results I am getting with EdgeR are not well correlated with the read counts data that I have.

e.g.: 
for gene Glyma12g14980.1
 fold change (lg Fc) = 33.21571049
Read Count: Control 1= 20, Control2= 136, COntrol 3= 15, P1= 0, P2=  0, P3=  0

Glyma10g03970.1
Fold Change= -37.80306427
Read Count: Control1=0, Control2=0, Control 3=0, P1=  815, P2=1590, P3=1649

even when I swap the groups the result remain the same i.e.:  group<- c(rep(\"Control\",3),rep(\"Treated\",3)) and  group<- c(rep(\"Treated\",3),rep(\"Control\",3)) the results are the same.

Please suggest where I am going wrong.

Regards
Saad Khan
Graduate Student
University of Missouri

 -- output of sessionInfo(): 

***

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list