[BioC] Limma: bad spots flagged out?

Ana Conesa aconesa at ivia.es
Mon Feb 20 12:58:15 CET 2006


   Dear list,
   I   have   a  doubt  about  the  real  use  of  spots  weights  during
   normalization in limma. According with the documentation spots weights
   equal  to  0  are  ignored  during  normalization (and other posterior
   analyses)  , and although these spots are not removed they do not have
   any  influence on the rest. I had observed some strange behavior on my
   normalized  data  and made a tried to make a check on this. What I did
   was  to  replace  weights==0 spots by NA and redo analysis. I have
   found  the  results  do  no   spots  are  adequately ignored by limma or I made a conceptual mistake
   in this check. This is exactly the code I used for checking:
   > RG.b <- backgroundCorrect(RG)
   > MA.1 <- normalizeWithinArrays(RG.b)
   > RG.b$G[RG.b$weights==0] <- NA
   > RG.b$R[RG.b$weights==0] <- NA
   > MA.2 <- normalizeWithinArrays(RG.b)
   > MA.1$M[2,]
            BFN33        S50.05        S65.02        S65.03        S71.02
   S75.01      S77.03     control     control     control     control
   -0.96794887     0.06715693    -0.08766477    -0.50161127   -1.25216169
   -0.80724650    -0.61351625    -0.80751427    -0.49960303   -0.66912129
   0.27447918
   > MA.2$M[2,]
            BFN33        S50.05        S65.02        S65.03        S71.02
   S75.01      S77.03     control     control     control     control
   -0.94108221     0.04803980    -0.08030398    -0.49226094            NA
   -0.79588160             NA    -0.77182854             NA   -0.61426784
   0.18076174
   Thank you
   Ana

      O @@@@@     Ana Conesa, PhD.
     @@@ O @@ O @   Centro de Genómica
      @  O  @@@@  O  @   Instituto Valenciano de Investigaciones Agrarias
   (IVIA)
     @@@ O @@@@        @@@@ O @     46113 Moncada (Valencia) SPAIN
        ||       Tel. +34 963424000 ext.70161; Fax. +34 963424001
        ||       email: aconesa at ivia.es


More information about the Bioconductor mailing list