Phylogenetic Inference and Congruence

Biological trees agree above chance, but designed systems can too. Common descent gains force from genomic insertions, chromosome order, and fossils, not congruence alone.

Draft:  August 2026

"Bootstrapping provides us with a confidence interval within which is contained not the true phylogeny, but the phylogeny that would be estimated on repeated sampling of many characters from the underlying pool of characters." — Joe Felsenstein, in the 1985 paper that brought the bootstrap into phylogenetics1

"the original tree of life concept is obsolete: it would not even be a 'tree of one percent'.  Nevertheless, there seems to be a discernible signal of consistency between the trees in the forest of life, down to the deepest branching levels." — Puigbò, Wolf and Koonin, 20092

"It would be easy to attribute their similarities to a common designer, but such an attribution would be trivial. ... what distinguishes a common designer from a common ancestor?" — Todd Wood, creationist biologist, asking this article's question in his own paper3

Summary

The strongest argument for common descent is not that living things are similar.  It is that they are similar in a particular pattern, and that independent sets of evidence keep recovering the same pattern.  Build a tree from one gene.  Build another from a different gene.  They agree far more than chance allows.  Do it with bones instead of DNA and you get much the same tree again.

This article asks how that argument holds up.  It covers how trees are built, what the numbers printed on their branches mean, how much gene trees really conflict, and one question nobody had answered:  do designed objects produce a nested hierarchy too?

The short answer is that they do — and that this matters less than either side expects.

Every number below marked "measured here" was computed for this article.  The scripts are in the research/sim/ folder beside it, they take a seed, and they run in seconds.

The conclusions, in order of confidence

  1. Two random trees share almost nothing, so agreement between trees really is strong evidence.  There are 2.2 × 1020 possible trees for 20 species.  Two trees drawn at random from that set share, on average, one-sixth of one grouping.  Not one grouping in six.  One-sixth of a single grouping, out of the seventeen a 20-species tree has.  Add species and it falls further, settling near one-eighth.4  That is the baseline.  It is brutally low, and it is why congruence carries weight.

  2. A support value of 95% does not mean 95% confidence, and at genome scale it means almost nothing.  Felsenstein said in 1985 what the bootstrap measures, and it is not the probability that a branch is real.1  Sudhir Kumar and colleagues explain why it fails on large data:  "For genome-scale data sets, the site-sampling variance reduces to virtually zero."5  Every branch goes to 100% whether it is right or not.  In one documented case, parsimony backed a known-wrong tree in 96% of bootstrap replicates.6  In another, hominin-style cladistics returned seven clades above the 70% threshold and every one of them was false, at 98%, 97%, 95%, 92%, 78%, 73% and 73%.7

  3. Gene trees conflict a lot, and the papers reporting it also report the correction.  Salichos and Rokas built 1,070 yeast gene trees and got 1,070 different answers.8  Their own figure gives the numbers that matter:  average distance from gene trees to the species tree 0.40, among the gene trees themselves 0.52, among random trees 0.99.8  Gene trees disagree.  They are also nowhere near random, and the gap between 0.52 and 0.99 is the whole argument.  Filter for supported conflict and the numbers drop by an order of magnitude — a median of 6% across the ninety-one land-plant branches Wickett and colleagues label, measured here from their own figures.9

  4. Designed objects can produce strong congruence.  This is measured here, and it goes against the popular form of the argument.  Linux distributions were compared by which software packages they ship.  Two independent halves of the package list produce trees agreeing at QS = 0.91 on the 24 distributions whose real family history is documented, and 0.84 on 24 drawn from the full set — on a scale where 0 is chance and 1 is identical.10  That beats several real biological datasets, including a mammal alignment at 0.60 measured on the same 24 species.

  5. But that congruence is not the family history — and here it can be checked.  Linux distributions have a written-down genealogy.  Somebody recorded that Ubuntu was forked from Debian.  The package-derived tree was compared against those documented forks.  It matches at QS = 0.15, against a random tree's 95th percentile of 0.11.10  Across eight subsamples it ranged 0.08 to 0.23 and fell below the random 95th percentile twice out of eight.  The tree is highly reproducible and it is almost entirely the wrong tree.  Package choice tracks what a distribution is for, not what it descended from.

  6. So congruence by itself proves less than the standard argument claims.  That point goes to the creationists, and it should be conceded plainly.  Reproducible structure is not the same thing as genealogy.  A system built by choosing modules to meet requirements can produce a stable, strongly congruent tree that is not its history.

  7. The design side gets no comfort either, because the two patterns are distinguishable.  Designed objects show high self-agreement and low genealogical accuracy.  Copying traditions show high self-agreement and recover the real history.  Hand-copied manuscripts — 17 scribes copying under observation, true stemma recorded — score 0.88 for self-agreement, and the tree built from them recovers the recorded copying history well above chance:  against the released stemma the neighbour-joining tree scores QS = 0.49 and the parsimony tree 0.85, where a random tree's 95th percentile is 0.07.  What it does not recover is the direction of copying — who copied whom — which an undirected tree cannot carry.11  Genomes behave like the copying case, not the design case, and the reason is that their tree is confirmed by evidence that has nothing to do with sequence similarity:  shared retrotransposon insertions,12 chromosome-scale gene linkage,13 and the order of fossils in rock.14

  8. Congruence keeps climbing with more data for genealogies, which is the test that survives every objection about dataset size.  Measured here at a fixed 24 taxa:  yeast goes 0.22 → 0.43 → 0.66 → 0.91 → 1.00 as characters increase from 20 to 40 to 160 to 1,280 to all 445,152.15  Manuscripts climb the same way.  So do Linux distributions — because they are forked, and forking is copying.  Small product feature tables never climb, because they run out of characters at a few dozen.

  9. The published claim that artifacts fit trees as well as organisms do is real but does not survive inspection.  Collard, Shennan and Tehrani compared 21 cultural datasets against 21 biological ones:  mean retention index 0.59 against 0.61, p = 0.900.16  But their biological set pools two different populations, their cultural datasets average 94 characters against 1,228 for the biological ones, and their stated reason for choosing that statistic — that it does not depend on dataset size — fails on their own data.  Nunn and colleagues then showed the statistic cannot detect borrowing:  across 5,000 simulations, the amount of horizontal transmission did not enter the regression tree for either index.  They add the other half, which cuts the other way — an RI above about 0.60 is, in their words, "usually indicative of a high degree of vertical[] transmission (phylogenesis) and a low degree of horizontal transmission," and Collard's cultural mean is 0.59.17

  10. The one clean natural experiment says functional design does not track descent.  Gray, Bryant and Greenhill took 134 Polynesian canoe traits and tested them against the independently known language tree for the same eleven societies.  Basic vocabulary scored −2.89.  Symbolic canoe traits −6.13.  Randomised canoe data −6.30.  Functional canoe traits −6.64.18  Functional design fit the true population tree worse than randomised data.

  11. The only serious design model fitted to real genomes wins on its own terms, and its author says why that is not enough.  Winston Ewert's dependency-graph model beats a tree model by enormous margins on nine genomic databases — reproduced here from his released files, and every published figure reproduces in direction.19  But his tree model excludes, in his own words, lineage sorting, gene flow, horizontal transfer, convergence and gene resurrection.20  Expressed as a share of the tree model's own error, the graph explains 0.9% to 14.6% more of the data — and his own control, where a dependency graph is literally the true answer, scores 5.6%, above four of the nine real databases.  Nobody has replied in print in eight years.

  12. No creationist model makes a risky prediction about congruence that differs from common descent.  Baraminology accepts nested hierarchy inside a created kind, and kinds are large — one holobaramin of about 10,000 grass species, roughly 150 fossil horse species in one monobaramin.21  Everything below about the family level is common ground.  The movement's own best statistical critique is internal:  Colin Reeves, published by Answers in Genesis, concludes the core statistic "commits a category error" and "really should be abandoned."22

  13. Several quotations in wide circulation are quote-mines, including in this project's own notes.  Ford Doolittle's "It is as if we have failed at the task that Darwin set for us" is followed immediately by "But in fact, our science is working just as it should."23  Full audit in quotes; the ones that could not be verified are in quarantine.

The verdict.  The nested hierarchy is real, it is strong, and it is evidence of copying with modification.  It is not, by itself, evidence that the copying was biological rather than intelligent — a designer who works by modifying previous versions produces the same pattern, and Linux proves it.  What the argument actually rests on is narrower and better:  the tree recovered from DNA is independently confirmed by rare genomic events, chromosome structure and fossil order, and no design model has yet explained that convergence of evidence.  The popular version of the argument overreaches.  The narrow version stands.

How a tree gets built

Everything here depends on one procedure, so it is worth being concrete.

You start with a table.  Rows are organisms.  Columns are characters — a position in a DNA sequence, or the presence of a bone.  Each cell holds that organism's state for that character.  You then ask which branching diagram best explains the table.

There are too many diagrams to check.  The number of unrooted trees connecting N things is (2N−5)!!, and it explodes:

Number of taxa Possible unrooted trees
10 2,027,025
11 34,459,425
20 2.216 × 1020
30 8.687 × 1036
50 2.838 × 1074
100 1.700 × 10182

So every method is a scoring rule plus a search.  Four families are in use.

Distance methods reduce each pair of organisms to one number, then fit a tree to those numbers.  Neighbour-joining is the standard.24  It is what this article's software uses.  Its predecessor UPGMA assumes every lineage changes at the same speed.  When that fails, it fails hard.  In Saitou and Nei's own Table 6, with rates varying between lineages, UPGMA recovered the true 8-species tree in 0 of 100 trials.  Neighbour-joining got 75 and 72 out of 100 on the same data.24

Parsimony picks the tree needing the fewest changes.  It is intuitive.  It is also provably wrong under conditions described below.

Maximum likelihood writes down a model of how sequences change, then picks the tree making the data most probable.25  This is the workhorse.

Bayesian methods use the same models but sample trees in proportion to their probability.  Every branch gets a probability as a by-product.

There is a fifth idea that matters for what follows.  Genes really can have different histories from the species carrying them.  Species-tree methods treat each gene tree as noisy evidence about the species tree rather than as a measurement of it.26 27  This matters because gluing all the genes into one long sequence can be statistically inconsistent:  more data then drives you more confidently to the wrong answer.  Kubatko and Degnan showed it in simulation;28 Roch and Steel proved it.29

Where the methods fail

The failure modes are published, named, and mostly found by the people who use the methods.

Long-branch attraction.  Felsenstein proved in 1978 that parsimony can be positively misleading.30  Two lineages that have each changed a lot get pulled together on the tree.  They have both had many chances to reach the same state by accident.  His condition is explicit:  parsimony is consistent only when P² < Q(1−Q), where P is the chance of change on the two long branches and Q on the short ones.30

Felsenstein also wrote the sentence that should discipline both sides of this debate:

"it does little good simply to point out that the evolutionary models employed here do not apply. ... That amounts to a confession of ignorance rather than validation of the inference method in question."30

Alignment.  Before comparing position 400 across ten species, you must decide which positions correspond.  That is a separate inference, and it moves answers.  Wong, Suchard and Huelsenbeck ran 1,502 yeast genes through seven alignment programs.  46.2% gave a different tree depending on the aligner alone.31  There were only seven species, so only 945 trees to disagree about.  The mitigating detail, which their own figure caption gives, is that the conflicting trees had low support.

Model choice matters less than people think.  This cuts against a common creationist charge — that researchers tune the model until the tree looks right.  Abadi and colleagues simulated 7,200 datasets.  Fitting the true model recovered the exact tree 50.17% of the time.  Fitting a general model got 50.82%.  Fitting a deliberately wrong model got 48.31%.32  Fitting the correct model was no better than fitting a generic one.  On real datasets AIC and BIC pick different models 62% of the time, yet still return the same tree 79% of the time; on simulated data no two criteria ever disagree on more than 17% of trees.32  If you want to attack this inference, this is the wrong place.

Rooting is the weakest step.  An unrooted tree has no direction, so no ancestors and descendants.  The root usually comes from an outgroup, and the outgroup is by construction the most distant thing in the analysis — the longest branch.  Shavit and colleagues found a single outgroup "frequently disrupts an otherwise correct ingroup phylogeny," that two outgroups beat one, and that "the ingroup was most accurately recovered when no outgroup was used."33

Is the data filtered until it agrees?

This is the most serious methodological charge creationists make.  One version of it is fair.

The charge usually cites a 2013 paper titled Pruning Rogue Taxa Improves Phylogenetic Accuracy.34  This project's own note file paraphrases it as "when we remove the data that doesn't fit well our phylogenetic tree looks better."35

That paraphrase does not survive reading the paper.  Aberer, Krompass and Stamatakis tested the method on 400 simulated datasets where the true tree is known by construction.  They ran the control that matters.  Pruning the same number of taxa at random made agreement with the true tree worse.  Informed pruning made it better.34  A procedure tested against a known answer, with a randomisation control, is not circular.  The paraphrase is withdrawn here.

