Motif histogram
1 2 |
1 2 3 4 5 6 7 8 9 10 11 | values<-loadHeightFile(heightFile)
data<-values$data
reg<-ascore(data,1,"top",3)
test<-readDNAStringSet(fastaFile,use.names = TRUE)
motifs<-as.matrix(read.table("data/normal_not_abnormal_motifs"))
mList<-unlist(lapply(c(motifs,addmotifs),IUPACtoBase))
cList<-unlist(lapply(lapply(c(motifs,"CGNNGC"),IUPACtoBase),compliment))
locationsM<-lapply(mList,grep,test)
locationsC<-lapply(cList,grep,test)
l<-length(cList)
motifHist(mList,cList,locationsM,locationsC,4,l,reg)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.