I'm reading smd file into limma and would like to give weight of zero
to
spots with FLAG != 0 (-50 in my file). I type the following:
* RG<-read.maimages("/home/jstraubh/mello/data/cex211.xls",
source="smd", wt.fun=wtflags(0), fill=T)
After entering this generates an error:
Error in "[<-"(*tmp*, , i, value = numeric(0)) :
nothing to replace with
Thank you for your help.
Juerg Straubhaar,
University of Massachusetts Medical School
[[alternative HTML version deleted]]
At 12:23 AM 13/09/2003, Straubhaar, Juerg wrote:
>I'm reading smd file into limma and would like to give weight of zero
to
>spots with FLAG != 0 (-50 in my file). I type the following:
>
>* RG<-read.maimages("/home/jstraubh/mello/data/cex211.xls",
>source="smd", wt.fun=wtflags(0), fill=T)
1. Have you checked that this command works without use of the
'wt.fun'
argument? The fact that your data file has an 'xls' extension suggests
that
it is an excel spreadsheet. If this is the case, you need to use excel
to
save your file as tab-delimited text. Like most read commands in R,
read.maimages() is designed to read text files, not proprietry binary
formats.
2. The function wtflags() is only for Genepix data, see ?wt.flags. If
you
have a column of your data file called "FLAG", then you need
wt.fun = function(x) x$FLAG >= 0
For example,
RG <- read.maimages(filenames, source="smd", wt.fun= function(x)
x$FLAG >= 0)
Gordon
>After entering this generates an error:
>
>Error in "[<-"(*tmp*, , i, value = numeric(0)) :
> nothing to replace with
>
>Thank you for your help.
>
>Juerg Straubhaar,
>University of Massachusetts Medical School