Skip to content

Commit e76c97b

Browse files
authored
Merge pull request #11 from QinLiu-UTas/master
Add robust phylogenetics tutorial to Advanced Tutorial
2 parents fff86ce + 1e2c39e commit e76c97b

1 file changed

Lines changed: 26 additions & 0 deletions

File tree

doc/Advanced-Tutorial.md

Lines changed: 26 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -520,6 +520,32 @@ This will print an output file `example.phy.mlrate` that looks like:
520520
9 2.44313
521521
10 0.00001
522522

523+
524+
Robust phylogenetics analysis using trimmed log-likelihood method
525+
----------------------
526+
<div class="hline"></div>
527+
528+
Phylogenetic inference can be highly sensitive to fast-evolving, saturated or erroneous sites in a sequence alignment. To address this issue, IQ-TREE implements the `trimmed log-likelihood` method - a robust and dynamic approach that improves tree inference by selectively down-weighting problematic sites.
529+
530+
This method works by dynamically excluding a user-defined proportion of sites with the lowest log-likelihood values during the tree search. As the search progresses, the likelihood of each site is recalculated at each step using current tree and model parameters. This ensures that site removal is always conditional on the current model, tree topology and branch lengths, avoiding circularity.
531+
532+
To use the trimmed log-likelihood method, please make sure that IQ-TREE version 3.0 or later is installed. In the command-line interface, the method is invoked using the option `--robust-phy`. Although it is referred to here as the trimmed log-likelihood method, IQ-TREE uses the name `--robust-phy` to reflect the broader goal of improving the robustness of phylogenetic inference.
533+
534+
You can run the trimmed log-likelihood method from the command line by specifying the alignment, a substitution model, and the proportion of sites to retain:
535+
536+
iqtree3 -s <MY_ALIGNMENT> --robust-phy <PROPORTION_TO_RETAIN> -m <MODEL>
537+
538+
Additional options are available to assist downstream analysis. For instance, IQ-TREE can write site log-likelihoods to a `.sitelh` file, allowing users to identify the excluded sites by examining their log-likelihood values.
539+
540+
For example, for a dataset `data.phy`, if users apply a `JC` model and trim `2%` of sites (i.e., retain `98%` of sites) and wish to generate an output that includes the site log-likelihoods, the corresponding command would be:
541+
542+
iqtree3 -s <data.phy> --robust-phy <0.98> -m <JC> -wsl
543+
544+
If you use the trimmed log-likelihood method in a publication, please cite:
545+
546+
> __Liu, Qin, Bui Quang Minh, Robert Lanfear, Michael A. Charleston, Shane A. Richards, and Barbara R. Holland__ Robust Phylogenetics. _bioRxiv_ (2025): 2025-04.
547+
<https://doi.org/10.1101/2025.04.01.646540>
548+
523549
Where to go from here?
524550
----------------------
525551
<div class="hline"></div>

0 commit comments

Comments
 (0)