library(motifmatchr) anno <- getPeakAnnotation(ArchRProj = proj, name = "Motif") motif_ix <- matchMotifs(anno$motifs, GR, genome = "hg38",out = "score") score <- assays(motif_ix)$motifScores apply(score,1,function(x) which(x==max(x))) colnames(score)[tail(sort(as.matrix(score)), 10)]
Preview:
downloadDownload PNG
downloadDownload JPEG
downloadDownload SVG
Tip: You can change the style, width & colours of the snippet with the inspect tool before clicking Download!
Click to optimize width for Twitter