p value after enrichPeakOverlap analysis in ChIPseeker
0
0
Entering edit mode
mwong043 • 0
@mwong043-23926
Last seen 3.5 years ago
Woodbridge

Hi. I am trying to see whether the peaks in my TF ChIP-seq experiment overlap with the peaks of different epigenetic marks. However, it gave me the same p value even though the number of overlapped peaks were different.

I used the enrichPeakOverlap function as below:

H3K4me2 overlap test

enrichPeakOverlap(queryPeak = "IntersectH3K4me2.narrowPeak", targetPeak = "TargetChIP_peaks.narrowPeak", nShuffle=5000, TxDb = TxDb.Hsapiens.UCSC.hg38.knownGene , pAdjustMethod="BH", chainFile=NULL)

permutation test of peak overlap... 2020-07-26 03:56:44 PM |======================================================================| 100% qSample tSample qLen tLen NOL 1 IntersectH3K4me2.narrowPeak TargetChIPpeaks.narrowPeak 407400 15559 5535

pvalue p.adjust 1 0.00019996 0.00019996

H3K4me3 overlap test

enrichPeakOverlap(queryPeak = "IntersectH3K4me3.narrowPeak", targetPeak = "TargetChIP_peaks.narrowPeak", nShuffle=5000, TxDb = TxDb.Hsapiens.UCSC.hg38.knownGene , pAdjustMethod="BH", chainFile=NULL)

permutation test of peak overlap... 2020-07-26 03:56:44 PM |======================================================================| 100% qSample tSample qLen tLen NOL 1 IntersectH3K4me3.narrowPeak TargetChIPpeaks.narrowPeak 145835 15559 766

pvalue p.adjust 1 0.00019996 0.00019996

The values are the same for other overlapping other epigenetic marks with my Target peaks. It either said 0.00019996 or 1 I am wondering is it normal? Thanks.

ChIPseeker • 1.1k views
ADD COMMENT
0
Entering edit mode

Hi mwong,

I have encountered the same thing. Do you have an idea of this after all?

Thanks!

Kylie

ADD REPLY

Login before adding your answer.

Traffic: 698 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