report.peak: Report peaks

Description Usage Arguments Details Value References See Also

View source: R/DBChIP.R

Description

Report most significant peaks

Usage

1
report.peak(test.res, FDR = NULL, FDR.method = "BH", n = 10, add.origin = TRUE, adaptive.threshold = c(0.05, 0.95))

Arguments

test.res

a list with the item test.stat, which is a data.frame of test statistics for testing non-differential binding at each site, include p-values and fold changes.

FDR

the desireable false discovery rate (FDR) level.

FDR.method

the method to control FDR. Default is "BH", the Benjamini and Hochberg (1995) method. Another option is "adaptive", which use the adaptive.threshold to estimate the number of true null hypotheses.

n

the top n differential peaks to return.

add.origin

logical. Whether to add peaks' origin information (how many and which conditions the consensus site are merged from) in the result.

adaptive.threshold

vector of two values between 0 and 1. The number of true null hypotheses is estimated as the number of p-values between these two values divided by the distance between these two values. Default is (0.05, 0.95). If it is set to (0.05, 1), it becomes the method recommended in Blanchard and Roquain (2009).

Details

The default is to return the top n differential peaks. If user specify FDR, a set of peaks under the threshold will be returned instead. The FDR is computed through the classical Benjamini & Hochberg 1995 method.

Value

a data.frame with with following fields

chr

chromosome.

pos

consensus binding position.

nsig

number of significant original binding sites that are merged into consensus site.

origin

the origin of merged binding sites.

ori.pos

the original positions of merged binding sites. The consensus position is a weighted average of original positions.

FC.condition name

the fold change comparing to the first condition.

pval

p-value for testing non-differential binding.

FDR

the q-value.

References

Benjamini, Y. and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B, 57, 289-300.

Blanchard, G. and Roquain, E. (2009). Adaptive false discovery rate control under independence and dependence. Journal of Machine Learning Research, 10, 2837-2871.

See Also

DBChIP


DBChIP documentation built on Nov. 8, 2020, 8:27 p.m.