in reply to Perl and Bioinformatics
GuysData Structures. So BioPerl gives you some pretty cool data structures that are easy to handle. Its when you run into custom structures that you get problems. For instance if I'm working with E.coli, I have ~5e6 bp of DNA - 1e7bp if I'm working on each nucleotide on both strands. How do I manage working with an analysis that needs to annotate every base, eg working with coverage from next gen analysis. Using arrays or hashes gets ugly because you will typically run out of memory. I'm not aware of an out of the box BioPerl solution, though I could stand to be corrected. You could use pack and unpack. You could use DB::File. You might even go to Berkeley DB. But the problem is general enough that it would be useful to see one or more tutorials on what to do for these larger analysis problems that are beyond simple scripts and not necessarily part of the BioPerl toolbox.
I enjoyed the node. Here is a suggested topic that would be worth pursuing for the biological crowd
MadraghRua
yet another biologist hacking perl....
|
---|
Replies are listed 'Best First'. | |
---|---|
Re^2: Perl and Bioinformatics
by BioLion (Curate) on Feb 19, 2010 at 15:27 UTC |
In Section
Tutorials