Limma and Quantarray files
2
0
Entering edit mode
@alessio-venier-1766
Last seen 9.6 years ago
Item Type: Appointment Begin time: Mon, 02 Oct 2006 05:00 pm +0000 Duration: 1 Hour Location: Hi, I've a problem using limma with quantarray output files. I've changed all output files to ANSI format, because the original Unicode format is not accepted by limma. I'm able to use one file per time without any trouble but when I try to load more than one file using a "target" file using, >target<-readTargets(file="targets.txt",path="C:\\Playground\\CrossTes t") >RG<-read.maimages(target,source="quantarray",path="C:\\Playground\\Cr ossTest") I receive an error message like: Read C:\Playground\CrossTest/PA52-Export.txt Error in if (skip > 0) readLines(file, skip) : argument is of length zero How can I solve this problem? Kind Regards Alessio
limma limma • 765 views
ADD COMMENT
0
Entering edit mode
Gavin Kelly ▴ 50
@gavin-kelly-1449
Last seen 9.6 years ago
> Hi, > > I've a problem using limma with quantarray output files. I've changed all output files to ANSI format, > because the original Unicode format is not accepted by limma. > I'm able to use one file per time without any trouble but when I try to load more than > one file using a "target" file using, > > >target<-readTargets(file="targets.txt",path="C:\\Playground\\CrossTes t" ) > > >RG<-read.maimages(target,source="quantarray",path="C:\\Playground\\Cr os sTest") > > I receive an error message like: > > Read C:\Playground\CrossTest/PA52-Export.txt > Error in if (skip > 0) readLines(file, skip) : > argument is of length zero > > How can I solve this problem? > > Kind Regards > > Alessio > Alessio, I had a similar problem in a recent version of limma (not sure which version you're using). If you edit the source for read.maimages, and look for the second instance of the line: skip <- grep("Begin Data", firstfield) and comment this out. Save this as your own function, and it should work okay. This will only work if all your quantarray files have the data (as flagged by the "Begin Data" string) starting at the same line. If this is not the case, then you'll need to find the instance of the phrase "nlines=100" just above the skip assignment, and remove it/customise it to your files (and uncomment the skip assignment) Hope this helps - Gavin -- Gavin Kelly (DPhil) Senior Statistician Bioinformatics & Biostatistics, Room 708 Cancer Research UK 44 Lincoln's Inn Fields London WC2A 3PX United Kingdom
ADD COMMENT
0
Entering edit mode
@alessio-venier-1766
Last seen 9.6 years ago
Dear Gordon and Gavin, many thanks for your prompt answer and support! Cheers Alessio -----Original Message----- From: "Gordon K Smyth" <smyth@wehi.edu.au> To: "Gavin Kelly" <gavin.kelly at="" cancer.org.uk="">, "Alessio Venier" <0509388V at student.gla.ac.uk> Date: Wed, 4 Oct 2006 21:56:03 +1000 (EST) Subject: [BioC] Limma and Quantarray files Dear Alessio and Gavin, Thanks for pointing out the problem with reading quantarray files. The bug has now been fixed in limma 2.9.1. limma now reads the annotation columns for quantarray files as well. Best wishes Gordon > Date: Wed, 4 Oct 2006 09:43:50 +0100 > From: "Gavin Kelly" <gavin.kelly at="" cancer.org.uk=""> > Subject: Re: [BioC] Limma and Quantarray files > To: <bioconductor at="" stat.math.ethz.ch=""> > Message-ID: > <1BFC6DAE48E6D4478BDFF41C012CD311FC4FB4 at uk-lif- lexc02.crwin.crnet.org> > Content-Type: text/plain; charset="us-ascii" > >> Hi, >> >> I've a problem using limma with quantarray output files. I've changed > all output files to ANSI format, >> because the original Unicode format is not accepted by limma. >> I'm able to use one file per time without any trouble but when I try > to load more than >> one file using a "target" file using, >> >> >>target<-readTargets(file="targets.txt",path="C:\\Playground\\CrossTe st" > ) >> >> >>RG<-read.maimages(target,source="quantarray",path="C:\\Playground\\C ros > sTest") >> >> I receive an error message like: >> >> Read C:\Playground\CrossTest/PA52-Export.txt >> Error in if (skip > 0) readLines(file, skip) : >> argument is of length zero >> >> How can I solve this problem? >> >> Kind Regards >> >> Alessio >> > > Alessio, > I had a similar problem in a recent version of limma (not sure which > version you're using). If you edit the source for read.maimages, and > look for the second instance of the line: > > skip <- grep("Begin Data", firstfield) > > and comment this out. Save this as your own function, and it should > work okay. This will only work if all your quantarray files have the > data (as flagged by the "Begin Data" string) starting at the same line. > If this is not the case, then you'll need to find the instance of the > phrase "nlines=100" just above the skip assignment, and remove > it/customise it to your files (and uncomment the skip assignment) > > Hope this helps - Gavin > > -- > Gavin Kelly (DPhil) Senior Statistician > Bioinformatics & Biostatistics, Room 708 > Cancer Research UK > 44 Lincoln's Inn Fields > London WC2A 3PX > United Kingdom ure which > version you're using). If you edit the source for read.maimages, and > look for the second instance of the line: > > skip <- grep("Begin Data", firstfield) > > and comment this out. Save this as your own function, and it should > work okay. This will only work if all your quantarray files have the > data (as flagged by the "Begin Data" string) starting at the same line. > If this is not the case, then you'll need to find the instance of the > phrase "nlines=100" just above the skip assignment, and remove > it/customise it to your files (and uncomment the skip ass
ADD COMMENT

Login before adding your answer.

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