[BioC] Spatial Correction

Rafael A. Irizarry ririzarr at jhsph.edu
Mon Jun 9 02:58:00 MEST 2003


all robust summaries ive used: li and wong's algo, median polish and rlm
within RMA appear to find these spatial artifacts as outliers (under a
typical definition of an outlier). an easy way to see this is to make an
image plot of the absolute value of the residuals from the model fit.



On Sun, 8 Jun 2003, Justin Borevitz wrote:

> Hi all great to be part of a highly active/responsive group!
> 
> Has anyone else thought about spatial correction for affy arrays?  Charles
> Berry and I have done some simple corrections using sliding windows, see
> Figure 1. http://www.genome.org/cgi/content/full/13/3/513/F1 or
> http://naturalvariation.org/sfp/Borevitz_etal_2003_GR.pdf  I know dCHIP can
> mask out scratches.  Could there be something that would take advantage of
> multichip analysis, kind of like quantile normalization but that would
> consider the data in its original image form (rectangular matrix of feature
> intensities)?  I often see smudges on arrays where useful data is below that
> needs to be locally corrected for.
> 
> FYI I have the cdf package for the Arabidopsis ATH1 array here.
> 
> http://natural.salk.edu/ath1121501cdf.tar.gz   ## source
> http://natural.salk.edu/ath1121501cdf.zip   ## win32
> 
> original cdf from affy
> http://natural.salk.edu/Ath1.cdf.zip
> 
> Justin
> http://naturalvariation.org 
> 
> -----Original Message-----
> From: bioconductor-bounces at stat.math.ethz.ch
> [mailto:bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Johannes Hüsing
> Sent: Saturday, June 07, 2003 1:47 AM
> To: Tapan Mehta
> Cc: bioconductor at stat.math.ethz.ch; John Zhang
> Subject: Re: [BioC] Bioconductor commands in Unix online documentation
> request
> 
> Tapan Mehta <tapmehta at yahoo.com> [Sat, Jun 07, 2003 at 09:07:27AM CEST]:
> > Thank you all for the help. It would be nice if I
> > could know about some documentation for R/bioconductor
> > commands that can be used in Unix through command
> > line.
> 
> There is no direct interface to internal R commands from any shell I
> know, if you mean that.
> 
> You can always type 
> 
> echo "sqrt(2)" | R CMD BATCH
> 
> and then inspect .R.Rout, which is the closest thing to what I think you
> mean.
> 
> Greetings
> 
> 
> Johannes
> -- 
> Johannes Hüsing   There is something fascinating about science. One gets
> hannes at ruhrau.de  such wholesale returns of conjecture from such a 
>                   trifling investment of fact.                Mark Twain
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>



More information about the Bioconductor mailing list