Skip to contents

Regulatory peak motif enrichment analysis

Usage

SingleMotifAnalysis(peakfile, Species)

Arguments

peakfile

ATAC-Seq expression matrix or bed file

Species

Select the species, Homo or Mus.

Value

Motif result

Examples

library(linkage)
library(LinkageData)
peakpath <- system.file("extdata","ENSG00000000419.rdata",package = "LinkageData")
load(peakpath)
motif <- SingleMotifAnalysis(PeakFile, "Homo")
head(motif)
#>           ID  name seqnames    start      end width strand     score
#> 1   MA0030.1 FOXF2    chr20 51077124 51077137    14      + 11.398136
#> 1.1 MA0030.1 FOXF2    chr20 51120047 51120060    14      - 13.015671
#> 2   MA0031.1 FOXD1    chr20 51077129 51077136     8      + 11.972521
#> 2.1 MA0031.1 FOXD1    chr20 51120048 51120055     8      - 12.972000
#> 2.2 MA0031.1 FOXD1    chr20 51359554 51359561     8      + 11.568211
#> 3   MA0051.1  IRF2    chr20 51117554 51117571    18      -  9.042872