diff hicPCA.xml @ 4:38a612a120a8 draft

planemo upload for repository https://github.com/maxplanck-ie/HiCExplorer/tree/master/galaxy/wrapper/ commit eec0a4d5a7c5ba4ec0fbd2ead8280c3d143bb9d8

--- a/hicPCA.xml	Wed Mar 07 03:26:15 2018 -0500
+++ b/hicPCA.xml	Fri Apr 27 03:31:20 2018 -0400
@@ -1,5 +1,5 @@
 <tool id="hicexplorer_hicpca" name="@BINARY@" version="@WRAPPER_VERSION@.0">
-    <description>computes the principal components for A / B compartment analysis</description>
+    <description>compute the principal components for A / B compartment analysis</description>
     <macros>
         <token name="@BINARY@">hicPCA</token>
         <import>macros.xml</import>
@@ -8,10 +8,10 @@
     <command detect_errors="exit_code"><![CDATA[
 
         hicPCA --matrix '$matrix_h5_cooler'
-       
+
         --outputFileName pca1.$outputFormat pca2.$outputFormat
         --format $outputFormat
-        
+
         && mv pca1.$outputFormat pca1
         && mv pca2.$outputFormat pca2
 
@@ -26,20 +26,20 @@
 
     </inputs>
     <outputs>
-       
-        <data name="pca1" from_work_dir="pca1" format="bigwig">
+
+        <data name="pca1" from_work_dir="pca1" format="bigwig" label="${tool.name} on ${matrix_h5_cooler.name} [${on_string}]: PC1">
             <filter>outputFormat == 'bigwig'</filter>
             <change_format>
                 <when input="outputFormat" value="bedgraph" format="bedgraph" />
             </change_format>
         </data>
-        <data name="pca2" from_work_dir="pca2" format="bigwig">
+        <data name="pca2" from_work_dir="pca2" format="bigwig"  label="${tool.name} on ${matrix_h5_cooler.name} [${on_string}]: PC2">
             <change_format>
                 <when input="outputFormat" value="bedgraph" format="bedgraph" />
             </change_format>
         </data>
 
-        
+
     </outputs>
     <tests>
         <test>
@@ -59,23 +59,38 @@
     <help><![CDATA[
 Principal component analysis
 ============================
-`hicPCA` computes two eigenvector files based on the input matrix for an A / B compartment analysis.
+
+`Lieberman-Aiden et al. (2009)`_ demonstrated that open and closed chromatin domains throughout the genome occupy different spatial compartments in the nucleus, defined as A (activate) and B (inactive) compartments.
 
-Input
+**hicPCA** computes two eigenvector files based on the input matrix for an A / B compartment analysis following the computation steps detailed by `Lieberman-Aiden et al. (2009)`_: the transformation of the contact matrix
+into an observed vs. expected matrix and consecutively a Pearson correlation matrix shows a plaid pattern. These plaid pattern are called A and B. Applying a PCA on the Pearson correlation matrix gives the eigenvectors
+and Lieberman-Aiden shows that the values of the eigenvectors correspond to the distribution of genes and with features of open and closed chromatin. In some cases the first principal component corresponds to the two
+chromosomes arms and the second eigenvector to the plaid pattern. Therefore always the first two principal components needs to be returned and investigated.
+
+_________________
+
+Usage
 -----
-- the matrix to be analysed
 
-Parameters
-__________
-- Output file format: bigwig or bedgraph
+This tool must be used on Hi-C contact matrices with large bins (over 20kb) using ``hicMergeMatrixBins`` and corrected with ``hicCorrectMatrix``. Using matrices with a too high resolution (small bins or at restriction enzyme resolution) might take several days to run (even with over 100 CPU) or will fail due to memory limitations.
+
+_________________
 
 Output
 ------
-Two files with the first and the second eigenvector.
+Two files are outputed by **hicPCA**, one with the first (pca1) and one with the second (pca2) eigenvector as bigwig or bedgraph. These files can be plotted alongside Hi-C heatmaps, gene density or external datasets such as open chromatin or histone marks enrichment using ``hicPlotTADs`` or ``hicPlotMatrix``.
+
+For example, below you can find a ``hicPlotMatrix`` of the Pearson correlation matrix derived from a contact matrix for chromosome 6 in mouse computed with ``hicTransform`` (which is part of A/B compartments computation). The optional data track at the bottom shows the first eigenvector for A/B compartment obtained using **hicPCA**.
+
+.. image:: $PATH_TO_IMAGES/hicPCA.png
+   :width: 60 %
+
+_________________
 
 | For more information about HiCExplorer please consider our documentation on readthedocs.io_
 
 .. _readthedocs.io: http://hicexplorer.readthedocs.io/en/latest/index.html
+.. _`Lieberman-Aiden et al. (2009)`: https://dx.doi.org/10.1126%2Fscience.1181369
 ]]></help>
     <expand macro="citations" />
 </tool>