A narrower criticism survives, and the paper states it.  Rogue taxa "are mostly simply excluded (pruned) in current studies."34  So the published tree has no answer for exactly the taxa whose position was disputed.  The uncertainty is removed from the picture rather than resolved.

Gene filtering is more mixed.  Shen and colleagues found in 2021 that removing conflicting genes helps when lineage sorting and estimation error are low.  At medium-to-high levels it "eliminated or extensively reduced incongruence, but the resulting congruent species phylogenies were not always topologically identical to the true species trees."36  That is manufactured agreement without accuracy.  It is exactly what critics allege, under conditions critics do not specify.

The most useful audit is by Salichos and Rokas, who ran the standard fixes on their own data and reported the failures.  Dropping gap-containing sites made things worse.  Using the 100 slowest-evolving genes made things much worse.  Selecting genes by average support improved 15 branches and worsened none.8  A filter chosen on a criterion independent of the answer, which improves the whole tree at once, is a different thing from a filter chosen to get a preferred result.  The difference is testable case by case.

What "95% support" actually means

Nearly every published tree has numbers on its branches.  Few readers know what they measure.  Several different things get confused.

The bootstrap resamples the character columns, rebuilds the tree, and reports how often each grouping comes back.1  Felsenstein said plainly what it is a confidence interval for.  Not the true tree — the tree you would estimate on repeated sampling from the same pool of characters.  His next sentence matters:  "As such it may be misleading if the method used to infer phylogenies is inconsistent."1  A biased method produces a confident bias.

The familiar rule that 70% means 95% confidence comes from one 1993 paper, and it is conditional in the original.  It holds "under conditions of equal rates of change, symmetric phylogenies, and internodal change of ≤20% of the characters."  The same abstract says that with unequal rates, "bootstrap proportions >50% are overestimates of accuracy."37  Those conditions are rarely checked.

Bayesian posterior probabilities are overconfident, and by how much depends on what is wrong with the model.  Erixon and colleagues give the cleanest comparison.  Among branches supported above 95%, 0.2% were wrong under likelihood bootstrap with the right model, and 2.0% with the model too simple.  1.7% were wrong under Bayesian probability with the right model.  15.6% were wrong when the model was too simple.38  There are two known causes.  When a true internal branch has zero length, the probabilities for the three possible answers do not settle at one-third each — they stay polarised.39 40  And the prior on branch lengths matters more than anyone expects:  Yang and Rannala showed that on one fixed three-species dataset, sweeping that prior moves the answer from (1/3, 1/3, 1/3) to (0.925, 0.052, 0.023).41

At genome scale the bootstrap stops carrying information.  This is the most important point in this section:

"For genome-scale data sets, the site-sampling variance reduces to virtually zero.  In contrast, the systematic error does not appear to be reduced (and may even increase) with the increasing size of data sets and the antiquity of the species included, so it becomes the central determinant of statistical support."5

In their own example, 992,000 base pairs drove every branch but two to 100%.  When a modern paper reports 100% support everywhere, that says the dataset is large.  It does not say the tree is right.

Phillips, Delsuc and Penny showed the consequence on real data.  On the same eight yeasts and 106 genes behind the famous Rokas result, a different method returned a different tree — also at 100% bootstrap.  The cause was GC-content bias:  the alternative branchings each cost 688 extra GC/AT changes, four and a half times the 152-change margin by which ordinary nucleotide distances favoured the branch they picked.42  Recoding the data to remove the bias restored the original tree.  Two contradictory answers, both at maximum confidence, from one dataset.

This is why the field moved on.  The replacements are concordance factors.  The gene concordance factor is the percentage of gene trees containing a branch.  The site concordance factor is the percentage of informative sites supporting it.43  Minh, Hahn and Lanfear put it in one sentence:

"resampled datasets may always return the same tree (i.e., 100% bootstrap support), even though incomplete lineage sorting or other processes that lead to genealogical discordance are at work (e.g., gCF and sCF values ≪100%)."43

On their own 3,220-locus lizard dataset, bootstrap is 100% on every branch but one and the smallest gene concordance factor is 4.5%.  Note also that a site concordance factor's no-information baseline is 33%, not zero.  Any sentence of the form "this branch has 95% support" that does not name the statistic has said nothing.

The calibration that settles it

Arguing about what support numbers mean is less useful than measuring it.  The measurement exists.  Build a tree in a laboratory where the answer is recorded, then hand the data to the methods blind.

Hillis, Bull, White, Badgett and Molineux did this in 1992 with bacteriophage T7.  Nine lineages were grown down a known tree under a mutagen, with the ancestor frozen at every node.44  All five methods recovered the correct tree.  Every internal node had 100% bootstrap support but one, which had 92%.  Parsimony reconstructed 97.3% of 1,407 ancestral character states correctly.

That result should be printed at full strength.  It should also carry the authors' own caveat, or it is being reported selectively.  They chose the easiest possible tree and said so.  The symmetric equal-distance shape "has proven especially amenable to accurate reconstruction in theoretical studies ... and may thus be regarded as a 'null' model against which other topologies may be compared."44  On the same known tree, every method returned a wrong answer from restriction-fragment data, and four of six were wrong on 1,091 base pairs containing only 63 variable sites.45

Then the result to quote whenever support is treated as proof.  Cunningham, Zhu and Hillis rebuilt four-species trees from the same phage lineages under long-branch conditions:

"When parsimony was applied, the incorrect phylogeny was supported in 96% of the bootstrap pseudoreplicates."6

A better model recovered the right tree at 86%.  That is 96% support for a demonstrably wrong tree, where the truth was written down in advance.

The morphological version is worse.  Collard and Wood ran hominin-style craniodental cladistics on apes and baboons whose molecular relationships are known independently.  Every one of the seven clades at or above the 70% threshold was false, at 98%, 97%, 95%, 92%, 78%, 73% and 73%.7  Their own words:

"craniodental data can return impressive levels of statistical support (e.g., 97%) for patterns of phylogenetic relationship that are most likely incorrect."7

Strait and Grine replied that adding fossils fixes it.  Collard and Wood replied that Strait and Grine's own support values do not clear 70%.46 47

Now the counterweight, because taking only the failures from this section would mislead.  Sayyari and Mirarab measured a newer support statistic against known simulated trees.  Branches at local posterior probability ≥ 0.99 were correct 100% of the time when the gene trees were error-free, and 98% to 99.8% when they were estimated.48  Support statistics are not worthless.  The ones designed for genome-scale data are well calibrated.  The one everybody still prints is not.

How congruence is measured, and against what

Now the real question.  Two datasets sharing no characters are used to build two trees of the same organisms.  How much do they agree, and how impressed should anyone be?

The second half is the part people skip, and it decides the argument.

The founding experiment

Penny, Foulds and Hendy designed the test in 1982.  They took five proteins — cytochrome c, haemoglobin α, haemoglobin β, fibrinopeptide A and fibrinopeptide B — from eleven mammals.  They built a tree from each, then measured how far apart the trees were against what randomly chosen trees would give.49  Their logic runs through this whole article.  Agreement is only evidence if disagreement was possible.

What random looks like

Two trees can be compared by counting the groupings in one but not the other.  That is the Robinson-Foulds distance.50  The question is what that count looks like when the trees have nothing to do with each other.

Measured here, by drawing pairs of trees at random:

Taxa Max possible RF Mean RF between random trees As % of max Expected shared groupings
10 14 13.52 96.6% 0.24
20 34 33.67 99.0% 0.16
30 54 53.71 99.5% 0.15
50 94 93.72 99.7% 0.14
100 194 193.73 99.9% 0.13

Two random trees of fifty organisms share, on average, one-seventh of one grouping out of forty-seven available.  Simulation puts the chance that they share even one internal branch at about one in eight.

That is the baseline.  It is why the congruence argument has force.

The metrics are worse than they look

Robinson-Foulds is the standard, and it is a poor measure of how close two trees are.  Martin Smith showed that 86% of 100,000 random 11-leaf trees sit at the maximum possible distance from a reference tree, and that the measure returns fewer than eight distinct values across all 100,000.51  Two trees differing only in where one species sits can be at maximum distance.

The better statistic, used for this article's own measurements, is the quartet statistic.  Take every set of four species and note which pairing each tree asserts.  Its null is known exactly:  two random trees agree on exactly one-third of quartets, because there are three possible answers.  So define

QS = (fraction of quartets both trees resolve the same way − 1/3) ÷ (2/3)

QS = 0 is chance.  QS = 1 means the two trees are identical.  Every measurement below uses it.

Two warnings belong on the record.  The ILD test, the standard significance test for "do these datasets disagree?", is weaker than its use implies.  Darlu and Lecointre's finding is mainly about power, not false alarms:  with 100 sites the test rejected congruence in no more than 5.8% of simulations against a nominal 5%, but it has little power to detect real conflict unless characters are numerous and rates are similar across datasets.52  And the homoplasy indices used throughout the artifact literature move with the number of taxa — in opposite directions.  Sanderson and Donoghue showed the consistency index falls as taxa increase.53  Steell, Hsiang and Field measured both and found that "CI detected more absolute homoplasy as number of taxa was increased, and RI and RHI both detected less relative homoplasy as number of taxa was increased."54  Adding taxa lowers CI and raises RI.  That matters below, though not in the direction it first appears:  the headline artifact-versus-biology comparison uses RI, and on taxon count its two halves are nearly matched — 25.0 taxa on average against 22.4.  The size problem in that comparison is character count, 1,228 against 94, and no published measurement says how RI moves with that.  Any comparison between a twelve-object artifact study and a two-hundred-species genomic one is invalid unless taxon counts are matched.  That control is applied below, and several published comparisons do not survive it.

How much conflict is there really

Both sides get this wrong, in opposite directions, and the same papers contain the correction.

The famous results are real.  Salichos and Rokas built 1,070 gene trees from 23 yeasts and got 1,070 different topologies, none matching the tree from the combined data.8  The avian genome project reported that not one of its 14,536 loci — 8,251 exons, 2,516 introns and 3,769 ultraconserved elements — matched its main tree.55

They mean less than they sound like.  "Incongruent" in both papers means not identical at every branch, with no support threshold.  One weak rearrangement anywhere disqualifies a gene tree.  Salichos and Rokas's own figure gives the numbers that matter:  average normalised distance from gene trees to the species tree 0.40, among gene trees themselves 0.52, among 1,000 random trees 0.99.8

That triple is the most useful thing in this section.  Gene trees disagree.  They are also nowhere near random.

The headline is partly about tree size.  In the same paper, an 18-species vertebrate dataset had 15 of 1,086 gene trees identical to the combined tree.8  "No gene tree ever matches" is a statement about how many species were used.

Where the conflict lives is not disputed.  It sits at short, deep branches.  The bird paper describes "a wall of many (25 of 45; ~55%) very short internal branches (0.0006 to 0.002 substitutions per site), almost all at deep divergences within Neoaves."55  All nine of its below-100% branches were among the shortest.

Filter for supported conflict and the numbers collapse.  Per branch in the bird data, 0–38% of gene trees conflict at above 75% support, while only 0–20% support the same branches at that level.  The rest is noise.55  Across the ninety-one branches Wickett and colleagues label with a conflict percentage — 41 on their concatenation tree and 50 on their coalescent tree — the median is 6% and the maximum 16%.  On the concatenation tree alone the median is 9%.56 9  Conflicting gene trees are systematically the worse-supported ones.57

Much of the rest is predicted by common descent itself.  Genes sort differently from species when ancestral populations were large and speciation events close together.  Song and colleagues found this accounts for 77% of gene-tree variation across 447 mammalian genes.58  That is a prediction of the theory, not an embarrassment to it.

But not all of it.  Reddy and colleagues found the disagreement between the two big bird studies tracks data type — coding against non-coding sequence from the same genomes — more than species sampling, and blamed model error "difficult to detect."59  The deepest branch in the animal tree flips with the model.60 61  Pease and colleagues traced about 8.7% genome exchange between wild tomato lineages to hybridisation, not sorting.62  More data does not fix either.

Single genes can flip published, maximally supported deep branches — but only where the branch was already contested.  Shen, Hittinger and Rokas showed that removing the single most influential gene switched the answer at 3 of 17 contested branches.  Removing the top 1% of sites destabilised 13 of 17.  Both halves belong in the report:  removing a random gene changed nothing anywhere, and none of the six control branches moved under any perturbation.63

The anomaly zone is real and narrow.  Degnan and Rosenberg proved that for five or more species, every tree shape has branch lengths where the most common gene tree is not the species tree.64  Real cases have been found.65  No survey exists of how much of the tree of life sits inside it, and nobody should imply one does.

The best example of a branch sequence data could not settle is also the best example of settling it another way.  Whether comb jellies or sponges branch first among animals ran from 2008 to 2021 with both answers strongly supported depending on the model.  Chromosome-scale gene linkage settled it in 2023 — seven ancestral linkage groups, four requiring irreversible fusions on one stem, with a stated null:  "Across one hundred million randomized Hormiphora genomes, we never found syntenic signals comparable to those observed."13

Do designed things make a nested hierarchy?

