REBAR

Removing the Effects of Bias through Analysis of Residuals

Code and data accompanying the preprint:

Improving the accuracy of bulk fitness assays by correcting barcode processing biases.
Ryan Seamus McGee, Grant Kinsler, Dmitri Petrov, Mikhail Tikhonov.
bioRxiv 2023.10.31.565047; doi: https://doi.org/10.1101/2023.10.31.565047

Code

rebar.py is a python module that implements our barcode bias correction method.

You can run this bias correction method using the command line or as part of your own python scripts, as seen in the brief examples below and in the case study `run-rebar.ipynb' notebooks.

Further documentation about code usage is forthcoming. Questions regarding this code can be directed to ryansmcgee@gmail.com.

Running the method from the command line

Run rebar.py as a script from the command line, specifying the paths of input files and the desired output directory as arguments.

For example:

python rebar.py --counts ./case-study/data/kinsler-2020-preprocessed-data/counts.csv --samples ./case-study/data/kinsler-2020-preprocessed-data/samples.csv --variants ./case-study/data/kinsler-2020-preprocessed-data/variants_withNeutralGroups.csv --config ./case-study/rebar-config.json -outdir case-study/results/ 

Interfacing with the method within python scripts

Our method is implemented using an object-oriented structure that makes it easy to interface with as a "black box" or in a more "hands on" manner as desired.

For example, executing the method within your own script can be as simple as:

# Set up bias correction object
debiaser = BiasCorrection(counts='./data/kinsler-2020-preprocessed-data/counts.csv', 
                          samples='./data/kinsler-2020-preprocessed-data/samples.csv', 
                          variants='./data/kinsler-2020-preprocessed-data/variants_withNeutralGroups.csv',
                          config='./rebar-config.json',
                          outdir='./results/')

# Perform bias correction
debiaser.run()

Case Study Results

We applied our bias inference and correction method to the bulk fitness assay data from Kinsler et al. (2020) (refer to our preprint for more information). Input data, jupyter notebooks that replicate the execution of our method, results data, and figures can be found in the case-study directory.