[BioC] Re: read.GenePix() function

[Gordon Smyth]

This can be done by providing a user-specified function to the wt.fun 
argument of read.maimages() in the limma package. Just use something like

mywtfun <- function(dat) as.numeric( (dat$RF>dat$RB+2*dat$RBSD) & 
(dat$GF>dat$GB+2*dat$GBSD))
RG <- read.maimages(files, 
columns=list(Gf="GF,Rf="RF",Gb="GB",Rb="RB"),wt.fun=mywtfun)

Gordon

> > I want to specify these two columns is because I want to write a
> > function to filter spots whose RF and GF are less than RB+2*RBSD and
> > GB+2*GBSD. Do you know an approach I can do that?