# HG changeset patch # User mvdbeek # Date 1456402506 18000 # Node ID ea84562a5115ce60b369d8417a4d12374bc468fc # Parent 02e88556ce1d6c9af033e6fdfc505867026832b8 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/goseq_1_22_0 commit fdd0811efc61c31f88ff17096fbe8ee8cfacd766-dirty diff -r 02e88556ce1d -r ea84562a5115 get_length_and_gc_content.r --- 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