X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=wiki%2FRIO.wiki;h=107decaf8baa19339a0162780e9836b82cc716a4;hb=9fa5e015e26fb376a032ea1b2c758bd2802af021;hp=38a8660e43c0f5093938acab75e0126e28bf4a38;hpb=162f375cc97383cd18534716e81eeb724db28307;p=jalview.git diff --git a/wiki/RIO.wiki b/wiki/RIO.wiki index 38a8660..107deca 100644 --- a/wiki/RIO.wiki +++ b/wiki/RIO.wiki @@ -6,49 +6,35 @@ RIO (Resampled Inference of Orthologs) is a method for automated phylogenomics based on explicit phylogenetic inference. RIO analyses are performed over resampled phylogenetic trees to estimate the reliability of orthology assignments. + == Usage == {{{ -java -Xmx2048m -cp forester.jar org.forester.application.rio [options] [outfile] +java -Xmx2048m -cp forester.jar org.forester.application.rio [options] [logfile] }}} -=== Options === - - * -co: cutoff for ortholog output (default: 50) - - * -t : file-name for output table - - * -q : name for query (sequence/node), if this is bing used, [outfile] is required as well - - * -s : sort (default: 2) - - * -u : to output ultra-paralogs (species specific expansions/paralogs) - * -cu: cutoff for ultra-paralog output (default: 50) - -==== Sort ==== +=== Options === - * 0: orthologies - * 1: orthologies > super orthologies - * 2: super orthologies > orthologies + * -b : to use SDIR instead of GSDIR (faster, but non-binary species trees are disallowed) + ==== Gene trees ==== -The gene trees ideally are in [http://www.biomedcentral.com/1471-2105/10/356/ phyloXML] format, with taxonomy and sequence data in appropriate fields; but can also be in New Hamphshire (Newick) or Nexus format, as long as species information can be extracted from the gene names (e.g. "HUMAN" from "BCL2_HUMAN"). +The gene trees ideally are in [http://www.biomedcentral.com/1471-2105/10/356/ phyloXML] format, with taxonomy and sequence data in appropriate fields; but can also be in New Hamphshire (Newick) or Nexus format, as long as species information can be extracted from the gene names (e.g. "HUMAN" from "BCL2_HUMAN") ([http://forester.googlecode.com/files/gene_trees_rio.nh example]). +All gene trees must be completely binary. -==== Species tree ==== -Must be in [http://www.biomedcentral.com/1471-2105/10/356/ phyloXML] format ([http://forester.googlecode.com/files/species.xml example]). +==== Species tree ==== +Must be in [http://www.biomedcentral.com/1471-2105/10/356/ phyloXML] format ([http://forester.googlecode.com/files/species_tree_rio.xml example]). +The species tree is allowed to have nodes with more than two descendents (polytomies) as long as GSDIR ([GSDI GSDI] re-rooting) algorithm is used. -=== Examples === -`rio gene_trees.nh species.xml outfile -q=BCL2_HUMAN -t=outtable -u -cu=60 -co=60` -`rio gene_trees.nh species.xml -t=outtable` +=== Example === +`rio gene_trees.nh species.xml outtable.tsv log.txt` - === Example files === - * [http://forester.googlecode.com/files/wnt_gene_tree.xml gene tree] - * [http://forester.googlecode.com/files/species.xml species tree] - * [http://forester.googlecode.com/files/wnt_gsdi_log.txt log file (output)] - + * [http://forester.googlecode.com/files/gene_trees_rio.nh gene trees file] + * [http://forester.googlecode.com/files/species_tree_rio.xml species tree file] + == References ==