[BioC] differntially expressed gene
Prasad Siddavatam
siddavatam at gmail.com
Wed Mar 9 17:56:04 CET 2011
Hi,
I sincerely thank Dr.Brown for his response to my earlier post on design matrix.
I have two questions....1. about the discrepancies in my DE results
2. B-values
These are my target files
FileName Cy3 Cy5 Original Names
HIDEN_1.gpr Ref HI_Inf Heat_Inactivated_1
HIDEN_2.gpr Ref HI_Inf Heat_Inactivated_2
HIDEN_3.gpr Ref HI_Inf Heat_Inactivated_3
infected_1.gpr Ref Infect Live_Infection_1
infected_2.gpr Ref Infect Live_Infection_2
infected_3.gpr Ref Infect Live_Infection_3
design:
HI_Inf Infect
1 0
1 0
1 0
0 1
0 1
0 1
contrast:
Contrasts
Levels HI_INF INF INFvsHI_INF
HI_Inf 1 0 -1
Infect 0 1 1
When I used the above matrix and contrasts, I found 270 and 2484 DE genes (for
HI_INF and INF, respectively).
But when I divided the data into two separate analyses
1. For HI_INF I found 608 DE genes
HIDEN_1.gpr Ref HI_Inf Heat_Inactivated_1
HIDEN_2.gpr Ref HI_Inf Heat_Inactivated_2
HIDEN_3.gpr Ref HI_Inf Heat_Inactivated_3
2. For INF I found 868 DE genes
infected_1.gpr Ref Infect Live_Infection_1
infected_2.gpr Ref Infect Live_Infection_2
infected_3.gpr Ref Infect Live_Infection_3
Why is this difference? technically those should be same because rest of the
steps were similar between the two.
------------------------------------------------------------------------
I also found that there are some genes with negative "B Values" but < 0.5
adjusted p.values and p.values...see below
logFC t P.Value adj.P.Val B
-0.6740520 -3.655211 0.006576695 0.04978582 -2.453619
-0.3866386 -3.655013 0.006578564 0.04978582 -2.453912
-0.6554844 -3.652845 0.006599049 0.04992410 -2.457116
In this case, can I delete the genes with negative B-values though adjusted
p.values and p.values are < 0.05?
Your suggestions are highly appreciated
-regards
Prasad
More information about the Bioconductor
mailing list