exportNarrowPeaks
exportNarrowPeaks.RdexportNarrowPeaks
Examples
# call some peaks:
bam <- system.file("extdata", "ex1.bam", package="Rsamtools")
peaks <- callPeaks(bam, paired=TRUE)
#> Reading signal and identifying candidate regions...
#> Identified 2 candidate regions
#> Computing significance...
#> (In the absence of a control, FDR is unlikely to be calibrated)
#> Reporting 2 regions, 2 with FDR<0.05
# save them:
filepath <- tempfile(fileext="narrowPeak")
exportNarrowPeaks(peaks, filepath)