Two suggestions - 1) The code I gave you before was written as if your
reference was in the green (Cy3) channel. However, based on the
your 'modelMatrix(targets, ref="gDNA")' command, your reference is in
red (Cy5) channel. Therefore, you would have to reverse some of the
commands where appropriate (e.g., use method "Rquantile" and replace M
values with G values).
2) Find a local statistician to consult about the analysis, because it
appears you have a 2 x 7 factorial design (2 strains and 7 timepoints
treatment groups total). There are variety of ways to analyze this
experimental design, depending on what all you want to know from the
If you really only wanted to know which genes were different between
wt at each time point independently, then you could analyze the
from each time point separately. However, there is so much more
to be gained from this data set, which is why I suggest you consult a
Best of luck,
At 10:54 AM 11/21/2006, yanju wrote:
>I adapted my code according to your suggestion. Then at some time
>the results showed that the most differently expressed genes are
>This is every werld.
>And It doesnt matter if I change -1 in the design matrix to 1 (my
>new design matrix=old design matrix* -1, old design matrix was
>from modelMatrix function) or not. I mean this didnt effect my
>Since I could not figure it out, I paste my code here. Hope you could
>me what's wrong with my program. Basic information of the data:
>Two-channle array, 7 time points from 16-72h, at each time point
>some repelicants. Aim is to detect the different expressed genes at
> From the very begin of the code:
> # read targets and genepix files
> # bacground correction
> #design was similar like follows.
> # wt16 wt20 wt24
> #[1,] -1 0 0
> #[2,] 0 -1 0
> #[3,] 0 0 -1
> #Then it was multiply by -1 to have the positive value.
> MUvsWT24=mu24-wt24, MUvsWT36=mu36-wt36,
> MUvsWT60=mu60-wt60, MUvsWT72=mu72-wt72,levels=design_revise)
> #fit the data to linear model and Bayes statistical summary
> #detecting the top20 differently expressed genes at time
> #But as I said, most of the top20 were markers or the "no
>Hope you could help me figure out the problems. I really appreciate
Jenny Drnevich, Ph.D.
Functional Genomics Bioinformatics Specialist
W.M. Keck Center for Comparative and Functional Genomics
Roy J. Carver Biotechnology Center
University of Illinois, Urbana-Champaign
1201 W. Gregory Dr.
Urbana, IL 61801
e-mail: drnevich at uiuc.edu