# HG changeset patch # User modencode-dcc # Date 1358548101 18000 # Node ID 740902fb1c5abe1a55c75deedd8f6ed9effdaa93 # Parent 82b54b51b6cfa15cc801af8de79bacac4cdb6812 Uploaded diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/.DS_Store Binary file genome_tables/.DS_Store has changed diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/._.DS_Store Binary file genome_tables/._.DS_Store has changed diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/genome_table.dmel.r5.32.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/genome_tables/genome_table.dmel.r5.32.txt Fri Jan 18 17:28:21 2013 -0500 @@ -0,0 +1,15 @@ +YHet 347038 +dmel_mitochondrion_genome 19517 +2L 23011544 +X 22422827 +3L 24543557 +4 1351857 +2R 21146708 +3R 27905053 +Uextra 29004656 +2RHet 3288761 +2LHet 368872 +3LHet 2555491 +3RHet 2517507 +U 10049037 +XHet 204112 diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/genome_table.human.hg18.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/genome_tables/genome_table.human.hg18.txt Fri Jan 18 17:28:21 2013 -0500 @@ -0,0 +1,25 @@ +chr1 247249719 +chr2 242951149 +chr3 199501827 +chr4 191273063 +chr5 180857866 +chr6 170899992 +chr7 158821424 +chr8 146274826 +chr9 140273252 +chr10 135374737 +chr11 134452384 +chr12 132349534 +chr13 114142980 +chr14 106368585 +chr15 100338915 +chr16 88827254 +chr17 78774742 +chr18 76117153 +chr19 63811651 +chr20 62435964 +chr21 46944323 +chr22 49691432 +chrX 154913754 +chrY 57772954 +chrM 16571 diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/genome_table.human.hg19.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/genome_tables/genome_table.human.hg19.txt Fri Jan 18 17:28:21 2013 -0500 @@ -0,0 +1,25 @@ +chr1 249250621 +chr2 243199373 +chr3 198022430 +chr4 191154276 +chr5 180915260 +chr6 171115067 +chr7 159138663 +chr8 146364022 +chr9 141213431 +chr10 135534747 +chr11 135006516 +chr12 133851895 +chr13 115169878 +chr14 107349540 +chr15 102531392 +chr16 90354753 +chr17 81195210 +chr18 78077248 +chr19 59128983 +chr20 63025520 +chr21 48129895 +chr22 51304566 +chrX 155270560 +chrY 59373566 +chrM 16571 diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/genome_table.mm9.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/genome_tables/genome_table.mm9.txt Fri Jan 18 17:28:21 2013 -0500 @@ -0,0 +1,22 @@ +chr1 197195432 +chr2 181748087 +chr3 159599783 +chr4 155630120 +chr5 152537259 +chr6 149517037 +chr7 152524553 +chr8 131738871 +chr9 124076172 +chr10 129993255 +chr11 121843856 +chr12 121257530 +chr13 120284312 +chr14 125194864 +chr15 103494974 +chr16 98319150 +chr17 95272651 +chr18 90772031 +chr19 61342430 +chrX 166650296 +chrY 15902555 +chrM 16299 diff -r 82b54b51b6cf -r 740902fb1c5a genome_tables/genome_table.worm.ws220.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/genome_tables/genome_table.worm.ws220.txt Fri Jan 18 17:28:21 2013 -0500 @@ -0,0 +1,7 @@ +I 15072423 +II 15279345 +III 13783700 +IV 17493793 +V 20924149 +X 17718866 +MtDNA 13794 \ No newline at end of file diff -r 82b54b51b6cf -r 740902fb1c5a idrPlotWrapper.sh --- a/idrPlotWrapper.sh Fri Jan 18 17:28:12 2013 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,39 +0,0 @@ -#!/bin/bash - -# idrPlotWrapper.sh -# OICR: Kar Ming Chu -# July 2012 - -# BASH wrapper for batch-consistency-plot.r (part of the IDR package) -# For use with Galaxy - -# Usage of batch-consistency-plot.r: Rscript batch-consistency-plot-merged.r [npairs] [output.dir] [input.file.prefix 1, 2, 3 ...] -# npairs - will be a constant, since Galaxy requires explicit control over input and output files - -# Usage of THIS SCRIPT: ./idrPlotWrapper.sh em uri outputfile -# em - em.sav file provided by Galaxy -# uri - uri.sav file provided by Galaxy -# outputfile - output file name specified by Galaxy - -main() { - EM="${1}" # absolute file path to em.sav file, provided by Galaxy - URI="${2}" # absolute file parth to uri.sav file, provided by Galaxy - OUTFILE="${3}" # name of desired output file - - cp "${EM}" ./idrPlot-em.sav # cp to this directory and rename so they can be found by idrPlot - cp "${URI}" ./idrPlot-uri.sav - - Rscript /mnt/galaxyTools/galaxy-central/tools/modENCODE_DCC_tools/idr/batch-consistency-plot.r 1 ./idrPlot idrPlot - - # convert post script to pdf file - ps2pdf ./idrPlot-plot.ps ./idrPlot-plot.pdf - - # rename to output file name - mv ./idrPlot-plot.pdf "${OUTFILE}" - - # clean up - rm idrPlot-em.sav - rm idrPlot-uri.sav -} - -main "${@}" diff -r 82b54b51b6cf -r 740902fb1c5a idrToolDef.xml --- a/idrToolDef.xml Fri Jan 18 17:28:12 2013 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,61 +0,0 @@ - - - - Consistency Analysis on a pair of narrowPeak files - batch-consistency-analysis.r $input1 $input2 $halfwidth $overlap $option $sigvalue $gtable $rout $aboveIDR $ratio $emSav $uriSav - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -Reproducibility is essential to reliable scientific discovery in high-throughput experiments. The IDR (Irreproducible Discovery Rate) framework is a unified approach to measure the reproducibility of findings identified from replicate experiments and provide highly stable thresholds based on reproducibility. Unlike the usual scalar measures of reproducibility, the IDR approach creates a curve, which quantitatively assesses when the findings are no longer consistent across replicates. In layman's terms, the IDR method compares a pair of ranked lists of identifications (such as ChIP-seq peaks). These ranked lists should not be pre-thresholded i.e. they should provide identifications across the entire spectrum of high confidence/enrichment (signal) and low confidence/enrichment (noise). The IDR method then fits the bivariate rank distributions over the replicates in order to separate signal from noise based on a defined confidence of rank consistency and reproducibility of identifications i.e the IDR threshold. For more information on IDR, see https://sites.google.com/site/anshulkundaje/projects/idr - - -