Mercurial > repos > mvdbeek > r_goseq_1_22_0

diff get_length_and_gc_content.r @ 13:ea84562a5115 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 07:15:06 -0500
parents 02e88556ce1d
children 8ce951313688
line wrap: on
line diff
--- a/get_length_and_gc_content.r	Thu Feb 25 06:44:44 2016 -0500
+++ b/get_length_and_gc_content.r	Thu Feb 25 07:15:06 2016 -0500
@@ -41,7 +41,7 @@
     c(width, nGCs/width)
 }
 output <- t(sapply(split(reducedGTF, elementMetadata(reducedGTF)$gene_id), calc_GC_length))
-output <- setDT(output, keep.rownames = TRUE)[]
+output <- setDT(data.frame(output), keep.rownames = TRUE)[]
 colnames(output) <- c("#gene_id", "length", "GC")
 
 write.table(output, file=output_file, row.names=FALSE, quote=FALSE, sep="\t")
\ No newline at end of file