18 lines
783 B
R
18 lines
783 B
R
library(Logolas)
|
|
sequence <- c("CTATTGT","CTCTTAT","CTATTAA","CTATTTA", "CTATTAT","CTTGAAT",
|
|
"CTTAGAT","CTATTAA","CTATTTA","CTATTAT", "CTTTTAT","CTATAGT",
|
|
"CTATTTT","CTTATAT","CTATATT","CTCATTT", "CTTATTT","CAATAGT",
|
|
"CATTTGA","CTCTTAT","CTATTAT","CTTTTAT", "CTATAAT","CTTAGGT",
|
|
"CTATTGT","CTCATGT","CTATAGT", "CTCGTTA","CTAGAAT","CAATGGT")
|
|
logomaker(sequence,type = "Logo")
|
|
logomaker(sequence,type = "EDLogo")
|
|
bg <- apply(sequence, 1, function(x) return(median(x)))
|
|
data("N_Glycosyl_sequences")
|
|
logomaker(sequence,type = "EDLogo")
|
|
|
|
bg <- apply(N_Glycosyl_sequences, 1, function(x) return(median(x)))
|
|
bg <- bg/sum(bg)
|
|
bg
|
|
|
|
logomaker(N_Glycosyl_sequences, type = "EDLogo", bg=bg)
|
|
logomaker(N_Glycosyl_sequences, type = "EDLogo")
|