Package: MSCquartets 3.0

MSCquartets: Analyzing Gene Tree Quartets under the Multi-Species Coalescent

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) <doi:10.1214/19-EJS1576>, simplex plots of quartet concordance factors as presented by Allman et al. (2020) <doi:10.1101/2020.02.13.948083>, species tree inference methods based on quartet distances of Rhodes (2019) <doi:10.1109/TCBB.2019.2917204> and Yourdkhani and Rhodes (2019) <doi:10.1007/s11538-020-00773-4>, the NANUQ algorithm for inference of level-1 species networks of Allman et al. (2019) <doi:10.1186/s13015-019-0159-2>, the TINNIK algorithm for inference of the tree of blobs of an arbitrary network of Allman et al.(2022) <doi:10.1007/s00285-022-01838-9>, and NANUQ+ routines for resolving multifurcations in the tree of blobs to cycles as in Allman et al.(2024) (forthcoming). Software announcement by Rhodes et al. (2020) <doi:10.1093/bioinformatics/btaa868>.

Authors:Elizabeth Allman [aut], Hector Banos [aut], Jonathan Mitchell [aut], Kristina Wicke [aut], John Rhodes [aut, cre]

MSCquartets_3.0.tar.gz
MSCquartets_3.0.zip(r-4.5)MSCquartets_3.0.zip(r-4.4)MSCquartets_3.0.zip(r-4.3)
MSCquartets_3.0.tgz(r-4.4-x86_64)MSCquartets_3.0.tgz(r-4.4-arm64)MSCquartets_3.0.tgz(r-4.3-x86_64)MSCquartets_3.0.tgz(r-4.3-arm64)
MSCquartets_3.0.tar.gz(r-4.5-noble)MSCquartets_3.0.tar.gz(r-4.4-noble)
MSCquartets_3.0.tgz(r-4.4-emscripten)MSCquartets_3.0.tgz(r-4.3-emscripten)
MSCquartets.pdf |MSCquartets.html
MSCquartets/json (API)

# Install 'MSCquartets' in R:
install.packages('MSCquartets', repos = c('https://johnarhodes.r-universe.dev', 'https://cloud.r-project.org'))

Peer review:

Uses libs:
  • c++– GNU Standard C++ Library v3
Datasets:

On CRAN:

This package does not link to any Github/Gitlab/R-forge repository. No issue tracker or development information is available.

2.82 score 11 scripts 657 downloads 67 exports 27 dependencies

Last updated 24 days agofrom:6fe05b867f. Checks:OK: 9. Indexed: yes.

TargetResultDate
Doc / VignettesOKNov 19 2024
R-4.5-win-x86_64OKNov 19 2024
R-4.5-linux-x86_64OKNov 19 2024
R-4.4-win-x86_64OKNov 19 2024
R-4.4-mac-x86_64OKNov 19 2024
R-4.4-mac-aarch64OKNov 19 2024
R-4.3-win-x86_64OKNov 19 2024
R-4.3-mac-x86_64OKNov 19 2024
R-4.3-mac-aarch64OKNov 19 2024

Exports:allPermsblobDistanceBQinferencecircHybOrderscollapseEdgescombineCycleResolutionscompatibleSplitscutDensityestimateEdgeLengthsexpectedCFsexpmodNANUQCycleDistexpNANUQCycleDistfitCycleOrdersHolmBonferroniinitBquartetslabelIntNodesM0NANUQNANUQdistnexusDistnodeGroupsordersHeuristicmodNANUQpowerDivStatpvalHistQDCQDSquartetBTinferencePlotquartetCutMLEquartetCutTestquartetCutTestIndquartetDistquartetNetworkDistquartetStarTestquartetStarTestIndquartetTablequartetTableCollapsequartetTableDominantquartetTableParallelquartetTablePrintquartetTableResolvedquartetTallyCppquartetTestPlotquartetTreeErrorProbquartetTreeTestquartetTreeTestIndquartetWeightedDistresolveCycleresolveLevel1simplexCoordssimplexLabelssimplexPointsimplexPreparesimplexSegmentsimplexTextsortQuartetTableRowsT1densityT3densitytaxonNamesTINNIKTINNIKdisttopDisttreeFromSplitstreeOfBlobsWQDCWQDCrecursiveWQDSWQDSAdjustLengths

Dependencies:apeclicodetoolscpp11digestdoParallelfastmatchforeachgenericsglueigraphiteratorslatticelifecyclemagrittrMatrixnlmephangornpkgconfigquadprogrbibutilsRcppRcppProgressRdpackrlangvctrszipfR

Inferring a Level-1 network with the NANUQ+ routines

Rendered fromNANUQplus.Rmdusingknitr::rmarkdownon Nov 19 2024.

Last update: 2024-10-30
Started: 2024-10-30

Inferring a Tree of Blobs with TINNIK

Rendered fromTINNIK.Rmdusingknitr::rmarkdownon Nov 19 2024.

Last update: 2024-08-19
Started: 2024-08-19

Readme and manuals

