process_candidate_sites: Obtains and filters candidate sites from FIMO result

In kevinlkx/TOP: Predict transcription factor occupancy using DNase- or ATAC-seq data

Obtains and filters candidate sites from FIMO result

Description

Gets candidate sites from FIMO motif matches and add flanking regions, and filters the candidate sites by different thresholds, and filters out sites in blacklist regions.

Usage

process_candidate_sites(
  fimo_file,
  flank = 100,
  thresh_pValue = 1e-05,
  thresh_pwmscore = 0,
  blacklist_file = NULL,
  mapability_file = NULL,
  thresh_mapability = 0.8,
  bigWigAverageOverBed_path = "bigWigAverageOverBed"
)

Arguments

fimo_file	FIMO result ‘.txt’ file
flank	Flanking region (bp) around motif matches (default: 100)
thresh_pValue	FIMO p-value threshold (default: 1e-5)
thresh_pwmscore	FIMO PWM score threshold (default: 0)
blacklist_file	Filename of the blacklist regions (default: NULL)
mapability_file	Filename of the mapability reference file
thresh_mapability	Mapability threshold (default: 0.8, candidate sites need to be mapable at least 80% positions).
bigWigAverageOverBed_path	Path to bigWigAverageOverBed executable (only needed when filtering mapability).

Value

A data frame of processed candidate binding sites passing the threshold filtering.

Examples

## Not run: 
sites <- process_candidate_sites(fimo_file='fimo.txt',
                                 thresh_pValue=1e-5,
                                 blacklist_file='blacklist.hg38.bed.gz')

## End(Not run)

kevinlkx/TOP documentation built on May 4, 2022, 6:30 p.m.