readJaspar: Read BED file of TFBS scores from JASPAR into GRange object

In GabrielHoffman/tfbsDB: Store and plot TFBS sites from HOCOMOCO v11 scored by FIMO

Description

Read BED file of TFBS scores from JASPAR into GRange object

Usage

readJaspar(grQuery, bedFile, maxP = 1e-04)

Arguments

grQuery	GRange object of query interval
bedFile	BED file that is bgzipped and tabix'd.
maxP	maximum p-value of TFBS returned. score corresponds to p-value: 200 -> 10e-2; 300 -> 10e-3

Value

GRange of TFBS locations where scores are transformed p-values

Examples

library(GenomicRanges)
library(tfbsDB)

# range of TFBS to be loaded from file
grQuery = GRanges("chr1", IRanges(10280, 10410))

# example data from package
bedFile = file.path(system.file("extdata/", package = "tfbsDB"), "JASPAR2018_hg19_example.bed.gz")

# read from file
gr = readJaspar( grQuery, bedFile )

# show TFBS locations
plotTFBSdb( gr )

# color locations by p-value of motif match
plotTFBSdb( gr, colorByP=TRUE)

GabrielHoffman/tfbsDB documentation built on June 23, 2020, 5:29 a.m.