This script generates volcano plots with a false-discovery rate cutoff from sgRNA-level phenotypes from CRISPR-based screens.
Required: Python2, Numpy, Pandas
input file(s): .l2es files
Usage:
python volcano_plot_l2es_FDR.py PATH_of_L2ES PATH_for_OUTPUT
The script will ask users to specify the counts threshold, FDR rate (typically 0.05), figure name, and file path for a list of genes to label (for no gene labels: just press enter).
Outputs:
a volcano plot and hits files with phenotype, pvalue and product calculated.
Created by Ruilin Tian, Kampmann lab (2019).