Hi all, I was wondering if anyone could explain the memory usage behaviour I see when subsetting a ShortReadQ object.
If I create a ShortReadQ with 500 reads, each of length 200bp it uses ~230Kb of memory. If I do the same, but store 50,000 reads the memory usage goes up to ~21.5Mb - this seems entirely reasonable. However, if I then create a third object by subsetting the first 500 reads from the larger ShortReadQ it comes in at ~19.5Mb, which seems less understandable. I was hoping for some insight into why this might be. I've put some code in a Gist here that can run this example: https://gist.github.com/grimbough/ca9e668a5772dc00cafb Thanks, Mike -- Mike Smith Research Associate Statistics & Computational Biology Laboratory Cambridge University [[alternative HTML version deleted]] _______________________________________________ Bioc-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/bioc-devel
