Wednesday, March 21, 2012

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가 이뿌게 나옵니다. ㅎㅎ

No comments: