About metadata columns in Grange objects required for ChIPseeker analysis
0
0
Entering edit mode
ahua217 ▴ 10
@ahua217-24323
Last seen 3.7 years ago

Hello! I am using ChIPseeker to analyze the BED data from SEACR (a CUT & RUN analysis tool). When I used readPeakFile() function in ChIPseeker to read the BED file, it gave out such a list:

GRanges object with 18965 range and 3 metadata columns:

seqnames ranges strand | V4 V5 V6

<rle> <iranges> <rle> | <numeric> <numeric> <factor>

[1] chr10 42527196-42538077 * | 3500650 3840.71 chr10:42529404-42529405

... According to the description from SEACR (https://github.com/FredHutch/SEACR), V4 is "Total signal contained within denoted coordinates", V5 is "Maximum bedgraph signal attained at any base pair within denoted coordinates". I wonder if V4 and V5 satisfy the requirements of Grange objects for ChIPseeker. And which column was taken as the real score for ChIPseeker, V4 or V5? I am puzzled since the example data from ChIPseeker gave out a list with different metadata columns.

GRanges object with 1331 ranges and 2 metadata columns:

seqnames ranges strand | V4 V5

<rle> <iranges> <rle> | <character> <numeric>

[1] chr1 815093-817883 * | MACSpeak1 295.76

So my last question is, if my V4 and V5 are not suitable, how can I convert it into right order?

Thank you very much!

ChIPseeker • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 820 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6