[Bioc-sig-seq] readFastq

Martin Morgan mtmorgan at fhcrc.org
Wed Sep 1 03:55:23 CEST 2010


On 08/31/2010 06:18 PM, joseph wrote:
> can you import two lanes with the readFastq function to create a single 
> ShortReadQ object? 

Hi Joseph

after

  example(readFastq)

try

  fl = list.files(analysisPath(sp), "s_1_sequence.txt")
  readFastq(fl)
  readFastq(c(fl, fl))

also, something like

  readFastq(analysisPath(sp), "s_[123]_sequence.txt")

would read all files in the directory analysisPath(sp), matching the
pattern (see, e.g., ?list.files) "s_[123]_sequence.txt" into a single
object. Note that this requires enough memory to store all the data.
Also it might be interesting to compare the input speed of something like

  srapply(c(fl, fl), readFastq,
          reduce=function(lst) Reduce(append, lst))

before and after

  library(multicore)

and with larger data files than those provided as examples in the
ShortRead package.

Martin



> 
> Thanks
> Joseph
> 
> 
> 
>       
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioc-sig-sequencing mailing list
> Bioc-sig-sequencing at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing


-- 
Martin Morgan
Computational Biology / Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N.
PO Box 19024 Seattle, WA 98109

Location: Arnold Building M1 B861
Phone: (206) 667-2793



More information about the Bioc-sig-sequencing mailing list