Mercurial > repos > mvdbeek > r_goseq_1_22_0

comparison get_length_and_gc_content.r @ 16:8ce951313688 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/goseq_1_22_0 commit fdd0811efc61c31f88ff17096fbe8ee8cfacd766-dirty
author mvdbeek
date Thu, 25 Feb 2016 08:49:20 -0500
parents ea84562a5115
children
comparison
equal deleted inserted replaced
15:e3b1dd2f7f70 16:8ce951313688
43 output <- t(sapply(split(reducedGTF, elementMetadata(reducedGTF)$gene_id), calc_GC_length)) 43 output <- t(sapply(split(reducedGTF, elementMetadata(reducedGTF)$gene_id), calc_GC_length))
44 output <- setDT(data.frame(output), keep.rownames = TRUE)[] 44 output <- setDT(data.frame(output), keep.rownames = TRUE)[]
45 colnames(output) <- c("#gene_id", "length", "GC") 45 colnames(output) <- c("#gene_id", "length", "GC")
46 46
47 write.table(output, file=output_file, row.names=FALSE, quote=FALSE, sep="\t") 47 write.table(output, file=output_file, row.names=FALSE, quote=FALSE, sep="\t")
48
49 sessionInfo()