From ad2079d0e164f2c4be4884e40e3be81ac16505bf Mon Sep 17 00:00:00 2001
From: RevBayes analysis <1695515+ms609@users.noreply.github.com>
Date: Mon, 17 Feb 2025 10:55:09 +0000
Subject: [PATCH] Check multifurcating trees

---
 tests/testthat/test-Concordance.R | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/tests/testthat/test-Concordance.R b/tests/testthat/test-Concordance.R
index 0030176aa..9e0b37e92 100644
--- a/tests/testthat/test-Concordance.R
+++ b/tests/testthat/test-Concordance.R
@@ -76,6 +76,12 @@ test_that("QuartetConcordance(method = minhq)", {
   expect_concordance(iq = "iq", 9, c(0, 0, 1 / 12, 0, 0))
   expect_equal(unname(QuartetConcordance(tree, dat, method = "iqtree")), 
                c(56.7, 0, 85.4, 0, 62.5) / 100, tolerance = 0.01)
+  
+  collapsed <- CollapseNode(tree, c(12, 13))
+  expect_equal(
+    QuartetConcordance(collapsed, dat, method = "iqtree"), 
+    QuartetConcordance(tree, dat, method = "iqtree")[-c(2, 3)]
+  )
 })
 
 test_that("QuartetConcordance(method = minh)", {
@@ -114,6 +120,11 @@ test_that("QuartetConcordance(method = minh)", {
     unname(QuartetConcordance(tree, dat, method = "minh", n = 1234)),
     c(56.7, 0, 85.4, 0, 62.5) / 100)
   
+  collapsed <- CollapseNode(tree, c(12, 13))
+  expect_equal(tolerance = 0.05,
+    QuartetConcordance(collapsed, dat, method = "minh", n = 1234),
+    QuartetConcordance(tree, dat, method = "minh", n = 1234)[-(2:3)]
+  )
 })
 
 test_that("QuartetConcordance() calculates correct values - weighting", {