[BioC] how to deal with a 30G fastq file
Steve Lianoglou
mailinglist.honeypot at gmail.com
Thu Oct 6 15:48:23 CEST 2011
Hi Martin,
Just wanted to say:
On Wed, Oct 5, 2011 at 11:39 PM, Martin Morgan <mtmorgan at fhcrc.org> wrote:
> fq = FastqStreamer(<...>)
> while (length(res <- yield(fq)))
> # work, e.g., filter
That's really cool!
Then some navel gazing:
Have you thought about "inverting" this flow? Like, run the while loop
in "C-land" but pass an R expression/block/something in and have it be
evaluated within each iteration of the C/while loop?
I'm guessing calling an R function from within C code is costly, but
"while" loops in R are also slow (compared to while loops in C), so I
wonder which would win in the long run.
Just curious -- sorry if I missed some previous discussion on this topic.
Anyway, like I said -- this is really cool already.
Thanks,
-steve
--
Steve Lianoglou
Graduate Student: Computational Systems Biology
| Memorial Sloan-Kettering Cancer Center
| Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact
More information about the Bioconductor
mailing list