[BioC] how to design matrix by edgeR

wang peter wng.peter at gmail.com
Tue May 22 20:48:55 CEST 2012

dear Dr.Gordon:
   thank  you for your patient reply, but i still donot understand
what problem in my factor design?i have 3 questions as below:

the 1st question:
   my experiments involves two factors, time and treatment.
   so i define:
 time=factor( c('0h','0h','0h','6h','6h','6h','6h','12h',

the 2nd question:
i want to get DE cross time considersing control.
and my model.matrix is
design <- model.matrix(~time+treatment)?
or design <- model.matrix(~treatment+time)?

and yours is design <- model.matrix(~time+time:treatment)
which one is correct?

the third question:
      i donot do what is the difference between test by glmLRT
and test by pairedwised test functions?

To test for treatment effect at time 0:

  lrt0 <- glmLRT(y,fit,coef=8)

To test for treatment effect at time 1:

  lrt1 <- glmLRT(y,fit,coef=9)

thank you very much
if possible, can you give me some detailed material
of glmLRT or glmFIT
now i can only read your paper published recently

shan gao

More information about the Bioconductor mailing list