--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/hyphy_cfel.xml	Tue Apr 20 10:26:31 2021 +0000
@@ -0,0 +1,67 @@
+<?xml version="1.0"?>
+<tool id="hyphy_cfel" name="HyPhy-CFEL" version="@VERSION@+galaxy0" profile="19.09">
+    <description>Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements"/>
+    <command detect_errors="exit_code"><![CDATA[
+        ln -s '$input_file' cfel_input.fa &&
+        ln -s '$input_nhx' cfel_input.nhx &&
+        hyphy contrast-fel
+            --code '$gencodeid'
+            #for branch in str($branch_set).split(','):
+                --branch-set '$branch'
+            #end for
+            --srv '$srv'
+            --permutations '$permutations'
+            --pvalue '$pvalue'
+            --qvalue '$qvalue'
+            --output '$cfel_output'
+            cfel_input.fa
+            cfel_input.nhx
+            > '$cfel_log'
+        @ERRORS@
+    ]]></command>
+    <inputs>
+        <expand macro="inputs"/>
+        <expand macro="gencode"/>
+        <param argument="--branch-set" type="select" multiple="true" min="2" max="2" label="Set of branches to test">
+            <option value="Internal branches">Internal branches</option>
+            <option value="Terminal branches">Terminal branches</option>
+            <option value="Random set of branches">Random set of branches</option>
+            <option value="Unlabeled branches">Unlabeled branches</option>
+        </param>
+        <param argument="--permutations" type="boolean" truevalue="Yes" falsevalue="No" label="Perform permutation significance tests"/>
+        <param argument="--pvalue" type="float" value=".05" min="0" max="1" label="Significance value for site-tests"/>
+        <param argument="--qvalue" type="float" value=".2" min="0" max="1" label="Significance value for FDR reporting"/>
+        <param argument="--srv" type="select" label="Include synonymous rate variation" help = "Allow synonymous rates to vary from site to site">
+            <option value="Yes">Yes (recommended)</option>
+            <option value="No">No</option>
+        </param>
+    </inputs>
+    <outputs>
+        <data name="cfel_log" format="txt"/>
+        <data name="cfel_output" format="hyphy_results.json" />
+    </outputs>
+    <tests>
+        <test>
+            <param name="input_file" ftype="fasta" value="absrel-in1.fa"/>
+            <param name="input_nhx" ftype="nhx" value="absrel-in1.nhx"/>
+            <param name="branch_set" value="Internal branches,Terminal branches" />
+            <output name="cfel_output" file="cfel-out1.json" compare="sim_size"/>
+        </test>
+    </tests>
+    <help><![CDATA[
+Contrast-FEL : A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches
+==================================================================================================================
+
+A number of evolutionary hypotheses can be tested by comparing selective pressures among sets of branches in a phylogenetic tree. When the question of interest is to identify specific sites within genes that may be evolving differently, a common approach is to perform separate analyses on subsets of sequences and compare parameter estimates in a post hoc fashion. This approach is statistically suboptimal and not always applicable. Here, we develop a simple extension of a popular fixed effects likelihood method in the context of codon-based evolutionary phylogenetic maximum likelihood testing, Contrast-FEL. It is suitable for identifying individual alignment sites where any among the K≥2 sets of branches in a phylogenetic tree have detectably different ω ratios, indicative of different selective regimes. Using extensive simulations, we show that Contrast-FEL delivers good power, exceeding 90% for sufficiently large differences, while maintaining tight control over false positive rates, when the model is correctly specified. We conclude by applying Contrast-FEL to data from five previously published studies spanning a diverse range of organisms and focusing on different evolutionary questions.
+]]>
+
+    </help>
+
+    <expand macro="citations">
+        <citation type="doi">10.1093/molbev/msaa263</citation>
+    </expand>
+</tool>