Help Manual

Help pageTopics
Multispecies Coalescent Model Quartet PackageMSCquartets-package MSCquartets
Generate permutationsallPerms
Compute empirical distance between taxon groups.blobDistance
Main loop of B-quartet inferenceBQinference
Generate all circular orders with designated hybridcircHybOrders
Collapse short tree edgescollapseEdges
Combine several cycle resolutions on a tree of blobs to create a networkcombineCycleResolutions
Extract compatible splitscompatibleSplits
Probability density function for Cut ModelcutDensity
Simulated gene tree dataset from species treedataGeneTreeSample
Papionini gene tree datasetdataPapioniniVanderpool
Yeast gene tree datasetdataYeastRokas
Estimate edge lengths on a species tree from gene tree quartet countsestimateEdgeLengths
Produce table of expected quartet concordance factors for a species treeexpectedCFs
Expected modNANUQ cycle distanceexpmodNANUQCycleDist
Expected NANUQ cycle distanceexpNANUQCycleDist
Compute fit of circular orders to distance with least squaresfitCycleOrders
Apply Holm-Bonferroni method to adjust for multiple testsHolmBonferroni
Initialize vector of B quartetsinitBquartets
Label internal nodes on treelabelIntNodes
Modified Struve functionM0
Apply NANUQ network inference algorithm to gene tree dataNANUQ
Compute NANUQ distance and write to fileNANUQdist
Write a distance table to a file in nexus formatnexusDist
Groups taxa by deleting a node in a treenodeGroups
Choose cycle orders heuristically from empirical modNANUQ distanceordersHeuristicmodNANUQ
Power divergence statistic of Cressie & ReadpowerDivStat
pTable for Yeast datasetpTableYeastRokas
Plot histogram of log p-values in tablepvalHist
Compute Quartet Distance Consensus tree from gene tree dataQDC
Compute Quartet Distance SupertreeQDS
Produce simplex plot with results of B/T-quartet inferencequartetBTinferencePlot
Maximum likelihood estimate of quartet tree of blobs topology and CF under Cut modelquartetCutMLE
Hypothesis test for quartet counts fitting a resolved quartet tree of blobs under NMSCquartetCutTest
Multiple independent hypothesis tests for quartet counts fitting the Cut model under the NMSCquartetCutTestInd
Compute quartet distance between taxaquartetDist
Compute network quartet distance between taxaquartetNetworkDist
Hypothesis test for quartet counts fitting a star tree under the MSCquartetStarTest
Multiple independent hypothesis tests for gene quartet counts fitting a species quartet star tree under the MSCquartetStarTestInd
Produce table of counts of quartets displayed on treesquartetTable
Reduce quartet table by combining some taxaquartetTableCollapse
Produce table of dominant quartets, with estimates of internal edge lengthsquartetTableDominant
Produce table of counts of quartets displayed on trees, in parallel for large data setsquartetTableParallel
Print a quartet table with nice formattingquartetTablePrint
Modify quartet table to show only resolved quartetsquartetTableResolved
Build quartet table from distancesquartetTallyCpp
Produce simplex plot with results of quartet hypothesis testsquartetTestPlot
Bayesian posterior probability of error in 4-taxon unrooted species tree topology estimatequartetTreeErrorProb
Hypothesis test for quartet counts fitting a tree under the MSCquartetTreeTest
Multiple independent hypothesis tests for quartet counts fitting a species tree under the MSCquartetTreeTestInd
Compute the Weighted Quartet Distance between taxaquartetWeightedDist
Resolve a node on a tree of blobs as a cycleresolveCycle
Resolve Tree of Blobs to Level-1 networkresolveLevel1
Convert 3-d coordinates to 2-d probability simplex coordinatessimplexCoords
Label vertices of 2-d probability simplexsimplexLabels
Plot point in 2-d probability simplexsimplexPoint
Draw 2-d probability simplex, with model lines for T3 or T1 modelsimplexPrepare
Plot line segment in 2-d probability simplexsimplexSegment
Add text at a point in 2-d probability simplexsimplexText
Sort quartet table rows by lex ordersortQuartetTableRows
Probability density function for Model T1T1density
Probability density function for Model T3T3density
Quartet table for Heliconius gene tree datasettableHeleconiusMartin tableHeliconiusMartin
Quartet table for Leopardus datasettableLeopardusLescroart
Get all taxon names from a collection of treestaxonNames
TINNIK algorithm to infer species tree of blobsTINNIK
Compute TINNIK distance from quartets and hypothesis test p-valuesTINNIKdist
Topological distances on a treetopDist
Produce tree from compatible splitstreeFromSplits
Tree of blobs for a networktreeOfBlobs
Compute Weighted Quartet Distance Consensus tree from gene tree dataWQDC
Compute the Recursive Weighted Quartet Distance Consensus tree from gene tree dataWQDCrecursive
Compute the Weighted Quartet Distance SupertreeWQDS
Adjust edge lengths on tree built from Weighted Quartet distance to estimate metric treeWQDSAdjustLengths