comparison hyphy_cfel.xml @ 0:3d4ea0544af1 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/hyphy/ commit 2742ee3b4e90f65352845265d2f85c4263e0eabb"
author iuc
date Tue, 20 Apr 2021 10:26:31 +0000
parents
children fb833dd994f2
comparison
equal deleted inserted replaced
-1:000000000000 0:3d4ea0544af1
1 <?xml version="1.0"?>
2 <tool id="hyphy_cfel" name="HyPhy-CFEL" version="@VERSION@+galaxy0" profile="19.09">
3 <description>Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches</description>
4 <macros>
5 <import>macros.xml</import>
6 </macros>
7 <expand macro="requirements"/>
8 <command detect_errors="exit_code"><![CDATA[
9 ln -s '$input_file' cfel_input.fa &&
10 ln -s '$input_nhx' cfel_input.nhx &&
11 hyphy contrast-fel
12 --code '$gencodeid'
13 #for branch in str($branch_set).split(','):
14 --branch-set '$branch'
15 #end for
16 --srv '$srv'
17 --permutations '$permutations'
18 --pvalue '$pvalue'
19 --qvalue '$qvalue'
20 --output '$cfel_output'
21 cfel_input.fa
22 cfel_input.nhx
23 > '$cfel_log'
24 @ERRORS@
25 ]]></command>
26 <inputs>
27 <expand macro="inputs"/>
28 <expand macro="gencode"/>
29 <param argument="--branch-set" type="select" multiple="true" min="2" max="2" label="Set of branches to test">
30 <option value="Internal branches">Internal branches</option>
31 <option value="Terminal branches">Terminal branches</option>
32 <option value="Random set of branches">Random set of branches</option>
33 <option value="Unlabeled branches">Unlabeled branches</option>
34 </param>
35 <param argument="--permutations" type="boolean" truevalue="Yes" falsevalue="No" label="Perform permutation significance tests"/>
36 <param argument="--pvalue" type="float" value=".05" min="0" max="1" label="Significance value for site-tests"/>
37 <param argument="--qvalue" type="float" value=".2" min="0" max="1" label="Significance value for FDR reporting"/>
38 <param argument="--srv" type="select" label="Include synonymous rate variation" help = "Allow synonymous rates to vary from site to site">
39 <option value="Yes">Yes (recommended)</option>
40 <option value="No">No</option>
41 </param>
42 </inputs>
43 <outputs>
44 <data name="cfel_log" format="txt"/>
45 <data name="cfel_output" format="hyphy_results.json" />
46 </outputs>
47 <tests>
48 <test>
49 <param name="input_file" ftype="fasta" value="absrel-in1.fa"/>
50 <param name="input_nhx" ftype="nhx" value="absrel-in1.nhx"/>
51 <param name="branch_set" value="Internal branches,Terminal branches" />
52 <output name="cfel_output" file="cfel-out1.json" compare="sim_size"/>
53 </test>
54 </tests>
55 <help><![CDATA[
56 Contrast-FEL : A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches
57 ==================================================================================================================
58
59 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.
60 ]]>
61
62 </help>
63
64 <expand macro="citations">
65 <citation type="doi">10.1093/molbev/msaa263</citation>
66 </expand>
67 </tool>