anl_deseq <- function(h,...){
insert(st,"Statistical analysis : DESeq2", do.newline = TRUE )
dir <- getwd()
if(!dir.exists(paste(dir, "/stat/plot", sep = ""))){
dir.create(paste(dir, "/stat/plot", sep = ""), recursive = TRUE)}
tcount <- read.delim(file.path(dir, "rc","filtered_readcount_trnas.txt"))
ccount <- read.delim(file.path(dir, "rc","filtered_readcount_isodecoders.txt"))
acount <- read.delim(file.path(dir, "rc","filtered_readcount_isoacceptors.txt"))
t_out <- deseq(tcount)
c_out <- deseq(ccount)
a_out <- deseq(acount)
stat_trna[] <- t_out
stat_codon[] <- c_out
stat_aa[] <- a_out
write.table(t_out, "./stat/stat_trna_list.txt", sep="\t", quote = FALSE, row.names = F)
write.table(c_out, "./stat/stat_isodecoder_list.txt", sep="\t", quote = FALSE, row.names = F)
write.table(a_out, "./stat/stat_isoacceptor_list.txt", sep="\t", quote = FALSE, row.names = F)
insert(st, " ", do.newline = TRUE)
insert(st,"Done : Statistical analysis. Set the threshold value.", do.newline = TRUE )
insert(st, ".", do.newline = TRUE)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.