One of the last steps in my viral metagenomic pipeline is to rarefy my shared files to ensure even sampling across samples. This process is slightly different from 16S pipelines because the read counts need to be normalized based on bin lengths after subsampling but before downstream modeling.
In this repo, I have provided a basic script that I've used for the subsampling process in addition to the expected input/output files. The script is divided up into three main parts as seen below:
For this code review, I need help with the following two areas:
Theme | Goals | Teams |
---|---|---|
Defensive programming | 1. Validate inputs 2. Add progress messages 3. Add warning messages |
1. Pat, Josh, and Samara |
General code improvement | 1. Increase speed 2. Increase efficiency 3. Reduce complexity |
1. Nick and Lucas 2. Kelly and Sarah |
Groups will be given 30 min to work on their tasks and the last 20 min (5 min/group) will be used for discussion. Please create a branch of this repo and submit your work as a pull request when finished.