[BioC] how to merge replicate spots?

Barbara Cegielska barceg at ibch.poznan.pl
Mon Jul 6 12:53:08 CEST 2009


Michael,

No, it wasn't deliberate. I've change it and it seems  to work. Is it 
correct?

################################################################

 > fit <- lmFit(a, design)
 > contrast.matrix<-makeContrasts(AML.M2-K,levels=design)
 > fit2<-contrasts.fit(fit,contrast.matrix)
 > fit2<-eBayes(fit2)
 > topTable(fit2,adjust="BH")
    Block Row Column    ID    Name Status     logFC  AveExpr         
t      P.Value  adj.P.Val        B
850    35   6      1  6767    ST13 Oligos  3.075192 7.234434  9.388662 
7.710111e-05 0.02516929 2.306532
284    11   5      4  6280  S100A9 Oligos -3.005805 8.759735 -9.303701 
8.122192e-05 0.02516929 2.257757
859    35   9      1 23166   STAB1 Oligos  2.929309 7.131428  9.017997 
9.706814e-05 0.02516929 2.089691
841    35   3      1  6726    SRP9 Oligos  2.847275 7.485457  8.970909 
1.000094e-04 0.02516929 2.061377
557    21   6      4  3717    JAK2 Oligos  2.315902 8.071873  8.413855 
1.439105e-04 0.02516929 1.712576
854    35   7      4  4772  NFATC1 Oligos  2.453980 7.217102  8.247170 
1.611193e-04 0.02516929 1.603018
785    33   2      4  6093   ROCK1 Oligos  2.256970 7.425680  8.150150 
1.722239e-04 0.02516929 1.538097
436    17   2      1  3563   IL3RA Oligos  2.554791 7.891780  7.911541 
2.034941e-04 0.02530803 1.374718
937    37  10     10   291 SLC25A4 Oligos  2.507008 8.085160  7.737719 
2.304150e-04 0.02530803 1.252276
972    37  17     10   369    ARAF Oligos  2.282989 6.985434  7.484587 
2.772916e-04 0.02530803 1.068606
 > dim(fit2)
[1] 1023    1

####################################################################
P.s Sorry for not keeping discussions "on the list".

Regards,
B.



More information about the Bioconductor mailing list