limma question
1
0
Entering edit mode
@straubhaar-juerg-391
Last seen 9.7 years ago
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]]
limma limma • 633 views
ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 2 minutes ago
WEHI, Melbourne, Australia
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
ADD COMMENT

Login before adding your answer.

Traffic: 562 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6