This is the question the whole argument turns on, and most discussions of it are confused because "design" names three different processes.

Copying with modification.  An engineer starts from last year's model and changes it.  This produces a nested hierarchy for exactly the reason biology does — the objects have a genealogy.  Cars, Linux distributions, manuscripts and languages are mostly this.

Modular reuse by requirement.  A designer picks parts from a common catalogue according to what each product needs.  Similarity then tracks function, not ancestry.  Two different subsystems can imply two different groupings of the same objects.

Independent design to a specification.  No shared parts, no shared history.  No structure at all.

The creationist claim is the second one.  Almost every published study applying cladistics to artifacts measures the first.  That is why finding tree-like structure in Turkmen textiles settles nothing — copied things form trees, and nobody disputes it.

What the artifact literature found

The most-cited result says artifacts fit trees about as well as organisms do.  Collard, Shennan and Tehrani measured retention indices on 21 biological and 21 cultural datasets.  Biological mean 0.61, cultural mean 0.59, Mann-Whitney U = 215.5, p = 0.900.16

At face value that breaks the popular argument.  It does not survive inspection, for three reasons.

Their biological set is two populations pooled.  Split it as their own table allows and DNA datasets average 0.512, morphology and behaviour datasets 0.691, with the cultural datasets at 0.586 sitting between them.  The pooled null result hides an ordering.

The size confound is present in their own data.  They state that the retention index "is not affected by number of taxa or number of characters."  Within their biological half the rank correlation between the index and character count is −0.66 (Spearman; Pearson −0.37).  Their biological datasets average 1,228 characters.  Their cultural ones average 94.

The statistic cannot do the job.  Nunn, Arnold, Matthews and Borgerhoff Mulder ran 5,000 simulations sweeping the amount of horizontal borrowing.  Borrowing did not enter the regression tree for either index.  The consistency index was driven by the number of societies, explaining 83% of variance; the retention index by rate of change, explaining 77%.  Their conclusion:  "while high values may indicate phylogenesis, low values are uninformative… these metrics provide very low power."17

They also state the thing that matters most here.  Traits borrowed as a package "will tend to show a strong tree-like structure, albeit one that differs from the history of other cultural or genetic traits."  A strong tree that is the wrong tree is a predicted outcome, not a surprise.

The two artifact studies that pull hardest in opposite directions should both be printed.  Tehrani and Collard measured Turkmen textiles at consistency index 0.68 before conquest and 0.61 after, then noted published hominid cladograms at 0.58 to 0.71 — numerically identical.66  Against that, Tëmkin and Eldredge's analysis of Baroque cornets is usually reported as finding no hierarchical structure.  That paper could not be obtained for this article, its abstract does not say that, and the characterisation comes from third parties.  No number from it is used here.

One more measurement deserves to be better known.  The consistency and retention indices contradict each other on the same artifact dataset.  On seventeen Palaeoindian point classes coded with the same eight characters, the consistency index is 0.59 and the retention index 0.70.67 16  The Turkmen textiles point the same way — consistency index 0.68, retention index 0.44 — though there the two statistics come from two different analyses of the same tradition, five assemblages against six.66 16  Any argument resting on one index can be reversed by choosing the other.

The canoe experiment

There is one clean natural experiment, and it is the strongest published answer to the question.

Rogers and Ehrlich recorded 134 design traits of Polynesian canoes across eleven island groups, classifying each as functional or symbolic.68  Gray, Bryant and Greenhill then asked how well those traits fit the independently known language tree for the same eleven societies — a genealogy established from vocabulary, with no reference to canoes.18

Site-specific log-likelihoods, higher being better:

Data Fit to the true population tree
Basic vocabulary −2.89
Symbolic canoe traits −6.13
Randomised canoe data −6.30
Functional canoe traits −6.64

Functional design fit the true history worse than randomised data did.  And the functional and symbolic trait matrices from the same eleven objects are uncorrelated with each other — Mantel r = 0.00108, p = 0.5029.68

That is what design by requirement looks like when you can check it against the real genealogy.

The experiment run for this article

The published numbers above use different statistics on different-sized datasets, so they cannot be laid side by side.  This section runs one procedure over both domains.

The method

For every dataset the steps are identical.  Split the characters into two disjoint halves.  Build a tree from each half by neighbour-joining.  Measure QS between the two trees.  Repeat over many random splits and report the median.

Three controls are applied, because without them the comparison means nothing.

Matched species counts.  Tree statistics vary with the number of taxa, so every dataset is subsampled to the same size.

Matched character counts.  Genomes carry hundreds of thousands of characters and product feature tables carry a few dozen.  Without this control, "genomes are more congruent" could just mean "genomes have more data."

An empirical floor.  Each dataset is also run with every character column shuffled independently across taxa.  That destroys all shared structure while preserving every column's composition.  It should read zero, and it does — between −0.06 and +0.03 on every dataset tested.

The scale is anchored at both ends

Before measuring anything real, the pipeline was checked against cases with known answers.  Simulated descent with modification, on a tree whose branches are all long enough to resolve, reads QS = 1.000.  Characters drawn independently for each object — design with no shared anything — reads QS = −0.003, scattering between −0.09 and +0.09 across eighteen conditions.  Both ends are pinned, so a number in between means something.69

Result 1: designed objects can be strongly congruent

At 24 taxa, using every character each dataset has:

Dataset Kind Characters QS
Yeast, 1,233 genes genome 425,645 1.000
Linux distributions × packages designed 1,748 0.838
Hand-copied manuscripts copying 625 0.754
Laurasiatherian mammals genome 1,122 0.604
Web browsers × features designed 65 0.288
Text editors × features designed 38 0.072
Email clients × features designed 36 0.066
File systems × features designed 35 0.066
Databases × features designed 41 0.056
Version control × features designed 14 0.001

Linux distributions score 0.838 — higher than a mammal alignment cut down to 24 taxa to match, which reads 0.604.  Designed objects can produce a nested hierarchy, and a strong one.  That is a straightforward loss for the popular form of the common-descent argument, and it is reported first for that reason.

The small feature tables sit near zero.  But they also have only 14 to 65 characters, so they are underpowered rather than structureless.  The honest reading is that they are uninformative, not that they refute anything.

The matched-character control makes the point sharper and less comfortable.  Thin every dataset to 40 characters and yeast drops to 0.291, while the laurasiatherian mammals drop to 0.078 — below web browsers at 0.288.  Much of biology's apparent advantage is data volume.

Result 2: the Linux tree is the wrong tree

Here is where the comparison becomes decisive, and it is possible only because of one property no biological dataset has.

Linux distributions have a written-down genealogy.  People recorded at the time that Ubuntu was forked from Debian and Mint from Ubuntu.  The GNU/Linux Distribution Timeline lists a parent for 450 distributions.  So the question that cannot be asked of genomes can be asked here:  when two independent halves of the evidence agree on a tree, is that tree the actual family history?

Thirty-seven distributions have both package data and a documented parent.  Measured at 24 of them, over eight independent subsamples:

Measurement QS
Package half against package half (self-agreement) 0.849 – 0.981, median 0.91
Package tree against the documented fork history 0.078 – 0.228, median 0.15
A random tree against the documented fork history 95th percentile 0.097 – 0.151, median 0.11

The package-derived tree is highly reproducible and it is very nearly the wrong tree.  It matches the real fork history at roughly the level a random tree reaches one time in twenty, and in two of eight subsamples it fell below that.

The reason is not mysterious.  Which packages a distribution ships tracks what it is for — desktop or server or embedded, minimal or comprehensive, which init system and which C library its maintainers prefer.  Those are requirements, not ancestry.  This is modular reuse by requirement, caught in the act, in the one system where the real answer was checkable.

Result 3: congruence keeps climbing for genealogies

The single-number comparison is vulnerable to the objection that genomes simply have more data.  The test that survives it asks about the shape of the curve rather than its height.

As characters are added, does agreement between independent halves climb toward 1, or flatten out below it?

A real genealogy predicts convergence.  Every character was marked by the same branching history, so more characters means a better estimate of one true tree, and the halves must agree in the limit.

Measured at 24 taxa, except the chloroplast alignment, which has only 19 sequences and is measured on all of them:

Characters Yeast genes Chloroplast DNA Linux distros Manuscripts Laurasiatherian
20 0.220 0.233 0.267 0.165 0.110
40 0.427 0.305 0.376 0.373 0.146
80 0.501 0.589 0.424 0.630 0.186
160 0.659 0.664 0.550 0.744 0.233
320 0.740 0.735 0.663 0.831 0.452
640 0.843 0.786 0.717 0.853 0.574
1,280 0.908 0.830 0.811 — —
all available 1.000 0.988 0.857 0.884 0.690

Every one of these climbs, including the designed one.  The Linux curve sits inside the biological range at every character budget and is above the mammal alignment for most of it.

So the convergence test does not separate the domains either.  What it separates is genealogies from non-genealogies — and Linux distributions have a genealogy, because forking is copying.  The pipeline is detecting that they were copied.  It is simply detecting the wrong copying: package sets record which maintainers borrowed from which, which is a real history, just not the fork history the timeline records.

What this does and does not show

Stated plainly, because it is easy to overclaim in either direction.

It shows that congruence alone does not establish genealogy.  A designed system produced QS = 0.93 self-agreement on a tree that was demonstrably not its family history.  Anyone arguing "independent datasets agree, therefore common descent" has to reckon with that.  This is a real point for the creationist side and it does not appear to have been made with a measurement before.

It does not show that biology's tree is wrong.  The Linux case shows the failure is detectable when an independent check exists.  For genomes, independent checks exist and they pass — which is the subject of the next section.

It does not show that designed objects generally look like genomes.  The one designed dataset that scored high is the one that is also a copying tradition.  The datasets that are closest to design-without-descent — product feature tables — scored 0.001 to 0.29, and the largest of those, web browsers, is itself full of forks.

The biggest limitation is the shortage of clean data.  There is no large, published, design-without-descent character matrix in existence that could be located for this article.  The product feature tables are crowd-sourced, unaudited, and carry 14–24% unknown cells.  That gap is the main reason this section is weaker than it should be, and it is stated rather than papered over.

Ewert's dependency graph

One design model has been fitted to real genomes and compared head-to-head with a tree.  Winston Ewert's 2018 paper models gene distribution as a software dependency graph — modules that products import — rather than as inheritance down a tree.20

His result is that the dependency graph wins by enormous margins on nine genomic databases and loses on five simulated-evolution controls.  That result was reproduced here from his own released files.  Every published figure reproduces in direction;  the reproduction runs 4% to 26% more favourable to the graph than his published numbers on eight of the nine, and 79% on UniRef-50 — the sign is expected from the estimation method;  the size on UniRef-50 is not explained.19  It is not a typo, a sign error, or a fabrication, and it appears to be the first independent reproduction by anyone.

Three things have to be said about it, and they cut different ways.

His tree model is fairer than critics assume on the axes usually attacked.  It allows gene loss with a free rate parameter. It is Dollo, so genes are gained once and lost many times.  And the tree receives no complexity penalty at all while the graph pays a real one — stated in the paper and verified in his released scoring code.20 19  He also takes the maximum over trees rather than an average, which he says himself "will bias the results in favor of common descent."

His tree model is a strawman on one axis, and he says so.  It excludes incomplete lineage sorting, gene flow, horizontal transfer, convergent evolution and gene resurrection, which "occur in nature, but are not included in this model."20  He denies this matters, on the same page:  "Mechanisms which produce deviations from the tree are not relevant to that prediction."  No gene-family birth-death model is fitted, though he cites that literature.  He is comparing against a tree model no working biologist holds, and he concedes the escape route in print:  "given the wide variety of mechanisms that might be invoked along with common descent, it would not be surprising that some combination of those could produce data which fit a dependency graph better than a tree."

The headline number is true and nearly useless.  "The data is over 103000 times more likely" is his least favourable database.  Expressed as a share of the tree model's own error, the graph explains 0.9% to 14.6% more of the data.  His own JavaScript control — where a dependency graph is literally the true generating process — scores 5.6%, higher than four of the nine real databases.  That denominator is the reproduced tree model's total surprisal.  Measured instead against Ewert's own published Bayes factors, the JavaScript control is the most graph-favourable dataset in the paper, not a middling one.  That reframing is not original here:  Joshua Swamidass made the point in 2018, in the only substantial public discussion the paper ever received, and the figures above are this article's own recomputation of it.70

There is also a control nobody has run.  Under his tree model the per-species gene-loss counts rank the two reference-quality genomes lowest — human 920, mouse 681 — and fragmentary assemblies highest, with macaque charged 3,252 losses since a 25-million-year split.  The largest non-taxonomic modules are "everything except C. elegans" and "everything except the two flies."19  That is what genome annotation error would look like.  His simulated control has no annotation error in it at all.

Nobody has replied in print in eight years.  The journal reserved a DOI for a critique; none exists.  Four genuine citations, none in phylogenetics, none critical.19  The paper was ignored, not refuted, and that cuts against the mainstream side.

