https://github.com/agusinac/run-dada2-batch
Science Score: 49.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
✓DOI references
Found 3 DOI reference(s) in README -
✓Academic publication links
Links to: zenodo.org -
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (8.0%) to scientific vocabulary
Repository
Basic Info
- Host: GitHub
- Owner: agusinac
- License: mit
- Language: R
- Default Branch: main
- Size: 37.1 KB
Statistics
- Stars: 2
- Watchers: 1
- Forks: 0
- Open Issues: 0
- Releases: 1
Metadata Files
README.md
Inspiration
Dada2 gets very slow when dealing with large sample size and read depth. The parallel_dada2.R tries to improve the denoise step when processing large sample size with large read depth on an HPC environment.
This is done by splitting the data into batches and using foreach R module to process samples in parallel with the same learning rates per batch. The user has the option to specify --parallel flag and also in case they use novaseq sequencing data via --novaseq.
The parallel_dada2 works only on single-end reads, since in our organization we merge the paired-end reads with PEAR. We find an enrichment of reads with PEAR versus dada2 denoise paired-end in Qiime2.
Mapping file format example (Required):
sample-id absolute-filepath
S103 path/to/S103.fastq.gz
S104 path/to/S104.fastq.gz
Example: Running parallel_dada2.R from command line
bash
Rscript parallel_dada2.R
--metadata mapping.tsv
--batch_n 500
--cpus 8 \
--p-trunc-q 2 \
--p-max-ee 6 \
--p-min-fold-parent-over-abundance 2 \
--p-chimera-method consensus \
> dada_report.txt
Docker
```bash docker pull agusinac/run-dada2-batch:0.0.1
docker run --rm -v $(pwd):/scripts agusinac/run-dada2-batch:0.0.1 Rscript scripts/paralleldada2.R \ --metadata mapping.tsv --batchn 500 --cpus 8 \ --p-trunc-q 2 \ --p-max-ee 6 \ --p-min-fold-parent-over-abundance 2 \ --p-chimera-method consensus \ > dada_report.txt ```
Owner
- Name: Alem Gusinac
- Login: agusinac
- Kind: user
- Repositories: 1
- Profile: https://github.com/agusinac
GitHub Events
Total
- Watch event: 2
- Push event: 9
Last Year
- Watch event: 2
- Push event: 9