peakcall: Peak calling
In SCATE: SCATE: Single-cell ATAC-seq Signal Extraction and Enhancement

Description Usage Arguments Details Value Author(s) Examples

View source: R/peakcall.R

Peak calling function

1	peakcall(res, flank = 1, fdrcut = 1e-05)

`res`	Result matrix returned by SCATE
`flank`	Numeric variable of the number of flanking bins for each bin. For each bin, an averaged signal of itself and the flanking bins will be calculated and compared to a background distribution.
`fdrcut`	Numeric variable of FDR cutoff. Bins passing the FDR cutoff will be peaks.

This function performs peak calling for signal generated by SCATE

A list with length equal to the number of clusters. Each element is a data frame with five columns: chromosome name, starting location, ending location, FDR and signal. The data frame is ordered by FDR then by signal.

Zhicheng Ji, Weiqiang Zhou, Wenpin Hou, Hongkai Ji* <whou10@jhu.edu>

1
2
3

gr <- GRanges(seqnames="chr1",IRanges(start=seq_len(100)+1e6,end=seq_len(100)+1e6))
scateout <- SCATE(gr,clunum=156,genome='mm10')[seq_len(1000000),,drop=FALSE]
peakcall(scateout)