# HG changeset patch # User rnateam # Date 1552460799 14400 # Node ID 3ff6333973aa9e7f27e048e5329022a5806292f9 # Parent 4fbd735886ddcc6a86d99feed2c23609e76629d6 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/chipseeker commit d83142dbe2432fcb0f56dcd6311a05c061628ecc diff -r 4fbd735886dd -r 3ff6333973aa chipseeker.R --- a/chipseeker.R Tue Dec 04 10:02:01 2018 -0500 +++ b/chipseeker.R Wed Mar 13 03:06:39 2019 -0400 @@ -17,6 +17,8 @@ make_option(c("-d","--downstream"), type="integer", help="TSS downstream region"), make_option(c("-F","--flankgeneinfo"), type="logical", help="Add flanking gene info"), make_option(c("-D","--flankgenedist"), type="integer", help="Flanking gene distance"), + make_option(c("-j","--ignoreUpstream"), type="logical", help="Ignore upstream"), + make_option(c("-k","--ignoreDownstream"), type="logical", help="Ignore downstream"), make_option(c("-f","--format"), type="character", help="Output format (interval or tabular)."), make_option(c("-p","--plots"), type="logical", help="PDF of plots."), make_option(c("-r","--rdata"), type="logical", help="Output RData file.") @@ -31,15 +33,36 @@ down = args$downstream format = args$format +if (!is.null(args$flankgeneinfo)) { + flankgeneinfo <- TRUE +} else { + flankgeneinfo <- FALSE +} + +if (!is.null(args$ignoreUpstream)) { + ignoreUpstream <- TRUE +} else { + ignoreUpstream <- FALSE +} + +if (!is.null(args$ignoreDownstream)) { + ignoreDownstream <- TRUE +} else { + ignoreDownstream <- FALSE +} + peaks <- readPeakFile(peaks) # Make TxDb from GTF txdb <- makeTxDbFromGFF(gtf, format="gtf") -if (!is.null(args$flankgeneinfo)) { - peakAnno <- annotatePeak(peaks, TxDb=txdb, tssRegion=c(-up, down), addFlankGeneInfo=args$flankgeneinfo, flankDistance=args$flankgenedist) -} else { - peakAnno <- annotatePeak(peaks, TxDb=txdb, tssRegion=c(-up, down)) -} + +# Annotate peaks +peakAnno <- annotatePeak(peaks, TxDb=txdb, + tssRegion=c(-up, down), + addFlankGeneInfo=flankgeneinfo, + flankDistance=args$flankgenedist, + ignoreUpstream=ignoreUpstream, + ignoreDownstream=ignoreDownstream) # Add gene name features <- import(gtf, format="gtf") diff -r 4fbd735886dd -r 3ff6333973aa chipseeker.xml --- a/chipseeker.xml Tue Dec 04 10:02:01 2018 -0500 +++ b/chipseeker.xml Wed Mar 13 03:06:39 2019 -0400 @@ -1,8 +1,8 @@ - + for ChIP peak annotation and visualization - bioconductor-chipseeker - r-optparse + bioconductor-chipseeker + r-optparse /dev/null | grep -v -i "WARNING: ")", optparse version" $(R --vanilla --slave -e "library(optparse); cat(sessionInfo()\$otherPkgs\$optparse\$Version)" 2> /dev/null | grep -v -i "WARNING: ") @@ -25,14 +25,27 @@ -G '$gtf' -u $adv.upstream -d $adv.downstream + #if $adv.flankgeneinfo: -F $adv.flankgeneinfo - -D $adv.flankgenedist + #end if + + -D $adv.flankgenedist + + #if $adv.ignoreUpstream: + -j $adv.ignoreUpstream #end if + + #if $adv.ignoreDownstream: + -k $adv.ignoreDownstream + #end if + -f $format + #if $pdf: -p $pdf #end if + -r $rdata ]]> @@ -70,6 +83,8 @@ + +