added logoplot_example.R and ga_customers.csv

This commit is contained in:
Tanushree Tunstall 2022-04-05 14:52:08 +01:00
parent c647773520
commit 67d9e6160a
2 changed files with 10015 additions and 0 deletions

18
logolpot_example.R Normal file
View file

@ -0,0 +1,18 @@
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")