R/summary_blast.R

In cstubben/trinotateR: Trinotate annotation report summaries

summary_blast <- function(x){
   uniq <-  c( UniProt = uniqueN(x$uniprot),  genes=  uniqueN(x$gene), transcripts = uniqueN(x$transcript), proteins = uniqueN(x$protein) )
   y <- x[,.( genes=uniqueN(gene), transcripts=uniqueN(transcript), proteins = uniqueN(protein) , total =.N ), by=.(uniprot, domain, genus, name) ]
 
   y <- y[order(-y$genes, tolower(y$name)),]

    annot <- colSums(y[, .(total, genes, transcripts, proteins) ])    
   z <-  cbind(unique=uniq, annotations= annot) 
   attr(y, "counts") <- z
   message(nrow(y), " rows")
   y
}