Skip to content

Commit 7c6f829

Browse files
committed
Updated Estimating gene, site, and quartet concordance vectors (markdown)
1 parent 795c9bf commit 7c6f829

1 file changed

Lines changed: 12 additions & 1 deletion

File tree

doc/Estimating-gene,-site,-and-quartet-concordance-vectors.md

Lines changed: 12 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -327,4 +327,15 @@ And the CSV looks like this:
327327
| quartet| 4 | 0 | 0 | 0 |
328328

329329
You'll notice that both include 95% confidence intervals for the concordance and discordance factors. These are calculated using 1000 bootstraps of the count data, and provide useful context for interpreting the values, and particularly for interpreting potential *differences* in the values.
330-
330+
331+
You'll notice that for node 642, there is an enormous amount of discordance. Indeed, the first three entries of the vector all have confidence intervals that overlap considerably for gene and quartet concordance factors. So although ASTRAL did what it is supposed to do and chose the node with the highest quartet concordance factor, it would be difficult to be extremely confident that this is the correct topology for the species tree. In support of that, examining the `concordance_vectors.csv` file shows that the posterior probability for this branch is just 0.5, which is extremely low.
332+
333+
In contrast, node 641 has a lot less discordance:
334+
335+
![concordance_table_641](https://github.com/iqtree/iqtree2/assets/895251/e2968f22-9201-433e-94ea-c0d5c7095134)
336+
337+
Almost all genes and quartets are concordance with this node. The sites have more discordance, and &#936;<sub>2</sub> seems a lot higher for sites than the other entries in the vector. This deserves further investigation, but *could* occur if a lot of the genes that support &#936;<sub>2</sub> are very informative, and/or if there's lot of homoplasy.
338+
339+
# Conclusion
340+
341+
We hope this recipe provides some useful guidance on calculating concordance vectors for your data!

0 commit comments

Comments
 (0)