Learn R Programming

⚠️There's a newer version (3.2) of this package.Take me there.

MSCquartets (version 2.0)

Analyzing Gene Tree Quartets under the Multi-Species Coalescent

Description

Methods for analyzing and using quartets displayed on a collection of gene trees, primarily to make inferences about the species tree or network under the multi-species coalescent model. These include quartet hypothesis tests for the model, as developed by Mitchell et al. (2019) , simplex plots of quartet concordance factors as presented by Allman et al. (2020) , species tree inference methods based on quartet distances of Rhodes (2019) and Yourdkhani and Rhodes (2019) , the NANUQ algorithm for inference of level-1 species networks of Allman et al. (2019) , and the TINNIK algorithm for inference of the tree of blobs of an arbitrary network of Allman et al.(2022) . Software announcement by Rhodes et al. (2020) .

Copy Link

Version

Install

install.packages('MSCquartets')

Monthly Downloads

665

Version

2.0

License

MIT + file LICENSE

Maintainer

John Rhodes

Last Published

March 23rd, 2024

Functions in MSCquartets (2.0)

WQDSAdjustLengths

Adjust edge lengths on tree built from Weighted Quartet distance to estimate metric tree
WQDCrecursive

Compute the Recursive Weighted Quartet Distance Consensus tree from gene tree data
compatibleSplits

Extract compatible splits
collapseEdges

Collapse short tree edges
TINNIKdist

Compute TINNIK distance from quartets and hypothesis test p-values
estimateEdgeLengths

Estimate edge lengths on a species tree from gene tree quartet counts
pTableYeastRokas

pTable for Yeast dataset
nodeGroups

Groups taxa by deleting a node in a tree
powerDivStat

Power divergence statistic of Cressie & Read
pvalHist

Plot histogram of log p-values in table
expectedCFs

Produce table of expected quartet concordance factors for a species tree
quartetTableCollapse

Reduce quartet table by combining some taxa
quartetTable

Produce table of counts of quartets displayed on trees
dataHeliconiusMartin

Heliconius gene tree dataset
dataYeastRokas

Yeast gene tree dataset
quartetStarTestInd

Multiple independent hypothesis tests for gene quartet counts fitting a species quartet star tree under the MSC
quartetStarTest

Hypothesis test for quartet counts fitting a star tree under the MSC
quartetCutTest

Hypothesis test for quartet counts fitting a resolved quartet tree of blobs under NMSC
quartetCutTestInd

Multiple independent hypothesis tests for quartet counts fitting the Cut model under the NMSC
quartetCutMLE

Maximum likelihood estimate of quartet tree of blobs topology and CF under Cut model
quartetTableDominant

Produce table of dominant quartets, with estimates of internal edge lengths
quartetTreeTestInd

Multiple independent hypothesis tests for quartet counts fitting a species tree under the MSC
quartetBTinferencePlot

Produce simplex plot with results of B/T-quartet inference
quartetTallyCpp

Build quartet table from distances
quartetTreeErrorProb

Bayesian posterior probability of error in 4-taxon unrooted species tree topology estimate
quartetDist

Compute quartet distance between taxa
nexusDist

Write a distance table to a file in nexus format
initBquartets

Initialize vector of B quartets
quartetNetworkDist

Compute network quartet distance between taxa
quartetTableParallel

Produce table of counts of quartets displayed on trees, in parallel for large data sets
quartetTablePrint

Print a quartet table with nice formatting
quartetTestPlot

Produce simplex plot with results of quartet hypothesis tests
quartetTreeTest

Hypothesis test for quartet counts fitting a tree under the MSC
quartetTableResolved

Modify quartet table to show only resolved quartets
simplexPrepare

Draw 2-d probability simplex, with model lines for T3 or T1 model
simplexPoint

Plot point in 2-d probability simplex
treeFromSplits

Produce tree from compatible splits
topDist

Topological distances on a tree
sortQuartetTableRows

Sort quartet table rows by lex order
taxonNames

Get all taxon names from a collection of trees
simplexLabels

Label vertices of 2-d probability simplex
simplexCoords

Convert 3-d coordinates to 2-d probability simplex coordinates
quartetWeightedDist

Compute the Weighted Quartet Distance between taxa
treeOfBlobs

Tree of blobs for a network
simplexText

Add text at a point in 2-d probability simplex
simplexSegment

Plot line segment in 2-d probability simplex
QDS

Compute Quartet Distance Supertree
T3density

Probability density function for Model T3
M0

Modified Struve function
QDC

Compute Quartet Distance Consensus tree from gene tree data
BQinference

Main loop of B-quartet inference
MSCquartets-package

Multispecies Coalescent Model Quartet Package
NANUQdist

Compute NANUQ distance and write to file
T1density

Probability density function for Model T1
NANUQ

Apply NANUQ network inference algorithm to gene tree data
HolmBonferroni

Apply Holm-Bonferroni method to adjust for multiple tests
WQDS

Compute the Weighted Quartet Distance Supertree
dataGeneTreeSample

Simulated gene tree dataset from species tree
TINNIK

TINNIK algorithm to infer species tree of blobs
cutDensity

Probability density function for Cut Model
WQDC

Compute Weighted Quartet Distance Consensus tree from gene tree data