Mercurial > repos > mvdbeek > r_goseq_1_22_0

changeset 6:d4b5942ed347 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 05:49:08 -0500
parents f9b964d1d386
children 15ce6435ab83
files get_length_and_gc_content.r
diffstat 1 files changed, 2 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/get_length_and_gc_content.r	Thu Feb 25 05:47:38 2016 -0500
+++ b/get_length_and_gc_content.r	Thu Feb 25 05:49:08 2016 -0500
@@ -1,4 +1,5 @@
 # originally by Devon Ryan, https://www.biostars.org/p/84467/
+sink(stdout(), type = "message")
 
 library(GenomicRanges)
 library(rtracklayer)
@@ -19,7 +20,7 @@
 output = args$output
 
 #Load the annotation and reduce it
-GTF <- import.gff(GTFfile, format="gtf", genome=NA, asRangedData=F, feature.type="exon")
+GTF <- import.gff(GTFfile, format="gtf", genome=NA, feature.type="exon")
 grl <- reduce(split(GTF, elementMetadata(GTF)$gene_name))
 reducedGTF <- unlist(grl, use.names=T)
 elementMetadata(reducedGTF)$gene_name <- rep(names(grl), elementLengths(grl))