11---
22layout : userdoc
33title : " Assessing Phylogenetic Assumptions"
4- author : M Bui
5- date : 2021-03-11
4+ author : Minh Bui
5+ date : 2021-03-15
66docid : 5
77icon : info-circle
88doctype : tutorial
99tags :
1010- tutorial
1111description : This guide is about evaluating the suitability of the data for phylogenetic analysis.
12+ sections :
13+ - name : Tests of symmetry
14+ url : tests-of-symmetry
15+ - name : Likelihood mapping
16+ url : likelihood-mapping
17+
1218---
1319
20+
21+ Assessing phylogenetic assumptions
22+ ==================================
23+
1424It is important to know that phylogenetic models rely on various simplifying assumptions to
1525ease computations. If your data severely violate these assumptions, it might
1626cause bias in phylogenetic estimates of tree topologies and other model
@@ -22,78 +32,11 @@ likely in both directions), and _homogeneity_ (substitution rates remain constan
2232This document shows several ways to check some of these assumptions that you
2333should perform before doing phylogenetic analysis.
2434
25- Likelihood mapping analysis
26- ---------------------------
27- <div class =" hline " ></div >
28-
29- Likelihood mapping ([ Strimmer and von Haeseler, 1997] ) is a visualisation method
30- to display the phylogenetic information of an alignment. It visualises the _ treelikeness_
31- of all quartets in a single triangular graph and therefore renders a quick
32- interpretation of the phylogenetic content.
33-
34- A simple likelihood mapping analysis can be conducted with:
35-
36- iqtree -s example.phy -lmap 2000 -n 0
37-
38- where ` -lmap ` option specify the number of quartets of taxa that will be drawn randomly
39- from the alignment. ` -n 0 ` tells IQ-TREE to stop the analysis right after running the
40- likelihood mapping. IQ-TREE will print the result in the ` .iqtree ` report file as well
41- as the likelihood mapping plot ` .lmap.svg ` (in SVG format) and ` .lmap.eps ` file (in EPS
42- figure format).
43-
44- You can now view the likelihood mapping plot file ` example.phy.lmap.svg ` , which looks like this:
45-
46- ![ Likelihood mapping plot.] ( images/example.phy.lmap.pdf )
47-
48- It shows phylogenetic information of the alignment ` example.phy ` .
49-
50- * Top sub-figure: distribution of quartets depicted by dots on the likelihood mapping plot.
51- * Left sub-figure: percentages of quartets falling in each of the three areas. The
52- three areas show support for one of the different groupings like (a,b)-(c,d).
53- * Right sub-figure: percentages of quartets falling in each of the seven areas.
54- Quartets falling into the three corners are informative and called fully-resolved quartets.
55- Those in three rectangles are partly informative (partly resolved quartets) and those in the center are uninformative
56- (unresolved quartets). A good data set should have high number of fully resolved quartets
57- and low number of unresolved quartets.
58-
59- The meanings can also be found in the ` LIKELIHOOD MAPPING STATISTICS ` section of the report file ` example.phy.iqtree ` :
60-
61-
62- LIKELIHOOD MAPPING STATISTICS
63- -----------------------------
64-
65- (a,b)-(c,d) (a,b)-(c,d)
66- /\ /\
67- / \ / \
68- / \ / 1 \
69- / a1 \ / \ / \
70- /\ /\ / \/ \
71- / \ / \ / /\ \
72- / \ / \ / 6 / \ 4 \
73- / \/ \ /\ / 7 \ /\
74- / | \ / \ /______\ / \
75- / a3 | a2 \ / 3 | 5 | 2 \
76- /__________|_________\ /_____|________|_____\
77- (a,d)-(b,c) (a,c)-(b,d) (a,d)-(b,c) (a,c)-(b,d)
78-
79- Division of the likelihood mapping plots into 3 or 7 areas.
80- On the left the areas show support for one of the different groupings
81- like (a,b|c,d).
82- On the right the right quartets falling into the areas 1, 2 and 3 are
83- informative. Those in the rectangles 4, 5 and 6 are partly informative
84- and those in the center (7) are not informative.
85- .....
86-
87-
88- The [ command reference] ( Command-Reference#likelihood-mapping-analysis ) will provide
89- more options and how to perform 2-, 3-, or 4-cluster likelihood mapping analysis.
90-
91-
9235Tests of symmetry
9336-----------------
9437
9538IQ-TREE provides three matched-pairs tests of symmetry ([ Naser-Khdour et al., 2019] ) to
96- test the three assumptions of stationarity, reversibility and homogeneity (SRH).
39+ test the two assumptions of stationarity and homogeneity (SRH).
9740A simple analysis:
9841
9942 iqtree2 -s example.phy -p example.nex --symtest-only
@@ -142,7 +85,7 @@ Now you may want to perform the phylogenetic analysis excluding all "bad" partit
14285
14386 iqtree2 -s example.phy -p example.nex --symtest-remove-bad
14487
145- that will remove all "bad" partitions with SymPval < 0.05 and continue the analysis with the
88+ that will remove all "bad" partitions where SymPval < 0.05 and continue the analysis with the
14689remaining "good" partitions. You may then compare the trees from "all" partitions
14790and from "good" only partitions to see if there is significant difference between them
14891with [ tree topology tests] ( Advanced-Tutorial#tree-topology-tests ) .
@@ -161,6 +104,73 @@ TEST OF SYMMETRY:
161104```
162105
163106
107+ Likelihood mapping
108+ ------------------
109+ <div class =" hline " ></div >
110+
111+ Likelihood mapping ([ Strimmer and von Haeseler, 1997] ) is a visualisation method
112+ to display the phylogenetic information of an alignment. It visualises the _ treelikeness_
113+ of all quartets in a single triangular graph and therefore renders a quick
114+ interpretation of the phylogenetic content.
115+
116+ A simple likelihood mapping analysis can be conducted with:
117+
118+ iqtree -s example.phy -lmap 2000 -n 0
119+
120+ where ` -lmap ` option specify the number of quartets of taxa that will be drawn randomly
121+ from the alignment. ` -n 0 ` tells IQ-TREE to stop the analysis right after running the
122+ likelihood mapping. IQ-TREE will print the result in the ` .iqtree ` report file as well
123+ as the likelihood mapping plot ` .lmap.svg ` (in SVG format) and ` .lmap.eps ` file (in EPS
124+ figure format).
125+
126+ You can now view the likelihood mapping plot file ` example.phy.lmap.svg ` , which looks like this:
127+
128+ ![ Likelihood mapping plot.] ( images/example.phy.lmap.pdf )
129+
130+ It shows phylogenetic information of the alignment ` example.phy ` .
131+
132+ * Top sub-figure: distribution of quartets depicted by dots on the likelihood mapping plot.
133+ * Left sub-figure: percentages of quartets falling in each of the three areas. The
134+ three areas show support for one of the different groupings like (a,b)-(c,d).
135+ * Right sub-figure: percentages of quartets falling in each of the seven areas.
136+ Quartets falling into the three corners are informative and called fully-resolved quartets.
137+ Those in three rectangles are partly informative (partly resolved quartets) and those in the center are uninformative
138+ (unresolved quartets). A good data set should have high number of fully resolved quartets
139+ and low number of unresolved quartets.
140+
141+ The meanings can also be found in the ` LIKELIHOOD MAPPING STATISTICS ` section of the report file ` example.phy.iqtree ` :
142+
143+
144+ LIKELIHOOD MAPPING STATISTICS
145+ -----------------------------
146+
147+ (a,b)-(c,d) (a,b)-(c,d)
148+ /\ /\
149+ / \ / \
150+ / \ / 1 \
151+ / a1 \ / \ / \
152+ /\ /\ / \/ \
153+ / \ / \ / /\ \
154+ / \ / \ / 6 / \ 4 \
155+ / \/ \ /\ / 7 \ /\
156+ / | \ / \ /______\ / \
157+ / a3 | a2 \ / 3 | 5 | 2 \
158+ /__________|_________\ /_____|________|_____\
159+ (a,d)-(b,c) (a,c)-(b,d) (a,d)-(b,c) (a,c)-(b,d)
160+
161+ Division of the likelihood mapping plots into 3 or 7 areas.
162+ On the left the areas show support for one of the different groupings
163+ like (a,b|c,d).
164+ On the right the right quartets falling into the areas 1, 2 and 3 are
165+ informative. Those in the rectangles 4, 5 and 6 are partly informative
166+ and those in the center (7) are not informative.
167+ .....
168+
169+
170+ The [ command reference] ( Command-Reference#likelihood-mapping-analysis ) will provide
171+ more options and how to perform 2-, 3-, or 4-cluster likelihood mapping analysis.
172+
173+
164174[ Strimmer and von Haeseler, 1997 ] : http://www.pnas.org/content/94/13/6815.long
165175[ Naser-Khdour et al., 2019 ] : https://doi.org/10.1093/gbe/evz193
166176
0 commit comments