# HG changeset patch # User gtrichard # Date 1565124513 14400 # Node ID a2022103f2a2221267561d5eed4c2a058530888a # Parent cbeae0f83ff84b4b95307957e4ce4e0f0ee12cfe Uploaded diff -r cbeae0f83ff8 -r a2022103f2a2 dsCompareCurves.xml --- a/dsCompareCurves.xml Mon Jul 15 17:54:31 2019 -0400 +++ b/dsCompareCurves.xml Tue Aug 06 16:48:33 2019 -0400 @@ -1,7 +1,6 @@ - compares multiple bigwigs over multiple region sets with bootstraps and - Wilcoxon rank-sum test using deepTools computeMatrix output + compares multiple genomic scores at multiple regions sets by bootstraps and per-bin distribution test dsCompareCurves @@ -178,16 +177,52 @@ - - ` + +dsCompareCurves: Comparing multiple genomic scores at multiple regions sets +=========================================================================== + + +**dsCompareCurves** assesses if multiple genomics signals (ChIP-seq, ATAC-seq, etc as .bigwig files) are significantly different or not for multiple region sets (.bed files). ``dsCompareCurves`` uses bootstraps and corrected Wilcoxon Rank-sum tests to do so. This tool uses the output of deepTools ``computeMatrix --outFileNameMatrix``, that thus must be used prior to ``dsCompareCurves``. Eitherway the regions are compared for each signal (``--comparison regions``), or all the signals are compared for each region set (``--comparison scores``). + + +_______________ + + +Example usage +------------- + - @REFERENCES@ - ]]> - +Below you can find an example of regions comparison. + +We consider two set of regions, conserved genes, and unconserved genes in Drosophila melanogaster. The want to know which histone marks from modENCODE are differentially enriched between these two regions sets. + +Prior to ``dsCompareCurves``, we used **deepTools** ``computeMatrix scale-regions`` with the two bed files as regions, and the three histone marks bigwig tracks as scores. We also added the optional parameter ``outFileNameMatrix`` in order to get the output as a tab-delimited file (.txt). + +The latter file is then used as input for the present tool in "regions", in order to compare the conserved and unconserved genes for each one of the three histone marks. You can find an image describing the output below using the said data set and the following galaxy parameters: + +.. class:: bold + +`Example parameters`_ + +.. _Example parameters: https://raw.githubusercontent.com/gtrichard/deepStats/master/docs/images/dsCompareCurves_galaxy.png + +:raw-html:`
` + +Output description +++++++++++++++++++ + +.. image:: https://raw.githubusercontent.com/gtrichard/deepStats/master/docs/images/dsCompareCurves_doc.png + :width: 800 + + ]]>