Hi Joseph, On Sun, Aug 22, 2010 at 1:27 PM, joseph <[email protected]> wrote: > can somebody show me how to determine the read length distribution from a > ShortReadQ or AlignedRead objects and plot them as histograms?
You can get the length of all the reads in your fastq/shortread/etc. object using the `width` function. Imagine `reads` is your ShortReadQ object, to plot the distribution of read lengths: R> hist(width(reads)) Hope that helps, -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 _______________________________________________ Bioc-sig-sequencing mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
