Select Model For Me with components
Phylogenetic inference with MrBayes. The the model of evolution is defined with the help of PartitionFinder, inference and test of MCMC convergence (GEOKS) and fit of the model on the data (Posterior Predictive test)
Preview
Run
Run this Workflow in the Taverna Workbench...
Option 1:
Copy and paste this link into File > 'Open workflow location...'
http://myexperiment.org/workflows/3404/download?version=21
[ More Info ]
Taverna is available from http://taverna.sourceforge.net/
If you are having problems downloading it in Taverna, you may need to provide your username and password in the URL so that Taverna can access the Workflow:
Replace http:// in the link above with http://yourusername:yourpassword@
Workflow Components
![]() | ![]() |
saverio vicario giacinto donvito bachir balech |
![]() | ![]() |
Select Model For Me with components |
![]() | ![]() |
Phylogenetic inference with MrBayes. The the model of evolution is defined with the help of PartitionFinder, inference and test of MCMC convergence (GEOKS) and fit of the model on the data (Posterior Predictive test) |
![]() | ![]() |
None
Inputs (8) | ![]() |
Name | Description |
---|---|
MSA_fasta | Multiple Sequence Alignment in Fasta format |
Possible_Partition | Partitioning of the MSA using nexus/raxml syntax. Define at least 2 partitions. Partition could define codon positions and/or gene boundaries and/or domains |
Statistical_Criterium | What statistical criterium should be used to compare different models? AIC = Aikake information Criterium AICc =Aikake information Criterium correct per sample size BIC = Bayesian Information Criterium 3 possible values: AIC, AICc, or BIC |
Number_of_runs | Number of independent runs to check convergence of the markovian integration necessary to perform bayesian phylogenetic inference. Minimum value 2. Always an integer |
Number_of_MCMCMC_generations | Number of Generation for each run of the markovian integration necessary to the phylogenetic inference. The convergence test (GeoKS) will tell post hoc if the number was sufficient. If is not the case try with a larger number following the educated guess of GeoKS |
File_name | Filename for nexus input file |
Linking_branches | Should models assume that branches in the tree are linked or unlinked across alignment parts? Two possible strings: linked unlinked |
SampleSize | Sample Size from the Posterior Distribution to test fit of model on Data. This integer should be smaller than number of generation divided by 100. |
Processors (14) | ![]() |
Name | Type | Description |
---|---|---|
SelectingEvolutionaryModelWithPartFinder | component | http://www.myexperiment.org/workflows/3943.html |
MrBayes | component |
ComponentMrBayes |
JSTWebDav3Input | component |
http://www.myexperiment.org/workflows/3942.html ComponentJSTWebDav3Input |
RandGen | localworker |
Scriptimport java.util.Random; Random r = new Random(); Integer a=r.nextInt(); if (a<0){ a=-1*a; } out=a ; |
ConVergenceTreeDiagnosticGeoKS | component |
http://www.myexperiment.org/workflows/3945.html ComponentConVergenceTreeDiagnosticGeoKS |
condition_checker | beanshell |
This service checks if the condition is "true" and if so outputs a value in the pass_flags list and
leaves the fail_flags empty. If the condition is not true then the service outputs a value in the
fail_flags list and leaves the pass_flags empty. Scriptimport java.util.ArrayList; ArrayList pass_flags = new ArrayList(); ArrayList fail_flags = new ArrayList(); if (condition.equals("True")) { pass_flags.add(condition); } else { fail_flags.add(condition); } |
Success_check | workflow | |
ConsensusTree | component |
http://www.myexperiment.org/workflows/3946.html ComponentConsensusTree |
TreeVisualizationItol | component |
http://www.myexperiment.org/workflows/3947.html ComponentTreeVisualizationItol |
MrBayesPPtest | component |
http://www.myexperiment.org/workflows/3952.html ComponentMrBayesPPtest |
NoSpace | beanshell |
Scriptout1=string.replaceAll(" ","_"); |
TargetOnly | beanshell |
Scriptout1=string.replaceAll(".+/",""); |
Failure_check | workflow | |
String_List_Union | localworker |
ScriptSet results = new HashSet(); results.addAll(list1); results.addAll(list2); List union = new ArrayList(); union.addAll(results); |
Beanshells (3) | ![]() |
Name | Description | Inputs | Outputs |
---|---|---|---|
condition_checker | This service checks if the condition is "true" and if so outputs a value in the pass_flags list and leaves the fail_flags empty. If the condition is not true then the service outputs a value in the fail_flags list and leaves the pass_flags empty. | condition |
pass_flags fail_flags |
NoSpace | string | out1 | |
TargetOnly | string | out1 |
Outputs (11) | ![]() |
Name | Description |
---|---|
Plot | Result of the Posterior Predictive test. Histogram of the distribution of the complexities (log of sum of sites entropy or maximal possible loglikelihood score) of the simulated data using the posterior distribution parameters compared with observed data complexity. The 1-alpha high posterior density of the distribution show the region where simulate data complexity match the observed one. Larger observed complexity indicates model too simplistic, while the contrary indicates overparametrization of the model. |
newickTree | small xml with tag res with one or more tag tree each one with a body that contain a tree in netwick format. each tree represent the consensus for a given partition of group of partition in the same order that are cited in the nexus input file of mrbayes |
Geoks | Result of the GeoKS test of convergence. XML format |
PhyloInferenceOutput | The output of the consensus service is a path were to obtain a zipped folder that includes all output from phylogenetic inference and the one of the consensus |
ViewTree | Link to visualize consensus Tree of phylogenetic inference. It could be more than one link if the model assume more than one tree (combination topology + branch length set). Use of the figure in the link below should be always accompagned by appropriate citation: Letunic I and Bork P (2011) Nucleic Acids Res doi: 10.1093/nar/gkr201 Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy |
MrBayesOutput | original output of mrbayes |
DetailsPPtest | Numeric Result of the Posterior Predictive test. |
PartitionFinderOutput | Main output of PartitionFinder, where the preffered partitioned model is described |
PartitionFinderDetails | Path were to retrieve all the detials of the partitionfinder software |
GeoKSDetails | Path were to get all details of convergence test calculation. |
log |
Datalinks (41) | ![]() |
Source | Sink |
---|---|
Statistical_Criterium | SelectingEvolutionaryModelWithPartFinder:StatCriterium |
Possible_Partition | SelectingEvolutionaryModelWithPartFinder:PossiblePartition |
MSA_fasta | SelectingEvolutionaryModelWithPartFinder:MSA_fasta |
Number_of_runs | SelectingEvolutionaryModelWithPartFinder:NumberRuns |
Number_of_MCMCMC_generations | SelectingEvolutionaryModelWithPartFinder:MCMCMCGenerationNumber |
Linking_branches | SelectingEvolutionaryModelWithPartFinder:LinkingBranches |
RandGen:out | SelectingEvolutionaryModelWithPartFinder:Token |
JSTWebDav3Input:REST_Service_redirection | MrBayes:Path2MrBayesInputFile |
SelectingEvolutionaryModelWithPartFinder:MrBayesInput | JSTWebDav3Input:InputBody |
RandGen:out | JSTWebDav3Input:Token |
NoSpace:out1 | JSTWebDav3Input:File_name |
MrBayes:PathOfResultsFolder | ConVergenceTreeDiagnosticGeoKS:MrBayesOuputPath |
ConVergenceTreeDiagnosticGeoKS:Success | condition_checker:condition |
condition_checker:pass_flags | Success_check:Sentinel |
MrBayes:PathOfResultsFolder | ConsensusTree:MrBayesFolder |
ConVergenceTreeDiagnosticGeoKS:Burnin | ConsensusTree:BurnIn |
ConsensusTree:Trees | TreeVisualizationItol:treefile |
ConVergenceTreeDiagnosticGeoKS:Burnin | MrBayesPPtest:Burnin |
SampleSize | MrBayesPPtest:SampleSize |
MrBayes:PathOfResultsFolder | MrBayesPPtest:MrBayesOutputFolder |
TargetOnly:out1 | MrBayesPPtest:NameRun |
File_name | NoSpace:string |
JSTWebDav3Input:REST_Service_redirection | TargetOnly:string |
condition_checker:fail_flags | Failure_check:Sentinel |
TreeVisualizationItol:Log | String_List_Union:list2 |
MrBayesPPtest:Log | String_List_Union:list1 |
ConsensusTree:Log | String_List_Union:list1 |
ConVergenceTreeDiagnosticGeoKS:Log | String_List_Union:list1 |
SelectingEvolutionaryModelWithPartFinder:Log | String_List_Union:list1 |
MrBayes:Log | String_List_Union:list1 |
MrBayesPPtest:Plot | Plot |
ConsensusTree:Trees | newickTree |
ConVergenceTreeDiagnosticGeoKS:NiceXML | Geoks |
ConsensusTree:ResultsFolder | PhyloInferenceOutput |
TreeVisualizationItol:html | ViewTree |
MrBayes:PathOfResultsFolder | MrBayesOutput |
MrBayesPPtest:PPtestDetails | DetailsPPtest |
SelectingEvolutionaryModelWithPartFinder:PartitionFinderOuput | PartitionFinderOutput |
SelectingEvolutionaryModelWithPartFinder:Path | PartitionFinderDetails |
ConVergenceTreeDiagnosticGeoKS:OutputPath | GeoKSDetails |
String_List_Union:union | log |
Coordinations (2) | ![]() |
Controller | Target |
---|---|
Success_check | MrBayesPPtest |
Success_check | ConsensusTree |
Workflow Type
Version 21 (latest)
(of 21)
None
Shared with Groups (0)
None
Log in to add to one of your Packs
Statistics
In chronological order:
-
Created by Saverio Vicario on Wednesday 20 February 2013 12:55:56 (UTC)
-
Created by Saverio Vicario on Wednesday 20 February 2013 12:57:51 (UTC)
-
Created by Saverio Vicario on Wednesday 20 February 2013 14:40:18 (UTC)
-
Created by Saverio Vicario on Monday 25 March 2013 16:29:14 (UTC)
Revision comment:R script do not ask any more boa library
-
Created by Saverio Vicario on Wednesday 27 March 2013 14:55:28 (UTC)
Revision comment:patched python script
-
Created by Saverio Vicario on Thursday 28 March 2013 10:04:56 (UTC)
Revision comment:adjusted python scripts such to handle non unix newline endindgs and to repair a bug in the model parsing of partitionfinder (recognize model without any parameters for change across sites)
-
Created by Saverio Vicario on Thursday 28 March 2013 15:43:09 (UTC)
Revision comment:other improvment on python script
-
Created by Saverio Vicario on Thursday 16 May 2013 13:16:49 (UTC)
-
Created by Saverio Vicario on Friday 17 May 2013 14:11:56 (UTC)
Revision comment:moved all but one interaction page to input, and made the parser of the partitiondefinition more robust to lack of newline where they should be
-
Created by Saverio Vicario on Tuesday 04 June 2013 12:41:24 (UTC)
-
Created by Saverio Vicario on Wednesday 05 June 2013 09:59:57 (UTC)
Revision comment:new address of JST
-
Created by Saverio Vicario on Thursday 06 June 2013 09:11:26 (UTC)
Revision comment:corrected JST interaction, add output to mrbayes in case of failure of convergence
-
Created by Saverio Vicario on Thursday 06 June 2013 13:16:01 (UTC)
-
Created by Saverio Vicario on Thursday 06 June 2013 14:31:10 (UTC)
Revision comment:adjusting output port name
-
Created by Saverio Vicario on Thursday 13 June 2013 09:37:23 (UTC)
Revision comment:add a new ouput port to access details posterior predictive test
-
Created by Saverio Vicario on Thursday 20 June 2013 14:53:44 (UTC)
Revision comment:Added log function
-
Created by Saverio Vicario on Thursday 20 June 2013 15:49:35 (UTC)
Revision comment:correct citation of itol
-
Created by Saverio Vicario on Friday 16 August 2013 12:09:45 (UTC)
-
Created by Saverio Vicario on Thursday 05 December 2013 17:06:18 (UTC)
Revision comment:sardine and retry on get of the webservices
-
Created by Saverio Vicario on Thursday 05 December 2013 17:10:13 (UTC)
Revision comment:retry and sardine
-
Created by Saverio Vicario on Friday 31 January 2014 11:28:33 (UTC)
Revision comment:All service are within component
Reviews
(0)
Other workflows that use similar services
(0)
There are no workflows in myExperiment that use similar services to this Workflow.
No comments yet
Log in to make a comment