SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Sample Prep / Library Generation (http://seqanswers.com/forums/forumdisplay.php?f=25)
-   -   bioanalyzeR: New R package for data from the Agilent Bioanalyzer and TapeStation (http://seqanswers.com/forums/showthread.php?t=93661)

jwfoley 03-31-2020 03:36 AM

bioanalyzeR: New R package for data from the Agilent Bioanalyzer and TapeStation
 
I wrote an R package for importing electrophoresis data from the Agilent Bioanalyzer and TapeStation and graphing the results with ggplot2 or analyzing them in more depth:

package on GitHub
vignette

One particular application is for sequencing library QC, where the variable of interest is the size distribution of your sample by molarity, but the Agilent software only shows you total fluorescence, which is proportional to mass instead of molarity and therefore gives you a misleading representation.

For example, here's the same library shown with fluorescence vs. log-scaled molecule length (analogous to the Agilent software's display) and then with molarity/bp vs. linear-scaled length (so area under the curve is molarity):

https://i.imgur.com/Zk3LcPF.png

And the package includes a function to automatically calculate the molar ratio of good inserts to bad adapter dimers in Illumina sequencing libraries.

This is a work in progress so I would be grateful for any feedback.

luc 03-31-2020 06:50 PM

Thanks!!!
This is a neat idea.

Marc_Jones 04-01-2020 04:26 AM

Made it to CRAN yet?

jwfoley 04-01-2020 09:24 AM

Quote:

Originally Posted by Marc_Jones (Post 232002)
Made it to CRAN yet?

Considering that or Bioconductor. It doesn't have any Bioconductor dependencies so CRAN might be easier for everyone involved, but on the other hand some of the homebrew data structures might be better implemented with Bioconductor tools and I haven't looked into that yet.


All times are GMT -8. The time now is 08:33 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.