[BioC] Plate or print-order effects - limma

Jenny Drnevich drnevich at uiuc.edu
Tue Jan 2 21:38:52 CET 2007


Hi Gordon,

FYI - I was playing around with plotting print order effects, and found a 
discrepancy between the help documentation and the actual behavior of 
'normalizeForPrintorder.rg '. The help file states "The function 
plotPrintorder or normalizeForPrintorder.rg with plot=TRUE returns no value 
but produces a plot as a side-effect." However, when I used 
'normalizeForPrintorder.rg' with plot=TRUE, it not only produced the plot 
but spit out the entire list!

I'm using this function instead of the wrapper 'plotPrintorder' because I'm 
having trouble making the necessary layout list. I could only find a little 
bit of information about the layout list in 'printorder' help page - is 
each component of the list (ngrid.r, ngrid.c, nspot.r and nspot.c) just a 
single number? In my case our printer has 48 tips and prints 12 x 4 blocks 
of 420 spots, one dip per spot, so would my layout just be a list of 
ngrid.r = 12, ngrid.c = 4, nspot.r = 420, and nspot.c = 1?

Thanks,
Jenny

 > sessionInfo()
R version 2.4.0 (2006-10-03)
i386-pc-mingw32

locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United 
States.1252;LC_MONETARY=English_United 
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

attached base packages:
[1] "splines"   "tools"     "methods"   "stats"     "graphics"  "grDevices"
[7] "utils"     "datasets"  "base"

other attached packages:
  affyQCReport    simpleaffy         made4 scatterplot3d          ade4
      "1.12.0"       "2.8.0"       "1.8.0"      "0.3-24"       "1.4-2"
       affyPLM      affydata affycoretools       annaffy        xtable
      "1.10.0"      "1.10.0"       "1.7.5"       "1.6.0"       "1.4-2"
         gcrma   matchprobes       biomaRt         RCurl           XML
       "2.6.0"       "1.6.0"       "1.8.0"       "0.7-0"       "1.2-0"
       GOstats      Category    genefilter      survival          KEGG
       "2.0.3"       "2.0.3"      "1.12.0"        "2.29"      "1.14.1"
          RBGL      annotate            GO         graph         limma
      "1.10.0"      "1.12.0"      "1.14.1"      "1.12.0"       "2.9.1"
          affy        affyio       Biobase       RWinEdt
      "1.12.1"       "1.2.0"      "1.12.2"       "1.7-5"


At 09:32 PM 12/31/2006, Gordon Smyth wrote:
>Dear Joao,
>
>See ?plotPrintorder and
>http://www.statsci.org/smyth/pubs/porder/porder.html. (I'm a bit
>puzzled why you haven't already tried this, considering that you're
>already using limma and are aware of printorder effects.)
>
>Note that GAL file order is not the same as printorder in general, so
>you can't just use plot().
>
>Best wishes
>Gordon
>
> >Date: Fri, 22 Dec 2006 19:33:11 +0000
> >From: J.delasHeras at ed.ac.uk
> >Subject: Re: [BioC] Plate or print-order effects - limma
> >To: bioconductor at stat.math.ethz.ch
> >
> >Just use 'plot':
> >
> >is MA is teh name of your MAlist containing your M and A values, then
> >
> >plot(MA$A[,1],MA$M[,1]) #prints the usual MA plot for slide 1
> >
> >plot(MA$M[,1]) #prints the M values, in the order given in the gal
> >file (by block)... is that what you wanted?
> >
> >just add your favourite parameters to 'plot' to make it look the way
> >you want it.
> >
> >Jose
> >
> >
> >Quoting Jo?o Fadista <Joao.Fadista at agrsci.dk>:
> >
> > > Dear all,
> > >
> > > I would like to visualize, for each array, the plate or print-order
> > > effects. To do that I want to be able to plot the M values in the
> > > vertical axis and the print order (i.e., the numerical order in
> > > which the spots were laid down during the printing of the array) in
> > > the horizontal axis. I am using genepix files.
> > >
> > > Thanks in advance,
> > >
> > > Best regards
> > >
> > > Jo?o Fadista
> > > Ph.d. student
> > >
> > >
> > >
> > >       Danish Institute of Agricultural Sciences
> > > Research Centre Foulum
> > > Dept. of Genetics and Biotechnology
> > > Blichers All? 20, P.O. BOX 50
> > > DK-8830 Tjele
> > >
> > > Phone:        +45 8999 1900
> > > Direct:       +45 8999 8999
> > > E-mail:       Joao.Fadista at agrsci.dk <mailto:Joao.Fadista at agrsci.dk>
> > > Web:  www.agrsci.org <http://www.agrsci.org/>
> > > ________________________________
> > >
> > > News and news media 
> <http://www.agrsci.org/navigation/nyheder_og_presse> .
> > >
> > > This email may contain information that is confidential. Any use or
> > > publication of this email without written permission from DIAS is
> > > not allowed. If you are not the intended recipient, please notify
> > > DIAS immediately and delete this email.
> > >
> > >
> > >       [[alternative HTML version deleted]]
> > >
> > >
> >
> >
> >
> >--
> >Dr. Jose I. de las Heras                      Email: J.delasHeras at ed.ac.uk
> >The Wellcome Trust Centre for Cell Biology    Phone: +44 (0)131 6513374
> >Institute for Cell & Molecular Biology        Fax:   +44 (0)131 6507360
> >Swann Building, Mayfield Road
> >University of Edinburgh
> >Edinburgh EH9 3JR
> >UK
>
>_______________________________________________
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch
>https://stat.ethz.ch/mailman/listinfo/bioconductor
>Search the archives: 
>http://news.gmane.org/gmane.science.biology.informatics.conductor

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

330 ERML
1201 W. Gregory Dr.
Urbana, IL 61801
USA

ph: 217-244-7355
fax: 217-265-5066
e-mail: drnevich at uiuc.edu



More information about the Bioconductor mailing list