Whether convergent traits share genes, substitutions and silent sites, which is the test that separates a tree from reuse case by case, is taken up in Convergent Evolution and Design Reuse.

Where the argument actually stands

If congruence alone does not establish genealogy, what does?

The answer is that the tree recovered from sequence is confirmed by evidence of a completely different kind — evidence where shared function cannot be the explanation, because the characters have no function.

Rare genomic changes.  A retrotransposon inserting at one specific position out of billions is very nearly impossible to repeat independently.  Fifty-three verified marsupial retroposon markers gave zero conflicts.12  That is a different kind of evidence from sequence similarity, and a designer reusing functional modules has no reason to reuse insertion addresses in a nested pattern.

The honest counterweight is that this class is not homoplasy-free.  In bears, 132,093 transposable-element insertions gave 100% bootstrap on every branch but one, and at one branch 278 and 265 insertions supported the two wrong topologies against 168 for the right one.71  The cause was lineage sorting, not repeated insertion — which is itself a prediction of common descent, but the raw counts still mislead.

Chromosome structure.  The comb-jelly result above rests on ancestral gene linkages requiring irreversible fusion events, with a permutation null over a hundred million randomised genomes.13

Fossil order.  Trees built with no stratigraphic input can be tested against the order fossils appear in rock.  Here the honest numbers are worse than usually quoted.  Of 384 cladograms, only 148 (39%) show a statistically significant rank correlation between branching order and fossil order.14  The commonly quoted 75% and 78% "pass rates" are threshold counts, not significance tests, and the stratigraphic consistency index sits at exactly 0.5 by construction on a fully balanced tree with distinct dates.72  This line of evidence is real but weaker than its reputation.

Morphology against molecules is the place the argument is most often overstated.  Springer and colleagues deleted the DNA of one living mammal order at a time and re-ran the tree on morphology alone.73  Re-run on a 4,541-character matrix twenty-four times larger, the majority of orders still moved.74  Brady and Springer count them:  13 of 18 orders (72%) landed in an interordinal position that contradicts the molecular scaffold, so it is not a sample-size problem.75  Even given fossils, reconstructed ancestors and a molecular scaffold, only 42% of orders kept their correct placement.75  The mechanism is not random noise:  "all three orders with highly specialized myrmecophages (Xenarthra, Pholidota, Tubulidentata) cluster with one of the other myrmecophagous orders when treated as pseudoextinct."74  Ant-eating builds the wrong tree.

That result matters for this site's own fossil articles, and it is carried into §6b rather than buried here.  Morphology and molecules agree at the ranks nobody disputes and disagree at the deep, short branches where the arguments live.

Baraminology on the same yardstick

Fairness requires putting the creationist alternative through the same audit.

The method's most damaging critique is creationist.  Colin Reeves, published by Answers in Genesis, concludes that baraminic distance correlation "commits a category error," that "the use of a t-test is completely invalid," that "we have no idea what the distribution of these values would be under either a null or an alternative hypothesis," and that it "really should be abandoned."22  Both halves belong in the report:  when Reeves replaced it with standard cluster analysis on Wood's own 30-species turtle dataset, he reproduced Wood's answer, including the detailed cluster composition.22

The answer depends on which species you include, and Wood showed it himself.  With the same 87 characters, a large sample gives significant correlation between all species in Homo; cutting to six species leaves one.  Even H. sapiens against H. neanderthalensis flips.76

There is no rule for choosing characters, and 38.4% of published analyses returned nothing.  Twenty-eight of 73 studies to 2009 gave "uncertain" results.  Wood's own words: "the criteria for selecting characters informative for baraminology remain unknown."3

Creationist analysts disagree flatly with each other on the same data.  Australopithecus sediba is human by four analyses, and a single reply paper by three authors calls that result a false positive without running an analysis of its own.  Homo naledi comes out human or not human depending on the analyst, and in one case the conclusion was preserved by downweighting the inconvenient characters by a factor of ten.77

The model's advocates have removed its falsifiability in print.  Wood in 2002:  failure to find discontinuity "should not [be considered] a failure."21  O'Micks in 2016:  "We cannot demand that God created each baramin to have significant morphological differences compared to other baramins."78

The orchard model concedes nested hierarchy inside a kind, and kinds are large.  One holobaramin of about 10,000 grass species.21  About 150 fossil horse species in one monobaramin.21  Chimpanzees, gorillas, Au. afarensis and Au. garhi in one kind.3  Everything below roughly the family level is common ground between the two models.

No published creationist model makes a risky prediction about congruence that differs from common descent.  The one genuine risky prediction on record — that a statistical bridge uniting humans and chimpanzees into one group would falsify the discontinuity hypothesis — is about a gap, not about congruence.  It has been tested, it did not trigger, and it does not discriminate, because mainstream phylogenetics predicts that gap too.

Todd Wood deserves the last word here, because he asked this article's question first and did not answer it:

"It would be easy to attribute their similarities to a common designer, but such an attribution would be trivial. ... what distinguishes a common designer from a common ancestor?"3

The measurement in this article is one attempt at an answer.  A common designer who copies is indistinguishable from a common ancestor by congruence alone.  A common designer who reuses modules by requirement is distinguishable, and produces the Linux pattern:  a stable tree that is not a history.  Genomes do not show the Linux pattern.

Common objections

"Trees can be built from any data, so a tree means nothing."  True and irrelevant as stated.  Every clustering algorithm returns a branching diagram.  The question is whether this diagram beats the alternatives by more than chance, and that is testable.  Shuffling each character column independently and re-running this article's pipeline gives QS between −0.06 and +0.03 on every dataset.  The four datasets with more than a few hundred characters give 0.60 to 1.00.  The difference is the evidence.  The small feature tables give 0.001 to 0.29, which is what too little data looks like.

"Different genes give different trees, so there is no tree of life."  The same papers reporting the conflict report distances of 0.40 to the species tree and 0.52 among gene trees, against 0.99 for random trees.  Filtered for supported conflict, the median across ninety-one land-plant branches is 6%.  The conflict is real, it is concentrated at short deep branches, and 77% of it in mammals is accounted for by ordinary population genetics.

"Horizontal gene transfer destroys the tree."  For prokaryotes it genuinely wrecks a single bifurcating representation, and Koonin's group says so:  the original concept "is obsolete."  The same paper measured the residual signal and found the nearly universal trees sit at 1.4% to 4.3% inconsistency against a random-tree baseline of about 80%.2  Both halves are true.  For animals and plants the transfer rate is orders of magnitude lower, and the headline paper claiming 128 foreign genes in the human genome was cut to a single-digit count of defensible candidates on re-analysis.

"Researchers filter out the data that disagrees."  Sometimes, and the practice divides.  Rogue-taxon pruning was validated against 400 datasets with known true trees, with a random-pruning control that made things worse.  Gene filtering helps under low conflict and manufactures false agreement under high conflict, and the paper showing that is a mainstream one.  Judge case by case; the criterion is whether the filter was chosen independently of the answer.

"A designer would reuse parts, which explains the pattern."  This is the best version of the objection and this article takes it seriously enough to measure.  It predicts the Linux pattern — high self-agreement, low genealogical accuracy, groupings that track function.  Genomes show high self-agreement and independent genealogical confirmation from characters that have no function to share.  The reuse hypothesis is not refuted, but it now has a signature, and genomes do not carry it.

"Bootstrap values of 100% prove the tree."  They do not.  At genome scale the bootstrap goes to 100% regardless, by a mechanism its own practitioners describe.  A documented case gives 96% support for a known-wrong tree.  Ask for concordance factors instead.

"Common descent is unfalsifiable because any conflict gets explained away as lineage sorting."  This has force.  Lineage sorting, transfer and model error are each real and each absorbs conflict.  The reply is that they are not free parameters — lineage sorting predicts how much conflict and where, and Song and colleagues tested that prediction quantitatively and got 77%.  But the objection is not empty, and this article does not pretend it is.

