quartetTreeErrorProb: Bayesian posterior probability of error in 4-taxon unrooted species tree topology estimate
Description
From a gene quartet count concordance factor (qcCF), computes Bayesian posterior probabilities
of the three 4-taxon species tree topologies and the Bayesian posterior probability that the
assumed topology is incorrect, under the assumption that the counts arise from the MSC on some species tree.
Usage
quartetTreeErrorProb(obs, model = "T3")
Value
(error.prob, top.probs) where error.prob is the species tree error probability
and top.probs is a vector of the three species tree topology probabilities in the order of obs;
for model "T1" the species tree used is the one
corresponding to the first count; for model "T3" the species
tree is the one corresponding to the largest count
Arguments
obs
vector of counts for 3 topologies
model
"T3" or "T1", for the models of MAR19;textualMSCquartets describing an unspecified species
tree topology ("T3"), or the topology whose count is the first entry of obs ("T1")
Details
The Jeffreys prior is used for internal branch length, along with the uniform prior
on the resolved topology.