# HG changeset patch
# User rico
# Date 1333652925 14400
# Node ID 0b4758202a611665e65c97b80816668004c1cbad
# Parent  af482edf0b712190a4eeb73bbf443763ad4dc822
Uploaded

diff -r af482edf0b71 -r 0b4758202a61 dpmix.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/dpmix.xml	Thu Apr 05 15:08:45 2012 -0400
@@ -0,0 +1,67 @@
+<tool id="gd_dpmix" name="Admixture" version="1.0.0">
+  <description>using dynamic programming</description>
+
+  <command interpreter="python">
+    dpmix.py "$input" "$data_source" "$switch_penalty" "$ap1_input" "$ap2_input" "$p_input" "$output" "$output2" "$output2.extra_files_path" "$input.dataset.metadata.dbkey" "$input.dataset.metadata.ref" "$GALAXY_DATA_INDEX_DIR"
+    #for $individual, $individual_col in zip($input.dataset.metadata.individual_names, $input.dataset.metadata.individual_columns)
+      #set $arg = '%s:%s' % ($individual_col, $individual)
+      "$arg"
+    #end for
+  </command>
+
+  <inputs>
+    <param name="input" type="data" format="wsf" label="Dataset">
+      <validator type="unspecified_build" message="This dataset does not have a reference species and cannot be used with this tool" />
+    </param>
+    <param name="ap1_input" type="data" format="ind" label="Ancestral population 1 individuals" />
+    <param name="ap2_input" type="data" format="ind" label="Ancestral population 2 individuals" />
+    <param name="p_input" type="data" format="ind" label="Potentially admixed individuals" />
+
+    <param name="data_source" type="select" format="integer" label="Data source">
+      <option value="0" selected="true">sequence coverage</option>
+      <option value="1">estimated genotype</option>
+    </param>
+
+    <param name="switch_penalty" type="integer" min="0" value="10" label="Switch penalty" />
+  </inputs>
+
+  <outputs>
+    <data name="output" format="tabular" />
+    <data name="output2" format="html" />
+  </outputs>
+
+  <tests>
+    <test>
+      <param name="input" value="genome_diversity/test_in/sample.wsf" ftype="wsf" />
+      <param name="ap1_input" value="genome_diversity/test_in/a.ind" ftype="ind" />
+      <param name="ap2_input" value="genome_diversity/test_in/b.ind" ftype="ind" />
+      <param name="p_input" value="genome_diversity/test_in/c.ind" ftype="ind" />
+      <param name="data_source" value="0" />
+      <param name="switch_penalty" value="10" />
+
+      <output name="output" file="genome_diversity/test_out/dpmix/dpmix.tabular" />
+
+      <output name="output2" file="genome_diversity/test_out/dpmix/dpmix.html" ftype="html" compare="diff" lines_diff="2">
+        <extra_files type="file" name='dpmix.pdf' value="genome_diversity/test_out/dpmix/dpmix.pdf" compare="sim_size" delta = "10000" />
+        <extra_files type="file" name='misc.txt' value="genome_diversity/test_out/dpmix/misc.txt" />
+      </output>
+    </test>
+  </tests>
+
+  <help>
+**What it does**
+
+The user specifies two "ancestral" populations (i.e., sources for
+chromosomes) and a set of potentially admixed individuals, and chooses
+between the sequence coverage or the estimated genotypes to measure
+the similarity of genomic intervals in admixed individuals to the two
+classes of ancestral chromosomes.  The user also picks a "switch penalty",
+typically between 10 and 100.  For each potentially admixed individual,
+the program divides the genome into three "genotypes": (0) homozygous
+for the second ancestral population (i.e., both chromosomes from that
+population), (1) heterozygous, or (2) homozygous for the second ancestral
+population.  Smaller values of the switch penalty (corresponding to more
+ancient admixture events) generally lead to the reconstruction of more
+frequent changes between genotypes.
+  </help>
+</tool>