[BioC] help with fit function

Alberto Goldoni alberto.goldoni1975 at gmail.com
Wed Sep 8 15:49:12 CEST 2010


Dear All,
thanks to the function: fit <- eBayes(fit) i can obtain the toptable like this:

> topTable(fit,coef=2)
      Row Col ProbeUID ControlType    ProbeName  GeneName SystematicName
5118   26  12     4558           1    (+)eQC-13       eQC            eQC
4520   23   1     4296           1    (+)eQC-02       eQC            eQC
6868   34  73     6435           0 A_42_P490236 409568_Rn      409568_Rn
4932   25  19     4381           1    (+)eQC-20       eQC            eQC
4729   24   2     4496           1    (+)eQC-03       eQC            eQC
4464   22 158     4240           0  A_43_P11468 NM_012571      NM_012571
4735   24  10     4502           1    (+)eQC-11       eQC            eQC
18777  92  42    17921           0 A_42_P781006 296921_Rn      296921_Rn
5110   26   4     4551           1    (+)eQC-05       eQC            eQC
4741   24  18     4508           1    (+)eQC-19       eQC            eQC

but from the "fit" table i would like to obtain a new table without
row in witch the value "fit$ControlType" would be "+1" and "-1".
like this:

      Row Col ProbeUID ControlType    ProbeName  GeneName SystematicName
6868   34  73     6435           0 A_42_P490236 409568_Rn      409568_Rn
4464   22 158     4240           0  A_43_P11468 NM_012571      NM_012571
18777  92  42    17921           0 A_42_P781006 296921_Rn      296921_Rn

is it possible, somebody can help me?

Best regards.





-- 
-----------------------------------------------------
Dr. Alberto Goldoni
Parma, Italy



More information about the Bioconductor mailing list