Sources


  1. 1.

    Felsenstein, Joseph.  “Confidence limits on phylogenies: an approach using the bootstrap.”  Evolution 39(4):783–791.  1985.  A methods paper that brought the bootstrap into phylogenetics: the characters are resampled with replacement, the tree is rebuilt on each sample, and a group that appears in 95% or more of the samples counts as statistically significant.  The paper says plainly that the bootstrap gives an interval for the tree that repeated sampling of characters would estimate, not the true tree.  The article uses it to say what a bootstrap value measures. ↩

    1. 1a.

      Page 786, right column:  “Bootstrapping provides us with a confidence interval within which is contained not the true phylogeny, but the phylogeny that would be estimated on repeated sampling of many characters from the underlying pool of characters.  As such it may be misleading if the method used to infer phylogenies is inconsistent.” ↩ ↩ ↩

  2. 2.

    Puigbò, Pere et al.  “Search for a 'Tree of Life' in the thicket of the phylogenetic forest.”  Journal of Biology 8:59.  2009.  A research article that examined a 'forest' of 6,901 gene trees from prokaryotes (bacteria and archaea) to see whether a central pattern survives horizontal gene transfer, the movement of genes between unrelated lineages.  The authors found that horizontal transfer is widespread, so very few gene trees agree fully, but 102 nearly universal trees still carry a consistent signal.  The article uses it to show both halves of the picture: transfer wrecks a single tree for prokaryotes, and a residual signal of inheritance remains. ↩

    1. 2a.

      “the original tree of life concept is obsolete: it would not even be a 'tree of one percent'.  Nevertheless, there seems to be a discernible signal of consistency between the trees in the forest of life, down to the deepest branching levels.”  Printed as 59.12, which is page 12 of article 59.  The study used 6,901 maximum-likelihood gene trees from 100 prokaryote genomes, and the 102 nearly universal trees have inconsistency scores of 1.4% to 4.3% against a random-tree baseline near 80%.  Consistency is high at shallow depths.  Page 11:  “abruptly drops, almost down to the level of random trees”.  This happens at the deepest branchings. ↩

  3. 3.

    Wood, Todd Charles.  “Baraminological Analysis Places Homo habilis, Homo rudolfensis, and Australopithecus sediba in the Human Holobaramin.”  Answers Research Journal 3:71–90.  2010.  A creationist paper that reran eight published cladistic studies of fossil and living hominids with statistical baraminology, a method that groups organisms into created kinds by how alike their characters are.  It found that the hominins split into as many as four holobaramins, with the genus Homo including Australopithecus sediba, and it cites Wood's own 2009 survey, in which 28 of 73 earlier studies (38.4%) gave "uncertain" results.  The article cites it for the author's admission that the criteria for choosing characters are unknown. ↩ ↩

    1. 3a.

      Wood raises the objection against his own position and does not answer it, then quotes Darwin at length without rebuttal. ↩ ↩

  4. 4.

    Measured for this article.  research/sim/rf-null.py, 20,000 replicates per row, seed 1.  Trees drawn uniformly over labelled unrooted binary topologies by sequential random-edge insertion.  The expected number of shared non-trivial splits between two random trees falls from 0.24 at 10 taxa to 0.13 at 100 and converges on 1/8.  Tree counts are (2N−5)!!, computed exactly. ↩

  5. 5.

    Kumar, Sudhir et al.  “Statistics and truth in phylogenomics.”  Molecular Biology and Evolution 29(2):457–472.  2012.  “For genome-scale data sets, the site-sampling variance reduces to virtually zero.”  “In contrast, the systematic error does not appear to be reduced (and may even increase) with the increasing size of data sets and the antiquity of the species included, so it becomes the central determinant of statistical support.”  A review of how statistical confidence behaves when trees are built from whole genomes.  With that much data the random error shrinks almost to nothing, so a branch can carry very high support while the result is wrong, because the systematic bias of the method is what remains.  The article uses it to explain why the bootstrap support value, a resampling test that gives each branch a percentage, fails on large data sets. ↩ ↩

  6. 6.

    Cunningham, C. W. et al.  “Best-fit maximum-likelihood models for phylogenetic inference: empirical tests with known phylogenies.”  Evolution 52(4):978–987.  1998.  A laboratory experiment with bacteriophage T7 lineages grown along a known tree, so the correct answer is known in advance.  When terminal branch lengths were very unequal, parsimony (choosing the tree that needs the fewest changes) supported the wrong tree in 96% of bootstrap replicates, while the best-fit model supported the correct tree with 86% bootstrap support.  The article uses it as its documented case of a method giving high support to a known-wrong tree. ↩ ↩

  7. 7.

    Collard, Mark, and Bernard Wood.  “How reliable are human phylogenetic hypotheses?”  PNAS 97(9):5003–5006.  2000.  A research paper in which two anthropologists ran the methods used on fossil human skulls and teeth on apes and monkeys whose DNA trees are well established.  The skull and tooth trees disagreed with the DNA trees, and the authors concluded that little confidence can be placed in phylogenies built from higher-primate skull and tooth evidence alone.  The article uses it to show that bootstrap support can be high for trees that are false. ↩ ↩

    1. 7a.

      Page 5005:  “More problematically, the tests suggest that such analyses can strongly support phylogenetic hypotheses that are misleading.  For example, the bootstrap-based tests indicate that craniodental data can return impressive levels of statistical support (e.g., 97%) for patterns of phylogenetic relationship that are most likely incorrect.  In other words, cladistic analyses of higher primate craniodental morphology may yield not only 'false-positive' results, but false-positive results that pass, by a substantial margin, the statistical test favored by many researchers.”  The seven clades at or above 70% were 98, 97, 95, 92, 78, 73 and 73. ↩

  8. 8.

    Salichos, Leonidas, and Antonis Rokas.  “Inferring ancient divergences requires genes with strong phylogenetic signals.”  Nature 497:327–331.  2013.  A phylogenomic study of 1,070 orthologous genes (genes that descend from one ancestral gene) from 23 yeast genomes, testing how gene trees conflict with the tree built by joining all the genes into one data set, a method called concatenation.  The authors found that all 1,070 gene trees differ from the concatenated tree, that the conflict is worst at short branches deep in the tree, and that choosing genes or branches with strong support makes the result more robust, with similar results in vertebrate and animal datasets.  The article uses it to show that gene trees conflict a lot, and that the authors also report which corrections work. ↩ ↩

    1. 8a.

      Figure 1b:  Average normalised Robinson-Foulds distance of the 1,070 yeast gene trees to the species tree 0.40; among the gene trees themselves 0.52; among 1,000 random trees 0.99.  This is the single most useful triple in the gene-conflict literature, and it sits in the same paper as the headline figure of 1,070 distinct topologies. ↩ ↩

    2. 8b.

      The 18-taxon vertebrate dataset had 15 gene trees of 1,086 identical to the concatenation tree.  The famous result that no gene tree matches is partly a function of taxon number. ↩

    3. 8c.

      Figure 2, page 329:  The authors' own audit of standard remedies on their own data.  Excluding gap-containing sites lowered internode certainty at 7 internodes and raised it at 0.  Using the 100 slowest-evolving genes lowered it at 9 and raised it at 2.  Selecting the 131 genes with highest average bootstrap improved gene support frequency at 15 internodes and worsened none. ↩

  9. 9.

    Measured for this article.  Wickett et al. 2014 never state a median.  Their Figs. 2 and 3 print each branch as bootstrap,conflict %, where conflict is the percentage of gene trees disagreeing at above 75% bootstrap.  Extracting those labels from the PDF text layer (pdftotext -layout, then the pattern (\*|\d{1,3}),(\d{1,2})) gives 41 labelled branches on Fig. 2 with median 9% and 50 on Fig. 3 with median 4% — the 91st label is 2,6 beside Oryza sativa, whose one-digit bootstrap a two-digit pattern misses;  pooled, 91 branches with median 6%, minimum 1% and maximum 16%.  Re-extracted by the orchestrator on 2026-08-22 with the same command (41 + 50 = 91, median 6).  The label order is inferred from the caption, which names both quantities and gives no other candidate. ↩ ↩

  10. 10.

    Measured for this article.  research/sim/compare-domains.py and research/sim/distro-truth-test.py, seeds 1-8.  Package data from Repology (17,911 projects × 278 repositories); documented genealogy from the GNU/Linux Distribution Timeline, 450 distributions with a curated parent.  Thirty-seven distributions have both.  Scored at 24 taxa on the 6,600-7,500 quartets the documented genealogy resolves; about 35% of quartets are unresolved because ten of the thirty-seven are direct children of Debian, and an earlier version of this test scored those as disagreements and understated the match.  Raw output in research/results/distro-truth-test.csv. ↩ ↩

  11. 11.

    Measured for this article.  The Computer-Assisted Stemmatology Challenge heinrichi tradition — Piispa Henrikin surmavirsi, hand-copied in pencil by 17 volunteers under controlled conditions, 37 manuscripts × 1,208 loci, with the true copying graph released.  Split-half QS 0.754 at 24 taxa and 0.884 using all 668 informative characters.  The comparison against the released stemma is research/sim/stemma_run_truth.py, written up in research/stemma-truth-test.md with raw output in research/results/stemma-truth-test.csv:  neighbour joining on Hamming distance over the 690 informative characters of the 34 witnesses that appear in the stemma scores QS +0.487 against it and parsimony +0.846, where the truth resolves 40,578 of the 46,376 quartets and a null of 2,000 random labelled trees has a 95th percentile of +0.066;  the 21 witnesses of the 0.8837 run that appear in the stemma score +0.78.  The dataset carries a licence restriction:  analysed in house, not redistributed.  Note that published accuracy figures on this corpus are close to meaningless without their floor — the no-information baseline of the challenge's own scoring rule is 56.6%, and the published neighbour-joining submission's "direction of copying" score of 58.467% is exactly what a graph with no edges scores, which is why nothing here is claimed about who copied whom. ↩

  12. 12.

    Nilsson, Maria A. et al.  “Tracking marsupial evolution using archaic genomic retroposon insertions.”  PLoS Biology 8(7):e1000436.  2010.  A genome study of marsupial mammals, such as kangaroos and opossums, which screened about 217,000 retroposon-containing loci (places where copies of jumping DNA sit) from opossum and kangaroo genomes and kept 53 informative markers.  The authors report that the four Australasian orders share a single origin, with Microbiotheria as their closest relative, and that the South American opossums branch off first.  The article uses the 53 markers as its example of rare genomic changes, a different kind of evidence from sequence similarity. ↩ ↩

  13. 13.

    Schultz, Darrin T. et al.  “Ancient gene linkages support ctenophores as sister to other animals.”  Nature 618(7963):110–117.  2023.  A genome study that compared which genes sit on the same chromosome across animals and their single-celled relatives, a pattern called synteny.  Chromosome fusions that could not be reversed placed comb jellies as the first branch off from all other animals, and a shuffling test on randomized genomes showed that chance does not produce the same pattern.  The article uses it as the clearest case of a branch that sequence data could not settle and chromosome structure did. ↩ ↩ ↩

  14. 14.

    Benton, Michael J., and Rebecca Hitchin.  “Congruence between phylogenetic and stratigraphic data on the history of life.”  Proceedings of the Royal Society B 264(1383):885–890.  1997.  A survey of 384 published family trees (cladograms) for echinoderms, fishes and tetrapods, comparing the branching order of each tree with the order in which its fossils appear in the rock.  It found good agreement on a variety of comparative metrics, which the authors read as a sign that most trees are broadly accurate and that the fossil record is a reasonably faithful sequence.  The article uses it as the main test of trees against fossil order. ↩ ↩

  15. 15.

    Measured for this article.  research/sim/convergence-curve.py, seed 20260815, 15 random half-splits per cell, taxon count held at 24 throughout, except the 19-sequence chloroplast alignment, so the only thing changing is how much evidence each half receives.  Raw output in research/results/convergence-curve.csv. ↩

  16. 16.

    Collard, Mark et al.  “Branching, blending, and the evolution of cultural similarities and differences among human populations.”  Evolution and Human Behavior 27(3):169–184.  2006.  A comparative study that fitted the branching tree model used in biology to 21 biological datasets and 21 cultural datasets, then compared how well the model fitted each group.  The cultural datasets were not significantly less tree-like than the biological ones, which the authors read as evidence against the claim that blending has always dominated cultural evolution.  The article uses it as the most-cited result that artifacts fit trees about as well as organisms do, and it examines where that result fails. ↩ ↩ ↩ ↩

  17. 17.

    Nunn, Charles L. et al.  “Simulating trait evolution for cross-cultural comparison.”  Philosophical Transactions of the Royal Society B 365:3807–3819.  2010.  Page 3817:  “while high values may indicate phylogenesis, low values are uninformative… these metrics provide very low power.”  A simulation study of cultural traits that can pass sideways between groups as well as from parent to child.  Across 5,000 simulations the authors tested whether the consistency index and retention index can detect that sideways passing, and found that faster change lowers these indexes more than borrowing does, so a low value does not mean borrowing.  The article uses it to show that the two indexes are weak tests, because borrowing does not show up in them. ↩ ↩

  18. 18.

    Gray, Russell D. et al.  “On the shape and fabric of human history.”  Philosophical Transactions of the Royal Society B 365:3923–3933.  2010.  A paper for a theme issue on cultural and linguistic diversity that sets out two debates, about how tree-like human history is (its shape) and how coherent its lineages are (its fabric), and then tests them with computational methods from biology.  It argues that these debates are better seen as positions on continuous dimensions, and it compares how tree-like vocabulary evolution is in different parts of the world.  The article uses it for the Polynesian canoe comparison, in which both the functional and the symbolic canoe traits lie close to the random distribution on the language tree. ↩ ↩

  19. 19.

    Measured for this article.  research/sim/ewert_repro.py against Ewert's own released supplementary archive (sources/data/ewert-2018c_dependency-graph/supplemental.zip, md5 b7741f8021252a5d470fa3eca71fb14e).  All fifteen datasets reproduce in direction — nine biological databases, five simulated-evolution controls and the JavaScript control;  the ratio of reproduced to published log Bayes factor runs 1.04-1.79 on the nine biological databases, 0.82-0.93 on the five simulated controls and 1.08 on the JavaScript one.04-1.79 on the nine biological databases and 0.82-0.93 on the five simulated controls.  Surprisal share — the fraction of the tree model's own error that the graph removes — is 0.94% (UniRef-50), 3.12% (OMA), 3.52% (Pfam), 5.35% (HomoloGene), 5.49% (TreeFam), 8.81% (EggNOG), 10.03% (OrthoDB), 13.76% (Ensembl) and 14.55% (Hogenom), with the JavaScript control at 5.64%.  Input matrices verified to agree cell-for-cell between his independently fitted tree and graph files.  Raw output in research/results/ewert-repro.csv. ↩ ↩ ↩ ↩ ↩

  20. 20.

    Ewert, Winston.  “The Dependency Graph of Life.”  BIO-Complexity 2018(3):1–27.  2018.  Page 21:  “In the case of the dependency graph model, there are many possible dependency graphs, and we need to assign a probability to each one.  This probability is only used in the case of the dependency graph model.”  A research article from the Biologic Institute by an author who writes as a defender of intelligent design, proposing that the nested pattern of life comes from shared software-style modules rather than from common descent.  It tests predictions of this model against biological and synthetic data, and the paper says that, as the first evaluation of the hypothesis, its tests are limited.  The article uses it as the one design model fitted to real genomes and compared head-to-head with a tree. ↩ ↩

    1. 20a.

      Page 9:  The paper lists incomplete lineage sorting, gene flow, horizontal gene transfer, convergent evolution and gene resurrection as mechanisms that.  “occur in nature, but are not included in this model.”  Pages 5–6:  The concession reads.  “given the wide variety of mechanisms that might be invoked along with common descent, it would not be surprising that some combination of those could produce data which fit a dependency graph better than a tree.” ↩ ↩

  21. 21.

    Wood, Todd Charles.  “A baraminology tutorial with examples from the grasses (Poaceae).”  TJ (Journal of Creation) 16(1):15–25.  2002.  A short tutorial on baraminology, the creationist method for finding "created kinds" (baramins) by how alike organisms are, using the grass family Poaceae as its worked example.  The article cites it for how large a created kind is taken to be: one holobaramin of about 10,000 grass species, and about 150 fossil horse species in one monobaramin. ↩

    1. 21a.

      On failing to find discontinuity, Wood wrote:  “should not [be considered] a failure.”  Paired with O'Micks 2016, who wrote: "We cannot demand that God created each baramin to have significant morphological differences compared to other baramins." This removes the model's falsifiability in its advocates' own words. ↩

    2. 21b.

      The grass holobaramin of roughly 10,000 species is the worked example the tutorial is built around, and it is the clearest published statement of how large a created kind is taken to be. ↩ ↩

  22. 22.

    Reeves, Colin R.  “A Critical Evaluation of Statistical Baraminology: Part 1 - Statistical Principles.”  Answers Research Journal 14:261–269.  2021.  Page 263:  “commits a category error”.  “the use of a t-test is completely invalid”.  “we have no idea what the distribution of these values would be under either a null or an alternative hypothesis”.  “really should be abandoned”.  A critique of baraminic distance correlation, the statistical method that creationist researchers led by Wood use to test whether species fall into created kinds (baramins).  Reeves, a mathematician at Coventry University, argues that the method counts matching traits between every pair of organisms, runs a t-test on the results, and rests on a shaky understanding of statistical principles, so its use should be abandoned.  The article uses it as the most damaging critique of the method, and it comes from inside the creationist camp, in a creationist journal. ↩ ↩

    1. 22a.

      The second part of the same critique, by the same author, on alternatives to the method and on its conceptual and practical problems, which the article cites for the alternatives it tests.  The citation is Reeves, C. R., A Critical Evaluation of Statistical Baraminology: Part 2-Alternatives and Conceptual and Practical Issues, Answers Research Journal 14:271-282, 2021, at https://assets.answersresearchjournal.org/doc/v14/statistical_baraminology_alternatives_practical_issues.pdf. ↩

  23. 23.

    Doolittle, W. Ford.  “Uprooting the Tree of Life.”  Scientific American 282(2):90–95.  2000.  A popular magazine feature in which Doolittle argues that gene sequences are reshaping the tree of life and that its early branches and root are in doubt.  The article uses it to show that a quotation in wide circulation is cut off from the sentence that follows it. ↩

  24. 24.

    Saitou, Naruya, and Masatoshi Nei.  “The Neighbor-joining Method: A New Method for Reconstructing Phylogenetic Trees.”  Molecular Biology and Evolution 4(4):406–425.  1987.  A methods paper that introduces neighbour-joining, a fast way to build a phylogenetic tree from a table of distances between organisms, by repeatedly joining the pair of taxa that gives the smallest total branch length.  The authors tested it by computer simulation against five other tree-making methods, and found that neighbour-joining and Sattath and Tversky's method were generally better than the rest at recovering the correct unrooted tree.  The article uses it as the standard distance method, the one its own software uses, and it sets its predecessor UPGMA against it. ↩

    1. 24a.

      Table 6, page 420:  With rates varying among lineages, UPGMA recovered the true 8-OTU unrooted topology in 0 of 100 replicates at 600 nucleotides, against 75/100 and 72/100 for neighbour-joining on the same simulated trees. ↩

  25. 25.

    Felsenstein, Joseph.  “Evolutionary Trees from DNA Sequences: A Maximum Likelihood Approach.”  Journal of Molecular Evolution 17:368–376.  1981.  A methods paper that gave a computationally feasible way to find the maximum likelihood tree from DNA sequences under a simple model of DNA change, with a computer program to do it.  Unlike parsimony, which can give misleading results when rates of change differ between lineages, the method can test whether rates are constant and gives a rough indication of the error in the tree.  The article uses it to define maximum likelihood, the workhorse method. ↩

  26. 26.

    Mirarab, S. et al.  “ASTRAL: genome-scale coalescent-based species tree estimation.”  Bioinformatics 30(17):i541–i548.  2014.  A methods paper that introduces ASTRAL, a fast program that estimates a species tree from many gene trees when the conflict between them comes from incomplete lineage sorting, where ancestral gene versions have not yet sorted into their species.  The authors report that it is statistically consistent, runs on datasets of thousands of genes, and beats the other summary methods in every simulated condition they tried, though it is often more accurate than concatenation only when incomplete lineage sorting is high and there are enough gene trees.  The article cites it as an example of a species-tree method, which treats each gene tree as noisy evidence about the species tree rather than as a measurement of it.

    1. 26a.

      ASTRAL is a species-tree method, so it treats each gene tree as noisy evidence about the species tree.  The authors describe it as statistically consistent and able to handle datasets with thousands of genes. ↩

  27. 27.

    Liu, Liang et al.  “A maximum pseudo-likelihood approach for estimating species trees under the coalescent model.”  BMC Evolutionary Biology 10:302.  2010.  A statistics paper that proposes MP-EST, a way to estimate a species tree from many gene trees under the multispecies coalescent model, which describes how gene lineages sort through ancestral populations.  The authors tested it on simulated data and on a placental mammal dataset, where its four major clades matched those of a Bayesian concatenation tree.  The article cites it with ASTRAL as a species-tree method, one that treats each gene tree as noisy evidence about the species tree.

    1. 27a.

      The method this paper introduces, MP-EST, estimates the species tree directly from many gene trees.  The authors show it gets the right answer more and more reliably as the number of genes grows, even when the species tree is in the anomaly zone. ↩

  28. 28.

    Kubatko, Laura Salter, and James H. Degnan.  “Inconsistency of Phylogenetic Estimates from Concatenated Data under Coalescence.”  Systematic Biology 56(1):17–24.  2007.  A simulation study of concatenation, which joins all the genes into one long sequence before building a tree, when the coalescent (a model of how gene lineages sort in populations) makes individual gene trees conflict.  It found that concatenation becomes statistically inconsistent as genes are added, so more data drive it toward the wrong tree, and that the bootstrap can give moderate to strong support to that wrong tree.  The article uses it as the simulation evidence that gluing genes together can mislead. ↩

  29. 29.

    Roch, Sebastien, and Mike Steel.  “Likelihood-based tree reconstruction on a concatenation of aligned sequence data sets can be statistically inconsistent.”  Theoretical Population Biology 100:56–62.  2015.  A mathematical proof that building a tree by maximum likelihood from genes joined into one long sequence (concatenation) is statistically inconsistent, which means it converges on the wrong species tree as the data grow.  The authors work under the multispecies coalescent, the model of how gene trees arise inside a species tree, and say the proof gives a formal basis to earlier simulation results by Kubatko and Degnan.  The article uses it as the proof that concatenation can keep returning the wrong tree as more genes are added.  Mirror:  arXiv ↩

  30. 30.

    Felsenstein, Joseph.  “Cases in which Parsimony or Compatibility Methods will be Positively Misleading.”  Systematic Zoology 27(4):401–410.  1978.  A theoretical paper showing that parsimony (choosing the tree that needs the fewest changes) and compatibility methods can be inconsistent: when rates of change differ between lineages, they converge on the wrong tree as more data are added.  It works out the exact conditions for simple three- and four-species cases.  The article uses it for long-branch attraction, in which two lineages that have each changed a lot are pulled together on the tree. ↩

    1. 30a.

      Page 407:  Equation (18), for unrooted Wagner parsimony on a four-taxon tree, consistency requires P2 < Q(1−Q), where P is the probability of change on the two long non-sister branches and Q on the short branches. ↩

    2. 30b.

      Page 408:  The sentence is aimed at defenders of parsimony, and it applies equally to anyone dismissing a phylogenetic result by naming an assumption without showing that violating it changes the answer. ↩

  31. 31.

    Wong, Karen M. et al.  “Alignment Uncertainty and Genomic Analysis.”  Science 319(5862):473–476.  2008.  A research report that tested how much the choice of sequence alignment (the step that lines up matching positions across genes) changes later results.  The authors aligned 1,502 sets of matching genes from seven yeast species with seven alignment programs, and 46.2% of the genes gave one or more different trees depending on the program used.  Trees from dissimilar alignments also had lower bootstrap support, a measure of confidence in a branch, and the article cites it for the point that alignment choice can change the tree. ↩

  32. 32.

    Abadi, Shiran et al.  “Model selection may not be a mandatory step for phylogeny reconstruction.”  Nature Communications 10:934.  2019.  A simulation and real-data study of the step where software picks the mathematical model of DNA change before it builds a tree.  The authors found that the standard selection criteria often pick different models yet give nearly the same tree, and that skipping the step and using the most parameter-rich model, GTR+I+G, gave similar trees.  The article uses it against the charge that researchers tune the model until the tree looks right. ↩

    1. 32a.

      Section Results:  AIC and BIC selected different models in 62% of empirical datasets, yet the minimum pairwise agreement on inferred topology between any two selection criteria was 83% over 7,200 simulated datasets, and all distances to the true tree were statistically indistinguishable (Wilcoxon signed-rank, Bonferroni corrected). ↩

  33. 33.

    Shavit, Liat et al.  “The Problem of Rooting Rapid Radiations.”  Molecular Biology and Evolution 24(11):2400–2411.  2007.  A computer simulation study that generated DNA sequences on trees shaped like a rapid radiation, a burst of many species splitting in a short time, and tested which methods rebuilt the right tree.  It found that several standard methods were biased toward specific tree shapes, and that rooting with one outgroup (a distant species used to set the direction of the tree) often broke an otherwise correct tree.  The article uses it to show that rooting is the weakest step in tree building. ↩

  34. 34.

    Aberer, Andre J. et al.  “Pruning Rogue Taxa Improves Phylogenetic Accuracy: An Efficient Algorithm and Webservice.”  Systematic Biology 62(1):162–166.  2013.  A methods paper that presents a faster algorithm for finding rogue taxa, which are species whose place in a set of trees is unstable, and a free web service that runs it.  On simulated datasets whose true tree is known, removing the rogues made the bootstrap consensus trees and the maximum-likelihood trees closer to the true trees.  The article uses it to test the charge that dropping badly fitting taxa makes a tree look better than the data allow. ↩

    1. 34a.

      Section Accuracy, page 165:  The study used 400 simulated datasets, described as.  “for which the 'true' tree is known.”  Pruning an equal number of taxa at random had a negative effect on congruence with the true tree, while informed pruning improved it, with support gain tracking truth gain linearly. ↩

    2. 34b.

      Page 162:  “are mostly simply excluded (pruned) in current studies.”  This is the limitation that survives:  the published tree carries no answer for the taxa whose placement was contested. ↩

  35. 35.

    Working note on the tree of life: Pruning Rogue data.  n.d.  A section of this project's own working notes on the tree of life.  Its paraphrase of the paper it discusses is withdrawn.  The article cites it for the paraphrase of the paper that the rogue-taxa charge usually cites. ↩

  36. 36.

    Shen, Xing-Xing et al.  “Dissecting Incongruence between Concatenation- and Quartet-Based Approaches in Phylogenomic Data.”  Systematic Biology 70(5):997–1014.  2021.  A study that compared two ways of building a species tree from genome-wide data, using three published datasets from animals, fungi and plants.  It found that 30–36% of the genes in each dataset point the two methods to different answers, and that dropping those genes helped when gene-tree error and lineage sorting (genes whose history differs from the species history because of ancestral variation) were low.  At higher levels the disagreement disappeared, but the trees that remained were not always the true ones, and the article cites this for agreement without accuracy. ↩

  37. 37.

    Hillis, David M., and James J. Bull.  “An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis.”  Systematic Biology 42(2):182–192.  1993.  A computer-simulation and laboratory study, using a laboratory-generated phylogeny, of whether bootstrap values from parsimony analyses measure what they are usually taken to measure.  It found bootstrap values unbiased but very imprecise as measures of repeatability, and found that a value of 70% or more usually meant at least a 95% chance that the clade is real only under equal rates of change, symmetric trees and at most 20% internodal change; with unequal rates, values above 50% overestimate accuracy.  The article uses it to show that the familiar 70% rule is conditional, and that its conditions are rarely checked. ↩

  38. 38.

    Erixon, Per et al.  “Reliability of Bayesian posterior probabilities and bootstrap frequencies in phylogenetics.”  Systematic Biology 52(5):665–673.  2003.  A simulation study that analysed the same simulated DNA data by maximum likelihood bootstrapping and by Bayesian analysis, using identical models, so the true tree is known.  Bayesian posterior probabilities were higher than bootstrap frequencies (how often a branch recurs when the data are resampled) for true clades, but they also led to more wrong conclusions, and both kinds of error grew sharply when the model was too simple.  The article uses it for the clearest comparison of how often well-supported branches turn out to be wrong. ↩

  39. 39.

    Lewis, Paul O. et al.  “Polytomies and Bayesian phylogenetic inference.”  Systematic Biology 54(2):241–253.  2005.  A simulation study of Bayesian tree-building, a method that gives each branch a posterior probability.  The authors simulated data from a star tree, which has no true internal branches, in 25 replicates of 1 site or 100,000 sites, and in replicate 11 one arbitrary resolution reached a posterior probability of 0.9892.  The article uses it for the case where a true internal branch has zero length, so that the three possible answers do not settle at one-third each. ↩

  40. 40.

    Yang, Ziheng.  “Fair-balance paradox, star-tree paradox, and Bayesian phylogenetics.”  Molecular Biology and Evolution 24(8):1639–1655.  2007.  A theory paper that approximates Bayesian posterior probabilities for three species with binary characters and checks the result by simulation.  It confirms the star-tree paradox: when the true tree is a star with no resolved internal branch, the probabilities for the three possible trees do not settle at one-third each as data grow, and they stay very sensitive to the prior.  The article cites it for why those probabilities stay polarised when a true branch has zero length. ↩

  41. 41.

    Yang, Ziheng, and Bruce Rannala.  “Branch-length prior influences Bayesian posterior probability of phylogeny.”  Systematic Biology 54(3):455–470.  2005.  A simulation and real-data study of how the prior (the starting assumption about how long each internal branch of the tree is) changes Bayesian posterior probabilities, which are scores for how likely a branch is to be true.  The authors found that priors assuming long internal branches gave high probabilities to trees, and that the bootstrap and the posterior probability are different measures of accuracy.  The article cites it for the effect of the branch-length prior on Bayesian probabilities. ↩

  42. 42.

    Phillips, Matthew J. et al.  “Genome-Scale Phylogeny and the Detection of Systematic Biases.”  Molecular Biology and Evolution 21(7):1455–1458.  2004.  A short letter that re-analysed a dataset of 106 genes from eight yeast species, which had given a tree with 100% bootstrap support on every branch.  The authors found that a different method, minimum evolution with a distance measure, gave a different tree, also at 100% bootstrap support, and traced the cause to a GC/AT bias, a skew in the mix of bases that can mislead tree-building.  Recoding the bases as purines and pyrimidines reduced the bias and reinforced the original tree.  The article uses it as a real-data case of two contradictory answers, both at maximum bootstrap support, from one dataset. ↩

  43. 43.

    Minh, Bui Quang et al.  “New methods to calculate concordance factors for phylogenomic datasets.”  Molecular Biology and Evolution 37(9):2727–2733.  2020.  A short methods paper that adds two agreement measures to the IQ-TREE software.  The gene concordance factor is the percentage of gene trees that contain a branch, and the site concordance factor is the percentage of informative sites that support it.  The authors tested both on ten empirical datasets and argue that they show what the bootstrap leaves out.  The article uses it to explain why the field moved on from the bootstrap to these measures. ↩

    1. 43a.

      “resampled datasets may always return the same tree (i.e., 100% bootstrap support), even though incomplete lineage sorting or other processes that lead to genealogical discordance are at work (e.g., gCF and sCF values ≪100%).”  On their 3,220-locus lizard dataset the bootstrap is 100% on every branch but one and the smallest gene concordance factor is 4.5%.  Note that the site concordance factor's no-information baseline is 33%, not 0%, because a random site picks one of three resolutions. ↩

  44. 44.

    Hillis, D. M. et al.  “Experimental Phylogenetics: Generation of a Known Phylogeny.”  Science 255(5044):589–592.  1992.  A laboratory experiment that generated a phylogeny whose true answer was known in advance, so that methods could be scored against it.  The article uses it as its main example of methods recovering a known tree, and it carries the authors' own caveat about the tree shape they chose. ↩

    1. 44a.

      Page 590:  “We chose to construct a symmetric phylogeny with equal distances among nodes (Fig. 1).  This topology (tree shape) has proven especially amenable to accurate reconstruction in theoretical studies and may thus be regarded as a 'null' model against which other topologies may be compared.”  Reporting the 1992 success without this sentence overstates it. ↩

  45. 45.

    Hillis, David M. et al.  “Application and Accuracy of Molecular Phylogenies.”  Science 264(5159):671–677.  1994.  Page 675:  “all methods estimated an incorrect phylogeny when using high-resolution restriction fragments, but they estimated the correct phylogeny when using restriction sites.”  “The sequence data consist of 1091 base pairs across four genes of T7 ... There are only 63 variable sites across the sequences.”  With those sequences.  “only parsimony and weighted parsimony estimate the correct tree”.  while maximum likelihood, neighbour-joining, Fitch-Margoliash and UPGMA.  “each estimate a single, incorrect tree”.  A Science article on how molecular phylogenies are tested, by computer simulation and by experiments on organisms whose history is known.  It argues that the two approaches are complementary, that simulations can be biased toward the method being tested, and that known phylogenies add a reality check.  The article uses it for the same known T7 tree, where the methods gave wrong answers from restriction-fragment data and four of six were wrong on 1,091 base pairs containing only 63 variable sites. ↩

  46. 46.

    Strait, D. S., and F. E. Grine.  “Inferring hominoid and early hominid phylogeny using craniodental characters: the role of fossil taxa.”  Journal of Human Evolution 47(6):399–452.  2004.  A cladistic study (a method that builds family trees from shared traits) of one hundred and ninety-eight craniodental characters (skull and teeth features) from apes and early hominids.  It tested whether adding fossil species changes the tree: trees from living apes alone were inconsistent with the molecular ape tree, while trees that included fossil hominids matched it, which the authors took to refute the claim that these characters are unreliable.  The article cites it as one side of a dispute with Collard and Wood, in which Strait and Grine replied that adding fossils fixes the problem. ↩

  47. 47.

    Collard, Mark, and Bernard Wood.  “Hominin homoiology: an assessment of the impact of phenotypic plasticity on phylogenetic analyses of humans and their fossil relatives.”  Journal of Human Evolution 52(5):573–584.  2007.  A study that tested the homoiology hypothesis, which says that shape changes from chewing strain in the hominin skull, not shared descent, mislead trees built from it.  Using hominoid (ape) skull measurements with a robust molecular tree, the authors found that the strain-affected measurements varied more, as predicted, but the trees built from them matched the molecular tree better, so homoiologies had no significant negative effect.  The article uses it in the back-and-forth with Strait and Grine over hominin trees. ↩

  48. 48.

    Sayyari, E., and S. Mirarab.  “Fast Coalescent-Based Computation of Local Branch Support from Quartet Frequencies.”  Molecular Biology and Evolution 33(7):1654–1668.  2016.  A methods paper that gives each branch of a species tree (the tree of species, as opposed to the trees of single genes) a score for how likely it is to be correct.  The authors tested the score on simulated trees with known answers and on real datasets, and found it very precise and better than multi-locus bootstrapping, a resampling method, at recovering true branches.  The article cites it as a well calibrated support score, unlike the bootstrap values most studies still print. ↩

  49. 49.

    Penny et al.  “Testing the theory of evolution by comparing phylogenetic trees constructed from five different protein sequences.”  Nature 297(5863):197–200.  1982.  A 1982 test of evolutionary theory that builds a phylogenetic tree, a family tree of species, from each of five protein sequences and compares the trees with one another.  The article uses it as the origin of the test it describes, in which agreement between trees counts as evidence only if disagreement was possible. ↩

  50. 50.

    Robinson, and Foulds.  “Comparison of phylogenetic trees.”  Mathematical Biosciences 53(1–2):131–147.  1981.  A 1981 paper on how to compare two phylogenetic trees.  The article uses it for the Robinson-Foulds distance, which counts the groupings that are present in one tree but not in the other. ↩

  51. 51.

    Smith, Martin R.  “Information theoretic generalized Robinson-Foulds metrics for comparing phylogenetic trees.”  Bioinformatics 36(20):5007–5013.  2020.  A methods paper that proposes new measures of how far apart two phylogenetic trees are, in bits of information, and tests them against the standard Robinson-Foulds measure.  In a test on 100,000 random 11-leaf trees compared with a pectinate (caterpillar-shaped) reference tree, 86% sat at the maximum Robinson-Foulds distance, and that measure reported fewer than eight distinct distance values.  The author concludes that mutual clustering information should be preferred to Robinson-Foulds, and the article uses the test to show that the standard measure is a poor gauge of how close two trees are. ↩

  52. 52.

    Darlu, Pierre, and Guillaume Lecointre.  “When does the incongruence length difference test fail?”  Molecular Biology and Evolution 19(4):432–437.  2002.  A simulation study of the incongruence length difference (ILD) test, a standard significance test for whether two sets of characters, such as DNA sequences, conflict with each other.  Under the same evolutionary conditions the test rejected agreement less often than its nominal 5% level, but it had little power to detect real conflict unless there were many informative sites and the rate of change was similar across sites.  The article uses it to show that the standard test for conflict is weaker than its use implies. ↩

  53. 53.

    Sanderson, Michael J., and Michael J. Donoghue.  “Patterns of variation in levels of homoplasy.”  Evolution 43(8):1781–1795.  1989.  A statistical study of 60 published cladistic analyses, which are trees built from several kinds of characters, such as anatomy or molecular data.  It found that the consistency index, a score of how well the data fit a tree, falls as more taxa are added, because homoplasy (the same trait evolving independently more than once) rises.  The article uses it for that effect of taxon number on the index. ↩

  54. 54.

    Steell, Elizabeth M. et al.  “Revealing patterns of homoplasy in discrete phylogenetic datasets with a cross-comparable index.”  Zoological Journal of the Linnean Society 204(1):zlaf024.  2025.  A methods paper that tests the consistency index (CI) and retention index (RI), two widely used measures of homoplasy (the same trait evolving more than once), on simulated and real morphological datasets, and proposes a new index, the relative homoplasy index (RHI).  It measures CI, RI and RHI against taxon number on the same datasets and finds that they move in opposite directions, whereas Sanderson and Donoghue measured only the consistency index.  The article cites it for the point that the indices used in the homoplasy literature move with taxon number in opposite directions. ↩

  55. 55.

    Jarvis, Erich D. et al.  “Whole-genome analyses resolve early branches in the tree of life of modern birds.”  Science 346(6215):1320–1331.  2014.  A genome study by a large international team that built a tree of modern birds from 48 avian genomes and then checked it against each locus separately.  The main tree was strongly supported, but no single intron, exon or ultraconserved element locus had a topology identical to it, and the least certain branches were short ones at deep divergences.  The article uses it to show where the conflict between genes sits in the bird tree. ↩

    1. 55a.

      “a wall of many (25 of 45; ~55%) very short internal branches (0.0006 to 0.002 substitutions per site), almost all at deep divergences within Neoaves.”  All nine branches below 100% bootstrap were among the shortest. ↩

    2. 55b.

      Per-branch statistics: 0-38% of 2,485 intron gene trees conflict at above 75% bootstrap; only 0-20% support the same branches at that level.  The remainder is uninformative rather than conflicting, and the authors state.  “the conflicting nodes largely reflected branches with low statistical support.” ↩

  56. 56.

    Wickett, Norman J. et al.  “Phylotranscriptomic analysis of the origin and early diversification of land plants.”  PNAS 111(45):E4859–E4868.  2014.  A transcriptome study that generated data for 92 streptophyte taxa (the group of green algae and land plants) and added 11 published plant genomes, using up to 852 nuclear genes.  The authors ran 69 analyses with different data and methods to test how robust the tree was, and found robust support for land plants sitting with one group of green algae, the Zygnematophyceae.  The article cites it for the conflict percentages reported for 91 branches of its trees. ↩

  57. 57.

    Smith, Stephen A. et al.  “Analysis of phylogenomic datasets reveals conflict, concordance, and gene duplications with examples from animals and plants.”  BMC Evolutionary Biology 15:150.  2015.  A methods article describing phyparts, a program that maps onto a species tree how many gene trees agree with each branch, conflict with it, or contain gene duplications.  The authors applied it to published transcriptome datasets from ants, bees and wasps (5,863 gene groups) and from flowering plants (Caryophyllales), and found significant conflict throughout both trees.  The article cites it for the point that conflicting gene trees are the more poorly supported ones. ↩

  58. 58.

    Song, Sen et al.  “Resolving conflict in eutherian mammal phylogeny using phylogenomics and the multispecies coalescent model.”  PNAS 109(37):14942–14947.  2012.  A genome-scale study of 447 nuclear genes from 33 placental (eutherian) mammal species and four outgroups, analyzed with coalescent methods, which model gene trees that can differ from the species tree.  The authors found that joining all genes into one data set gave strongly conflicting trees when they subsampled the data, while the coalescent methods gave a stable tree.  The article cites it for the finding that the coalescent model accounts for 77% of the variation in gene trees, a prediction of common descent rather than a failure of it. ↩

  59. 59.

    Reddy, S. et al.  “Why do phylogenomic data sets yield conflicting trees? Data type influences the avian tree of life more than taxon sampling.”  Systematic Biology 66(5):857–879.  2017.  A study of the early divergences in Neoaves, the largest group of living birds, which built a new matrix of 235 species that is 90% non-coding DNA to find out why two earlier avian phylogenomic studies disagree.  The authors found that data type mattered more than the number of taxa: the new trees agreed more with the earlier study built mostly on non-coding DNA, and they blame model fit, a problem that may be difficult to detect.  The article uses it to show that the disagreement tracks the kind of data more than the number of species, so more data does not settle it. ↩

  60. 60.

    Pisani, Davide et al.  “Genomic data do not support comb jellies as the sister group to all other animals.”  PNAS 112(50):15402–15407.  2015.  A reanalysis of three genome datasets that had been used to argue that comb jellies (ctenophores) are the closest living relatives of all other animals.  The authors say the earlier studies used simple models of how proteins change and distant outgroups to root the tree, and they corrected for both, along with a gene-loss bias.  With better-fitting models they found no support for comb jellies in that position, and they conclude the earlier support came from uncorrected bias.  The article uses it as one side of the dispute over the deepest branch of the animal tree, where the choice of model decides the answer. ↩

  61. 61.

    Li, Yuanning et al.  “Rooting the animal tree of life.”  Molecular Biology and Evolution 38(10):4322–4333.  2021.  A reanalysis of whether sponges or comb jellies (ctenophores) are the closest living relatives of all other animals.  The authors gathered 15 earlier phylogenomic studies and re-ran them under the same conditions with more models.  Comb jellies came out as the closest relatives across the whole range of conditions, while sponges appeared only in narrow settings, with most outgroups removed and a complex model with hundreds of extra categories, which fit no better than simpler models.  The article uses it as an example of the deepest branch of the animal tree flipping with the model. ↩

  62. 62.

    Pease, James B. et al.  “Phylogenomics reveals three sources of adaptive variation during a rapid radiation.”  PLoS Biology 14(2):e1002379.  2016.  A genome-wide study of all 13 species of wild tomato, using whole transcriptomes (the sets of genes switched on in a plant's tissues) to reconstruct how the group diversified.  The authors date the radiation to about 2.5 million years ago and find three sources of adaptive variation: hybridisation between species after they split, new genes evolved within lineages, and ancestral variation sorted differently across populations.  The article uses it for gene exchange between wild tomato lineages, which it traces to hybridisation rather than to the sorting of ancestral variation. ↩

  63. 63.

    Shen, Xing-Xing et al.  “Contentious relationships in phylogenomic studies can be driven by a handful of genes.”  Nature Ecology & Evolution 1:0126.  2017.  A reanalysis of published phylogenomic datasets from plants, animals and fungi, measuring how much each gene and site favours each branch of the tree.  The authors examined 17 contested branches and 6 well-established control branches, and found that removing the single most influential gene flipped the best-supported answer on 3 of the contested branches.  The article uses it to show that a single gene can flip a branch only where the branch was already contested, since removing random genes changed nothing and removing the top gene from a control branch did not switch its answer. ↩

  64. 64.

    Degnan, James H., and Noah A. Rosenberg.  “Discordance of species trees with their most likely gene trees.”  PLoS Genetics 2(5):e68.  2006.  A theoretical study that uses the coalescent model, a population-genetics model of how gene lineages merge backward in time, to prove a counterintuitive result.  For any species tree with five or more species there are branch lengths at which the most likely gene tree is not the species tree, so taking the most common gene tree as the species tree can converge on a wrong answer as more genes are added.  The article uses it for the proof that this anomaly zone exists. ↩

  65. 65.

    Linkem, Charles W. et al.  “Detecting the anomaly zone in species trees and evidence for a misleading signal in higher-level skink phylogeny.”  Systematic Biology 65(3):465–477.  2016.  A phylogenomic study of 16 lizard species in the skink family (Scincidae), using 429 new loci of DNA, that tests whether a hard-to-resolve part of the tree lies in the anomaly zone.  That is a set of conditions in which a wrong gene tree is more probable than the true species tree, and the authors checked for it by estimating ancestral population sizes and how long each species persisted.  They found at least three regions of the skink tree with the pattern the anomaly zone predicts, and say this helps reconcile conflicting earlier studies.  The article uses it as a real case showing that the anomaly zone exists in practice. ↩

  66. 66.

    Tehrani, Jamshid, and Mark Collard.  “Investigating cultural evolution through biological phylogenetic analyses of Turkmen textiles.”  Journal of Anthropological Archaeology 21:443–463.  2002.  A study that applied cladistic tree-building to decorative designs on woven Turkmen textiles, to test whether cultural change works more like branching descent or like borrowing between neighbouring groups.  Before Tsarist Russia's defeat of the Turkmen in 1881 the consistency index was 0.68 and branching accounted for about 70% of the resemblances; after 1881 it was 0.61, and the index sits in the same range as published hominid cladograms (family trees of the great apes and humans), which score 0.58 to 0.71.  The article cites it as one of the two artifact studies that point in opposite directions. ↩ ↩

  67. 67.

    O'Brien, Michael J. et al.  “Cladistics Is Useful for Reconstructing Archaeological Phylogenies: Palaeoindian Points from the Southeastern United States.”  Journal of Archaeological Science 28(10):1115–1136.  2001.  Figure 11:  “The shortest tree we found ... has a length of 22 and a CI of 0.59”.  Collard et al. 2006, Table 3, gives the retention index of 0.70 for the same 17 taxa and 8 characters.  A cladistics study, a method that sorts objects into a nested family tree by the features they inherited, applied to seventeen Palaeoindian stone projectile-point classes from the southeastern United States, coded for eight characters.  The authors built the shortest tree, which needs the fewest changes, and argue that such trees are testable hypotheses about how the point styles descended from one another.  The article uses it for the consistency index of 0.59 on this dataset, and sets against it the retention index of 0.70 that Collard and colleagues report for the same data. ↩

  68. 68.

    Rogers, Deborah S., and Paul R. Ehrlich.  “Natural selection and cultural rates of change.”  PNAS 105(9):3416–3420.  2008.  A study of design features on Polynesian canoes, split into functional features, which are tested against the environment, and symbolic features, which are not.  The authors found that natural selection seems to slow the change of functional designs, while symbolic designs differentiate more rapidly, which they take to show that cultural change can follow theoretically derived patterns.  The article uses the canoe traits as its test case, in which functional design fits the true history worse than randomised data does. ↩ ↩

  69. 69.

    Measured for this article.  research/sim/calibrate.py --plan rate --seed 20260815 for the ceiling:  process S1, Mk characters on a random tree with equal branch lengths at substitution rates 0.01 to 0.05 and 1,000 characters, reads exactly 1.0000 at 12 taxa on all three seeds and at 24 taxa at rate 0.02 on all three.  Give the same tree exponentially distributed branch lengths, so that some branches are too short to resolve, and it reads 0.75 to 0.98 instead — which is what the qualification in the sentence above is doing.  research/sim/calibrate.py --plan design --seed 20260815 for the floor:  process S5, characters drawn independently per object, gives a median of −0.003 over eighteen conditions (12, 24 and 48 taxa × two character budgets × three seeds), range −0.092 to +0.091.  Raw output in research/results/calibration-summary-rate.csv and research/results/calibration-summary-design.csv. ↩

  70. 70.

    Swamidass, S. J.  “Winston Ewert: The Dependency Graph of Life.”  Peaceful Science discussion thread.  Jul 21, 2018.  An online discussion thread, not a peer-reviewed paper, in which Swamidass replied to Ewert's paper; 23 posts survive, numbered up to 35, and Ewert took part for one day and then withdrew.  Post #25 makes the scale argument that a Bayes factor must be read against the size of the model's total surprisal, not in isolation, and Swamidass quoted the UniRef-50 row with its column labels transposed, though the arithmetic and the point stand.  The article cites it as the source of the scale argument it reframes, and this thread and its 2023 sequel are the entire published response to the paper. ↩

  71. 71.

    Lammers, Fritjof et al.  “Phylogenetic Conflict in Bears Identified by Automated Discovery of Transposable Element Insertions in Low-Coverage Genomes.”  Genome Biology and Evolution 9(10):2862–2878.  2017.  A genome study of seven bear species that built a pipeline, TeddyPi, to find transposable elements (pieces of DNA that copy and move themselves) in low-coverage genomes, and used them to reconstruct the bear family tree.  The authors got a well-resolved species tree, and they trace the conflict around the Asiatic black, sun and sloth bears to incomplete lineage sorting, where ancestral gene versions have not yet sorted out between species.  The article uses it as a counterweight: this marker class is treated as nearly free of repeated change, and it still gives conflicting counts. ↩

  72. 72.

    O'Connor, Anne, and Matthew A. Wills.  “Measuring Stratigraphic Congruence Across Trees, Higher Taxa, and Time.”  Systematic Biology 65(5):792–811.  2016.  A study of how well the order of branching in a tree matches the order in which fossils first appear, using 647 published animal and plant cladograms (family trees) spanning the entire Phanerozoic eon, the span of the animal and plant fossil record.  Eleven biasing variables account for 2% of the deviance (the share of variation in the scores that they explain) in the gap excess ratio (GER), against 26% or more for the other indices, so GER is the one to use.  The article cites it for the warning that published comparisons using the stratigraphic consistency index (SCI) and the modified Manhattan stratigraphic measure (MSM*) should be distrusted. ↩

  73. 73.

    Springer, Mark S. et al.  “The Adequacy of Morphology for Reconstructing the Early History of Placental Mammals.”  Systematic Biology 56(4):673–684.  2007.  A Points of View article that tested whether anatomy alone can rebuild the early history of placental mammals.  The authors ran pseudoextinction analyses, deleting the DNA of one living order at a time and re-running the tree on morphology alone.  The article cites it as the first of three studies on how well morphology recovers the molecular tree. ↩

  74. 74.

    Springer, Mark S. et al.  “Evolutionary Models for the Diversification of Placental Mammals Across the KPg Boundary.”  Frontiers in Genetics 10:1241.  2019.  A paper that re-ran the pseudoextinction experiment on O'Leary et al. 2013's 4,541-character matrix, which is twenty-four times larger than the 2007 dataset.  The majority of orders still moved, which kills the sample-size defence of morphology.  The article cites it for the re-run on the larger matrix. ↩

    1. 74a.

      “all three orders with highly specialized myrmecophages (Xenarthra, Pholidota, Tubulidentata) cluster with one of the other myrmecophagous orders when treated as pseudoextinct.”  The failure is driven by diet, not by noise, which is why more characters do not fix it. ↩

  75. 75.

    Brady, Peggy L., and Mark S. Springer.  “The effects of fossil taxa, hypothetical predicted ancestors, and a molecular scaffold on pseudoextinction analyses of extant placental orders.”  PLoS ONE 16(9):e0257338.  2021.  “only 8 of 19 pseudoextinct placental orders (42%) retained the same interordinal placement as on the molecular scaffold.”  A study in which each placental order was treated as extinct in turn and placed by morphology, with fossils, reconstructed ancestors and a molecular scaffold (a fixed frame taken from DNA trees).  The authors' figure of 42% is the share of orders that kept the same placement with all three; without fossils or a scaffold, the figure was 21%.  The article cites it for how often morphology still gets the placement wrong, even with fossils and reconstructed ancestors. ↩ ↩

  76. 76.

    Wood, Todd Charles.  “Taxon Sample Size in Hominin Baraminology: A Response to O'Micks.”  Answers Research Journal 9:369–372.  2016.  A reply to O'Micks, who had concluded that Homo naledi was probably not part of the human kind, which re-runs O'Micks's own 87-character matrix at two taxon samples.  With the larger sample every pair within Homo correlated significantly, and with six species only one pair did, including Homo sapiens against Homo neanderthalensis, which lost significance.  The article cites it to show that the answer depends on which species are included. ↩

  77. 77.

    DeWitt, David A. et al.  “Baraminological Analysis Places Homo habilis, Homo rudolfensis, and Australopithecus sediba in the Human Holobaramin: Discussion.”  Answers Research Journal 3:153–158.  2010.  “it would appear that the statistical baraminology results placing Homo habilis, Homo rudolfensis, and Australopithecus sediba within the human holobaramin is a false positive result.”  The nine hominin analyses behind the article's sentence that creationist analysts disagree.  Au. sediba human: Wood 2010 (ARJ 3:71-90), Wood 2011 (JCTS B 1:6-14), Wood 2016 (JCTS B 6:1-8) and Sinclair & Wood 2021 (ARJ 14:451-462).  Called a false positive: DeWitt, Habermehl & Menton 2010 (ARJ 3:153-158), a reply paper that ran no baraminic distance analysis of its own.  H. naledi human: O'Micks 2016 (JCTS B 6), Wood 2016 (JCTS B 6:14-30), Sinclair & Wood 2021.  Not human: O'Micks 2016 (ARJ 9:263-272) and O'Micks 2017 (ARJ 10:103-113).  The factor-of-ten downweighting is in O'Micks's 2016 reply on taxon sample size.  A reply to Todd Wood's 2010 paper on the same fossils, which calls the result placing Homo habilis, Homo rudolfensis and Australopithecus sediba in the human kind a false positive, without running a baraminic distance analysis of its own.  The article cites it for the flat disagreement among creationist analysts on the same data. ↩

  78. 78.

    O'Micks, Jean.  “Reply to 'Taxon Sample in Hominin Baraminology: A Response to O'Micks'.”  Answers Research Journal 9:373–375.  2016.  Page 373:  “We cannot demand that God created each baramin to have significant morphological differences compared to other baramins.”  The same page reports the factor-of-ten downweighting of the craniodental characters and concedes that.  “the weighted BDIST program does not calculate bootstrap values.”  The earlier citation gave the title as.  “Reply to 'Taxon Sample Size in Hominin Baraminology: A Response to O'Micks'”.  A reply by Jean O'Micks, an independent scholar, to a critic of his baraminology study of Homo naledi, a fossil species.  Baraminology is a creationist method that sorts species into groups thought to descend from one created kind, and O'Micks argues that a weighted version of the BDIST method places H. naledi with the australopiths.  The article uses it to show a baraminologist conceding that these groups need not differ in form. ↩

  79. 79.

    Penny, David, and Michael Hendy.  “Estimating the reliability of evolutionary trees.”  Molecular Biology and Evolution 3(5):403–417.  1986.  A 1986 paper that tests how reliable an evolutionary tree is, using six real protein sequences from the same 11 mammals rather than computer simulations.  The authors measure how trees built from larger random subsets of nucleotide positions become more alike, and report that the rodent-primate-lagomorph trichotomy is the least certain part of the final tree, followed by the position of the dog.  The article takes the design and headline figures of the 1982 test from this follow-up.

  80. 80.

    Theobald, D.  “29+ Evidences for Macroevolution, Part 1: the unique universal phylogenetic tree.”  TalkOrigins Archive.  2012.  A web article in the TalkOrigins Archive, a long-running online encyclopedia on evolution and creation, that sets out predictions of common descent and includes a section on how many trees are possible for a set of species.  The article cites it for the author's counts of possible trees, which it checks against his own formula.

    1. 80a.

      “nearly 14 trillion different possible phylogenies”.  The prose figure for 12 organisms, which is wrong: his own formula gives 13,749,310,575 rooted trees, about 13.7 billion.  Table 1.3.1:  Table 1.3.1 gives 34,459,425 unrooted trees for 11 taxa, which is correct.