Use the script "Preprocessing_selectRegion.py"
Please add bed files into the sub-folder "bed_file" and edit "bed_info.txt" before running this script.
Help:
usage: Preprocessing_selectRegion.py [-h] [-i INFO_BED] [-g GENOME]
[-p PARAMETERS [PARAMETERS ...]]
[-n NORMALIZE] [-o OUTPUT]
Data preprocessing and region selction for GATE
optional arguments:
-h, --help show this help message and exit
-i INFO_BED, --info_bed INFO_BED
information of aligned bed files to be inputed,must be
stored in the folder "bed_file"
-g GENOME, --genome GENOME
genome information for input data, can be
mm9/hg19/hg18/
-p PARAMETERS [PARAMETERS ...], --parameters PARAMETERS [PARAMETERS ...]
parameters [m,n] to select regions. those regions with
normalized signal no less than 'm' in at least one
time point for at least 'n' epi-marks were selected,
default:[5,5]
-n NORMALIZE, --normalize NORMALIZE
normalize total read counts into this number,
default:1E8
-o OUTPUT, --output OUTPUT
output txt file
None
The BED files in folder "/bed_file/" should be the ones with names in the third column of "bed_info.txt".
The output file of this script can be the input for GATE software.
Use the output of step 1 as input file and run the program following the steps here.
- In the Preprocessing_selectRegion.py, one of the module xplib.Annotation is missing.
Answer: xplib is from the BAM2X tools developed by my previous lab mate Xiaopeng Zhu. I have a forked version of BAM2X here: https://github.com/yu68/bam2x.
You can clone this bam2x repository and then add "/the_folder_you_put_bam2x/lib/" to your python path.