https://github.com/Zhang-lab/BeCorrect

https://github.com/Zhang-lab/BeCorrect

Last synced: about 1 month ago
JSON representation

• Host: GitHub
• URL: https://github.com/Zhang-lab/BeCorrect
• Owner: Zhang-lab
• Created: 2019-05-16T16:41:43.000Z (about 5 years ago)
• Default Branch: master
• Last Pushed: 2021-03-12T16:25:40.000Z (over 3 years ago)
• Last Synced: 2024-02-12T23:48:27.689Z (5 months ago)
• Language: C++
• Size: 10.7 KB
• Stars: 13
• Watchers: 3
• Forks: 3
• Open Issues: 3
• Metadata Files:
  • Readme: README.md

Lists

• awesome-atac-analysis - BeCorrect - seq signal density for visualization purposes (Uncategorized)

README

        
# BeCorrect

*New in v1.1.0: BeCorrect now accepts decimal values for "counts" table to work with edgeR normalization, quartile normalization, etc.

BeCorrect (v.1.1.0) is distributed as a standalone .cpp file.  BeCorrect can be compiled in either linux or Mac OS X using the following command:

g++ BeCorrect.cpp -o BeCorrect

Usage for BeCorrect can be displayed by running BeCorrect with no parameters supplied:

BeCorrect (v1.1.0): A program to rescale bedgraph files to visually match batch correction using RUVSeq.

  Usage: BeCorrect    1.bg ... n.bg

Creates batch corrected bedgraph files based on batch correction performed by RUVseq on counts table for CHIPseq/ATACseq peaks.  For CHIP/ATAC, raw counts, adjusted counts, and bedgraph files should be sorted by chromosome then start location.

Sample counts file for CHIP/ATAC should be as follows:

chr	start	end	name1	...	nameN

Name1	start1	end1	Val1.1	...	ValN.1

Name2	start2	end2	Val1.2	...	ValN.2

A full testing data set is available at http://regmedsrv1.wustl.edu/Public_SPACE/pgontarz/Public_html/SPACE/pgontarz/BeCorrect/.  To use this testing dataset, download and decompress using the command:

tar -xzvf BeCorrect_Test_Data.tar.gz

The decompressed directory contains 11 files: 8 bedgraph (.bg) files, a raw counts file, an adjusted counts file, and a chromosome names file.  From the directory containing the files, these can be used to run BeCorrect with the following command:

BeCorrect Raw_Counts.txt Adjusted_Counts_k3.txt rn6.names.txt 14_AH_raw.bg 15_AH_raw.bg 21_AH_raw.bg 22_AH_raw.bg 14_DG_raw.bg 15_DG_raw.bg 21_DG_raw.bg 22_DG_raw.bg

*Note: The .bg files must be supplied in the same order as the columns in the counts table.  BeCorrect does not check that files are supplied in the correct order.

It has been our experience that for a full-sized bedgraph file on a mammalian genome, BeCorrect requires ~1 min per file.