R에서 Seq Logo 사용

출처: Manual 사이트

여기저기 기웃기웃거리면서 pdf 보던중

R에서 Seq Logo를 그릴 수 있어서 끄적끄적... ㅎㅎ

library(seqLogo)

#seqLogo란 라이브러리 설치는 해주시고..

pwm <- PWM(DNAStringSet(c("GGT","GGT","GCA","GGA")))

seqLogo(t(t(pwm)*1/colSums(pwm)))

그러면 seqLogo가 이뿌게 나옵니다. ㅎㅎ