[BioC] Reading Imagene files with read.imagene in limma

[Gordon Smyth]

To exclude any set of flag values, define

mywtfun <- function(exclude.flags=c(1,2,3)) function(obj) 1-(obj$Flag %in% 
exclude.flags)

Now you can use

RG <- read.maimages(files, source="imagene", wt.fun=mywtfun(c(1,2)))

to remove only 1 and 2, or

RG <- read.maimages(files, source="imagene", wt.fun=mywtfun(c(1,2,4)))

to remove 1, 2 and 4, etc, etc.

Gordon

At 06:04 PM 8/06/2004, STKH (Steen Krogsgaard) wrote:
>Hi Gordon,
>
>I really don't want to reinvent the wheel (or any other part of the
>excellent functionality of limma). What I want to do is to save the flag
>columns from the data files in the RGList dataframe as RGList$flags and
>then set the weights to either 0 or 1 depending on the values of the
>flags. Since the flags say that the spot is either OK or unusable (with
>nothing in between) it doesn't really make sense (to me at least) to use
>weights between 0 and 1.
>
>I couldn't figure out how to add the flags "column" to the RGList within
>the wt.fun function, so that's why I wrapped read.imagene instead. The
>pseudocode is like this:
>read.my.imagene=function(bla bla bla, exclude.flags=c(1,2,3))
>{
>         obj=read.imagene(bla bla bla, wt.fun=wt.myimagene)
>         obj$flags=obj$weights
>         obj$weights=rep(1,length(obj$weights))
>         obj$weigths[obj$flags %in% exclude.flags]=0
>         obj
>}
>
>wt.myimagene=function(ima)
>{
>         ima[,"Flags"]
>}
>
>so, during read.imagene the Flags column from the datafiles are put in
>the weigths list (it's a list, right?), and this is subsequently copied
>to a $flags list. Then $weights are set to all 1's, and the $weights
>where $flags is in an exclude-list (exclude.flags) is set to 0. This
>mimicks the data preparation step used in GeneSight from BioDiscovery
>where "bad" flags are filterede out. The reason for not putting the
>filtering into wt.myimagene is that I want to let the user specify which
>flags to exclude. To exclude only flags 1 and 2, for instance, the call
>would be
>RG=read.my.imagene(files=files, exclude.flags=c(1,2))
>
>This is probably not the best way to accomplish this, but it's the one I
>could think up! I can see that for instance wtarea is defined by two
>function-statements and that this allow for a parameter to be specified
>in the call to read.imagene, but I don't really understand this
>notation. Besides, I still would like to have the original flag values
>stored in the RGList object so that I can re-do the filtering later or
>plot the flagges spots or whatever.
>
>Thanks for all the feedback
>Steen
>
>-----Original Message-----
>From: Gordon Smyth [mailto:smyth at wehi.edu.au]
>Sent: 8. juni 2004 00:38
>To: STKH (Steen Krogsgaard)
>Cc: bioconductor at stat.math.ethz.ch
>Subject: Re: [BioC] Reading Imagene files with read.imagene in limma
>
>
>At 08:14 PM 7/06/2004, STKH (Steen Krogsgaard) wrote:
> >Hi,
> >
> >I have a bunch of Imagene data files that I want to read into limma
> >using read.imagene. It all works very nicely like this:
> >
> >RG<-read.imagene(files), where files is a matrix with two columns and
> >one row for each array. The names of the two dye-files are in the two
> >columns.
> >
> >I want to wrap this call into another function in order to do some
> >processing of the quality flags that Imagene uses (I want to change the
>
> >flags into weights). The idea is:
>
>Why re-invent the wheel? read.maimages() already gives you complete
>facility to do this. The argument wt.fun given to read.maimages() can be
>
>any function you like which takes an imagene file as a data frame and
>converts flags into weights. The help says:
>
>'wt.fun' may be any user-supplied function which accepts
>       a data.frame argument and returns a vector of non-negative
>       weights. The columns of the data.frame are as in the image
>       analysis output files. See 'QualityWeights' for provided weight
>       functions.
>
>Gordon
>
> >read.my.imagene <-function(files, path = NULL, ext = NULL, names =
> >NULL, columns = NULL,
> >     wt.fun = NULL, verbose = TRUE, sep = "      ", quote = "\"",
> >exclude.flags=c(1,2,3), ...)
> >{
> >         obj=read.imagene(files=files, path=path, ext=ext, names=names,
>
> >columns=columns, wt.fun=wt.fun, verbose=verbose, sep=sep, quote=quote,
> >...)
> >         // do some more stuff
> >}
> >
> >but when I call read.my.imagene I get an error "Error in scan(file =
> >file, what = what, sep = sep, quote = quote, dec = dec,  :
> >         line 1 did not have 27 elements". I understand that scan is
> >complaining about the number of columns in the data file, but I'm
> >really puzzled why it works when I call read.imagene directly but not
> >when I call read.my.imagene which then in turn calls read.imagene.
> >
> >I also tried to copy the source code for read.imagene into a new
> >function (which I called test.read), and this produced the exact same
> >error. There were no code changes done, just a copy of the source code
> >into Notepad (yes, I'm on a Windows box) preceeded by "test.read = ".
> >
> >I use R 1.9.0 and the latest Bioconductor (downloaded last week).
> >
> >thanks in advance for any suggestions
> >Steen