peak_merge: Merge adjacent peaks within H3K4me3 or Pol II data.
In primirTSS: Prediction of pri-miRNA Transcription Start Site

Description Usage Arguments Value Examples

View source: R/peak_merge.R

Merge the adjacent segments provided as GRange object from original data. This function will merge adjacent peaks the distance between which is less than n base pairs apart and then return the merged segments.

1	peak_merge(peak, n = 250)

`peak`	A `GRange` object. The peaks to be merged from one certain ChIP-seq data, such as H3K4me3 data or Pol II data.
`n`	A number. `n` stipulates the distance(bp, base pair) between two seperate peaks within which they should be merged.

A GRanges object. The merged peaks for the following analysis to search for TSSs.

peak_df <- data.frame(chrom = c("chr1", "chr2", "chr1"),
                      chromStart = c(450, 460, 680),
                      chromEnd = c(470, 480, 710),
                      stringsAsFactors = FALSE)
peak <-  as(peak_df, "GRanges")

peak_merge(peak, n =250)