motifs2View: Motifs Two View

Description Usage Examples

Description

Visualize the similarity between two motifs, Very ineficent,

Usage

1
motifs2View(m1, m2, reg, Sequences, nearHeights = FALSE)

Examples

1
2
3
4
5
6
7
FastaFile<-"~/Dropbox/UTX-Alex/jan/combined.fasta"
Sequences <- readDNAStringSet(FastaFile, "fasta")
data<-loadHeightFile("~/Dropbox/UTX-Alex/jan/combined_heights.bed")$data
reg<-as.matrix(ascore(data,c(1,1,1),c("top","bottom","middle"),c(1,1,1)))
m<-motifs2View("CANNTG","TGACCT",reg[,3],Sequences)
k<-motifs2View("CANNTG","GATAAG",reg[,1],Sequences)
t0<-unlist(lapply(seq(min(k),max(k)),function(i)length(which(k==i))))

alexjgriffith/CCCA documentation built on May 10, 2019, 8:52 a.m.