# HG changeset patch # User iuc # Date 1562712369 14400 # Node ID 194e3f2c1d860ca978569d4747cdf280d58c8fa2 # Parent c97a786e8fb5df7780b61d641ad82ec508b63213 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/diffbind commit 515da67fc6c289aa7ea75a0819b5da3a902f87a2 diff -r c97a786e8fb5 -r 194e3f2c1d86 diffbind.R --- a/diffbind.R Thu Jun 07 05:21:22 2018 -0400 +++ b/diffbind.R Tue Jul 09 18:46:09 2019 -0400 @@ -57,16 +57,15 @@ bamReads=bams, bamControl=ctrls, Peaks=peaks, - Tissue=samples, # using "Tissue" column to display ids as labels in PCA plot - stringsAsFactors=FALSE) + Tissue=samples) # using "Tissue" column to display ids as labels in PCA plot } else { + sampleTable <- data.frame(SampleID=samples, Replicate=samples, Condition=groups, bamReads=bams, Peaks=peaks, - Tissue=samples, - stringsAsFactors=FALSE) + Tissue=samples) } sample = dba(sampleSheet=sampleTable, peakFormat='bed', scoreCol=opt$scorecol, bLowerScoreBetter=opt$lowerbetter) @@ -126,7 +125,7 @@ Strand=gsub("\\*", ".", strand(resSorted)), mcols(resSorted)) } -write.table(resSorted, file = opt$outfile, sep="\t", quote = FALSE, append=TRUE, row.names = FALSE) +write.table(resSorted, file = opt$outfile, sep="\t", quote = FALSE, row.names = FALSE) # Output binding affinity scores if (!is.null(opt$bmatrix)) { diff -r c97a786e8fb5 -r 194e3f2c1d86 diffbind.xml --- a/diffbind.xml Thu Jun 07 05:21:22 2018 -0400 +++ b/diffbind.xml Tue Jul 09 18:46:09 2019 -0400 @@ -1,9 +1,9 @@ - + differential binding analysis of ChIP-Seq peak data - bioconductor-diffbind - r-getopt - r-rjson + bioconductor-diffbind + r-getopt + r-rjson