1, web software: WebLogo: http://weblogo.berkeley.edu/
2, R package: "seqLogo": http://bioconductor.org/packages/2.2/bioc/html/seqLogo.html
ex:
source("http://bioconductor.org/biocLite.R")
biocLite("seqLogo")
library("seqLogo")
mFile <- system.file("Exfiles/pwm1", package="seqLogo")
m <- read.table(mFile)
pwm <- makePWM(m)
seqLogo(pwm)
沒有留言:
張貼留言