A bare-bones interface to the trained LanceOTron (LoT) model from the command line.
pip install lanceotron
- Clone the repository.
- Install dependencies with pip.
- Install the package.
- Run tests to ensure that everything is working.
git clone git@github.com:Chris1221/lanceotron.git; cd lanceotron # Step 1
pip install -r requirements.txt # Step 2
pip install -e . # Step 3
python -m unittest
To see available commands, use the --help
flag.
lanceotron --help
To call peaks from a bigWig track, use the callPeaks
command.
Option | Description | Default |
---|---|---|
file | BigWig Track to analyse | |
-t, --threshold | Threshold for selecting candidate peaks | 4 |
-w, --window | Window size for rolling mean to select candidate peaks | 400 |
-f, --folder | Output folder | "./" |
--skipheader | Skip writing the header | False |
To call peaks from a bigWig track with an input file, use the callPeaks_Input
command.
Option | Description | Default |
---|---|---|
file | BigWig track to analyse | |
-i, --input | Control input track to calculate significance of peaks | |
-t, --threshold | Threshold for selecting candidate peaks | 4 |
-w, --window | Window size for rolling mean to select candidate peaks | 400 |
-f, --folder | Output folder | "./" |
--skipheader | Skip writing the header | False |
To score the peaks in an existing Bed file, use the scoreBed
command.
Option | Description | Default |
---|---|---|
file | BigWig Track to analyse | |
-b, --bed | Bed file of regions to be scored | |
-f, --folder | Output folder | "./" |
--skipheader | Skip writing the header | False |
There is a basic bigWig file included in the test
subdirectory. To try out the caller, execute it on this file.
lanceotron callPeaks test/chr22.bw -f output_folder
@article {Hentges2021.01.25.428108,
author = {Hentges, Lance D. and Sergeant, Martin J. and Downes, Damien J. and Hughes, Jim R. and Taylor, Stephen},
title = {LanceOtron: a deep learning peak caller for ATAC-seq, ChIP-seq, and DNase-seq},
year = {2021},
doi = {10.1101/2021.01.25.428108},
publisher = {Cold Spring Harbor Laboratory},
URL = {https://www.biorxiv.org/content/early/2021/01/27/2021.01.25.428108},
journal = {bioRxiv}
}
To serve the documentation locally, use
python -m mkdocs serve
Please raise an issue if there is anything you wish to ask or contribute.