
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.
quartetTreeErrorProb(obs, model = "T3")
(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
vector of counts for 3 topologies
"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"
)
The Jeffreys prior is used for internal branch length, along with the uniform prior on the resolved topology.
MAR19MSCquartets
obs <- c(28,32,30)
quartetTreeErrorProb(obs,model="T1")
quartetTreeErrorProb(obs,model="T3")
Run the code above in your browser using DataLab