warning (condition has length>1) when reading in Agilent hyb text files
1
0
Entering edit mode
@john-fernandes-1557
Last seen 10.2 years ago
I'm reading Agilent text files generated from Feature Extractor 9.5.3 into limma with the following: myfunFE <- function(x) { okAboveBG <- x[,"rIsWellAboveBG"]==1 & x[,"gIsWellAboveBG"]==1 okSaturated <- x[,"rIsSaturated"]==0 & x[,"gIsSaturated"]==0 okPopnOutlier <- x[,"rIsFeatPopnOL"]==0 & x[,"gIsFeatPopnOL"]==0 okNonUnifOutlier <- x[,"rIsFeatNonUnifOL"]==0 & x[,"gIsFeatNonUnifOL"]==0 as.numeric(okAboveBG & okSaturated & okPopnOutlier & okNonUnifOutlier) } RG <- read.maimages(targets$FileName,path=hybpath,source="agilent", columns=list(Rf="rMedianSignal",Gf="gMedianSignal",Rb="rBGMedianSi gnal",Gb="gBGMedianSignal"), other.columns=othercols,wt.fun=myfunFE) I am getting a warning for each file: "the condition has length > 1 and only the first element will be used in: if (text.to.search != "") for (i in 1:ncn)" However, I am getting results from the weight function. Does the warning indicate a problem? I am running R 2.4.1 on Win XP with limma 2.9.17. Thanks, John Fernandes
limma limma • 1.5k views
ADD COMMENT
0
Entering edit mode
@john-fernandes-1557
Last seen 10.2 years ago
Anyone have an idea what could be causing the warning I'm getting? Do you need more info from me? Can I look up the error somewhere? Thanks, John > I'm reading Agilent text files generated from Feature Extractor 9.5.3 into > limma with the following: > > myfunFE <- function(x) { > okAboveBG <- x[,"rIsWellAboveBG"]==1 & x[,"gIsWellAboveBG"]==1 > okSaturated <- x[,"rIsSaturated"]==0 & x[,"gIsSaturated"]==0 > okPopnOutlier <- x[,"rIsFeatPopnOL"]==0 & x[,"gIsFeatPopnOL"]==0 > okNonUnifOutlier <- x[,"rIsFeatNonUnifOL"]==0 & > x[,"gIsFeatNonUnifOL"]==0 > as.numeric(okAboveBG & okSaturated & okPopnOutlier & > okNonUnifOutlier) > } > > RG <- read.maimages(targets$FileName,path=hybpath,source="agilent", > > columns=list(Rf="rMedianSignal",Gf="gMedianSignal",Rb="rBGMedianSign al",Gb="gBGMedianSignal"), > other.columns=othercols,wt.fun=myfunFE) > > I am getting a warning for each file: > > "the condition has length > 1 and only the first element will be used > in: if (text.to.search != "") for (i in 1:ncn)" > > > However, I am getting results from the weight function. Does the warning > indicate a problem? > > I am running R 2.4.1 on Win XP with limma 2.9.17. > > Thanks, > John Fernandes > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor >
ADD COMMENT
0
Entering edit mode
Hi John, limma 2.9.17 is old and this problem seems to have been fixed in more recent versions of limma: hpages at lamb1:~/svn/Rpacks-2.1/limma/R> svn log -v -r 24408 -------------------------------------------------------------------- ---- r24408 | smyth at wehi.edu.au | 2007-04-30 21:17:31 -0700 (Mon, 30 Apr 2007) | 5 lines Changed paths: M /trunk/madman/Rpacks/limma/DESCRIPTION M /trunk/madman/Rpacks/limma/R/read-maimages.R M /trunk/madman/Rpacks/limma/inst/doc/changelog.txt 1 May 2007: limma 2.11.1 - Bug fix to read.columns() to stop spurious warning message when text.to.search has length greater than one. -------------------------------------------------------------------- ---- I recommend you update your limma package. Also you could update to R 2.5.0 and Bioconductor 2.0: many things have changed between BioC 1.9 and BioC 2.0 and a lot of bugs have been fixed ;-) Cheers, H. John Fernandes wrote: > Anyone have an idea what could be causing the warning I'm getting? Do > you need more info from me? Can I look up the error somewhere? > > Thanks, > John > > >> I'm reading Agilent text files generated from Feature Extractor 9.5.3 into >> limma with the following: >> >> myfunFE <- function(x) { >> okAboveBG <- x[,"rIsWellAboveBG"]==1 & x[,"gIsWellAboveBG"]==1 >> okSaturated <- x[,"rIsSaturated"]==0 & x[,"gIsSaturated"]==0 >> okPopnOutlier <- x[,"rIsFeatPopnOL"]==0 & x[,"gIsFeatPopnOL"]==0 >> okNonUnifOutlier <- x[,"rIsFeatNonUnifOL"]==0 & >> x[,"gIsFeatNonUnifOL"]==0 >> as.numeric(okAboveBG & okSaturated & okPopnOutlier & >> okNonUnifOutlier) >> } >> >> RG <- read.maimages(targets$FileName,path=hybpath,source="agilent", >> >> columns=list(Rf="rMedianSignal",Gf="gMedianSignal",Rb="rBGMedianSig nal",Gb="gBGMedianSignal"), >> other.columns=othercols,wt.fun=myfunFE) >> >> I am getting a warning for each file: >> >> "the condition has length > 1 and only the first element will be used >> in: if (text.to.search != "") for (i in 1:ncn)" >> >> >> However, I am getting results from the weight function. Does the warning >> indicate a problem? >> >> I am running R 2.4.1 on Win XP with limma 2.9.17. >> >> Thanks, >> John Fernandes >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at stat.math.ethz.ch >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD REPLY

Login before adding your answer.

Traffic: 712 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