Learn R Programming

RPANDA

R: Phylogenetic ANalyses of DiversificAtion

Implements fits of diversification and phenotypic evolution models to phylogenetic data. See Morlon et al. PLoSB (2010), Morlon et al. PNAS (2011), Condamine et al. Eco Lett (2013), Morlon et al. Eco Lett (2014), Manceau et al. Eco Lett (2015), Lewitus & Morlon Syst Biol (2016), Drury et al. Syst Biol (2016), Manceau et al. Syst Biol (2016), Clavel & Morlon PNAS (2017), Drury et al. Syst Biol (2017), Lewitus & Morlon Syst Biol (2017), Drury et al. PLoSB (2018), Clavel et al. Syst Biol (2019), Maliet et al. Nature Ecol Evol (2019), Condamine et al. Eco Lett (2019), Aristide & Morlon Eco Lett (2019), Billaud et al. Syst Biol (2019), Lewitus & Morlon Syst Biol (2019), Maliet et al. Eco Lett (2020), Drury et al. PloS B (2021), Perez-Lamarque & Morlon Mol Ecol (2022), Perez-Lamarque et al. Peer Comm J (2022), Mazet et al. Methods Ecol Evol (2023), Drury et al. Current Bio (2024).

More information on the RPANDA package and worked examples can be found in Morlon et al. (2016)

The current stable version of the RPANDA package (2.2) is available on the CRAN repository. https://cran.r-project.org/package=RPANDA

Package Installation

From gitHub You can install RPANDA directly from gitHub with devtools:

library(devtools)

install_github("hmorlon/PANDA", dependencies = TRUE)

From the binaries

You can download the pre-released binaries for Windows and Mac OS X from the release page

From the source

Otherwise, you can install it directly from the source. Download the RPANDA folder and then, from the terminal console (linux, windows or mac):

R CMD build RPANDA

This will produce the RPANDA tarball.

Then, for compiling the binary:

R CMD INSTALL --build RPANDA_2.X.tar.gz

Report an issue

Any bugs encountered when using the package can be reported here

Copy Link

Version

Install

install.packages('RPANDA')

Monthly Downloads

496

Version

2.4

License

GPL-2

Issues

Pull Requests

Stars

Forks

Maintainer

H<c3><a9>l<c3><a8>ne Morlon

Last Published

January 14th, 2025

Functions in RPANDA (2.4)

InfTemp

Paleotemperature data across the Cenozoic
MPhiFFT

A class used internally to compute ClaDS's likelihood
JSDtree

Jensen-Shannon distance between phylogenies
JSDt_cluster

Clustering on the Jensen-Shannon distance between phylogenetic trait data
GIC.fit_pl.rpanda

Generalized Information Criterion (GIC) to compare models fit by Maximum Likelihood (ML) or Penalized Likelihood (PL).
JSDtree_cluster

Clustering of phylogenies
CreateGeoObject_BioGeoBEARS

Create biogeography object using a stochastic map from BioGeoBEARS
CreateGeobyClassObject

Create merged biogeography-by-class object
CreateClassObject

Create class object
CreateGeoObject

Create biogeography object
Phocoenidae

Phocoenidae phylogeny
Phyllostomidae

Phyllostomidae phylogeny
PhenotypicBM-class

Class "PhenotypicBM"
PhenotypicACDC-class

Class "PhenotypicACDC"
PhenotypicGMM-class

Class "PhenotypicGMM"
PhenotypicModel-class

Class "PhenotypicModel"
PhenotypicDD-class

Class "PhenotypicDD"
add.gts

Geological time scale
apply_prob_dtt

Calculates paleodiversity dynamics with the probabilistic approach.
RPANDA-package

RPANDA
Posdef

Positive definite symmetric matrices
ancestral

Estimation of traits ancestral states.
PhenotypicPM-class

Class "PhenotypicPM"
coccolithophore

Coccolithophore diversity since the Jurassic
PhenotypicOU-class

Class "PhenotypicOU"
co2_res

co2 data since the beginning of the Cenozoic
Phyllostomidae_genera

Phylogenies of Phyllostomidae genera
build_network.BipartiteEvol

Build the interaction network in BipartiteEvol
createModel

Creation of a PhenotypicModel
fitTipData

Maximum likelihood estimators of a model's parameters
d13c

d13c data since the Jurassic
PhenotypicADiag-class

Class "PhenotypicADiag"
co2

co2 data since the Jurassic
createModelCoevolution

Creation of a PhenotypicGMM
fit_ClaDS

Fit ClaDS to a phylogeny
define_species.BipartiteEvol

Build the phylogenies for BipartiteEvol
fitTipData-methods

~~ Methods for Function fitTipData ~~
div.rates

Diversification rates through time
div.models

Model comparison of diversification models
delineate_phylotypes

Automatic phylotypes delineation
fit_coal_var

Fit birth-death model using a coalescent approch
fit_t_comp

Fits models of trait evolution incorporating competitive interactions
fit_bd_backbone_c

Maximum likelihood fit of the general birth-death model (backbone and constraints)
fit_bd

Maximum likelihood fit of the general birth-death model
fit_env_in_past

Maximum likelihood fit of the environmental birth-death model excluding the recent past
fit_bd_backbone

Maximum likelihood fit of the general birth-death model (backbone)
fit_ClaDS0

Infer ClaDS0's parameter on a phylogeny
fit_bd_in_past

Maximum likelihood fit of the general birth-death model excluding the recent past
fit_coal_cst

Maximum likelihood fit of the equilibrium model
fit_env

Maximum likelihood fit of the environmental birth-death model
fit_t_env_ou

Maximum likelihood fit of the OU environmental model of trait evolution
get.comb.shift

Combinations of shifts of diversification.
getDataLikelihood-methods

~~ Methods for Function getDataLikelihood ~~
fit_t_comp_subgroup

Fits models of trait evolution incorporating competitive interactions, restricting competition to occur only between members of a subgroup
fit_t_standard

Fits standard models of trait evolution incorporating known and nuisance measurement error
fit_t_env

Maximum likelihood fit of the environmental model of trait evolution
fit_t_pl

High-dimensional phylogenetic models of trait evolution
getDataLikelihood

Likelihood of tip trait values.
get.sampling.fractions

Sampling fractions of subclades
foraminifera

Foraminifera diversity since the Jurassic
gic_criterion

Generalized Information Criterion (GIC) to compare models fit by Maximum Likelihood (ML) or Penalized Likelihood (PL).
greenalgae

Green algae diversity since the Jurassic
getMAPS_ClaDS

Gets the Maximum A Posteriori for each ClaDS parameter
getTipDistribution-methods

Distribution of tip trait values.
getTipDistribution

Distribution of tip trait values.
likelihood_coal_cst

Likelihood of a phylogeny under the equilibrium diversity model
likelihood_bd_backbone

Likelihood of a phylogeny under the general birth-death model (backbone)
getMAPS_ClaDS0

Gets the Maximum A Posteriori for each ClaDS0 parameter
likelihood_bd

Likelihood of a phylogeny under the general birth-death model
landplant

Land plant diversity since the Jurassic
likelihood_t_env

Likelihood of a dataset under environmental models of trait evolution.
lines.fit_t.env.ou

Add to a plot line segments joining the phenotypic evolutionary optimum through time estimated by the fit_t_env_ou function
likelihood_coal_var

Likelihood of a birth-death model using a coalescent approch
lines.fit_t.env

Add to a plot line segments joining the phenotypic evolutionary rate through time estimated by the fit_t_env function
likelihood_t_MC

Likelihood of a dataset under the matching competition model.
likelihood_t_DD

Likelihood of a dataset under diversity-dependent models.
likelihood_t_DD_geog

Likelihood of a dataset under diversity-dependent models with biogeography.
make_gen.BipartiteEvol

Compute the genealogies for BipartiteEvol
likelihood_subgroup_model

Likelihood of a dataset under models with biogeography fit to a subgroup.
modelSelection-methods

~~ Methods for Function modelSelection ~~
mantel_test

Compute Mantel test
modelSelection

Phenotypic model selection from tip trait data.
mycorrhizal_network

Mycorrhizal network from La Réunion island
plot.fit_t.env.ou

Plot the phenotypic evolutionary optimum through time estimated by the fit_t_env_ou function
mantel_test_nbpartners

Compute Mantel test
paleodiv

Paleodiversity through time
phyl.pca_pl

Regularized Phylogenetic Principal Component Analysis (PCA).
likelihood_t_MC_geog

Likelihood of a dataset under the matching competition model with biogeography.
plot_ClaDS0_chains

Plot the MCMC chains obtained when infering ClaDS0 parameters
ostracoda

Ostracod diversity since the Jurassic
phylosignal_sub_network

Compute clade-specific phylogenetic signals in a bipartite interaction network
plot_ClaDS_chains

Plot the MCMC chains obtained when infering ClaDS parameters
phylosignal_network

Compute phylogenetic signal in a bipartite interaction network
pi_estimator

Compute nucleotidic diversity (Pi estimator)
plot_dtt

Plot diversity through time
plot_spectR

Spectral density plot of a phylogeny.
plot_prob_dtt

Plot diversity through time with confidence intervals.
plot_phylosignal_sub_network

Plot clade-specific phylogenetic signals in a bipartite interaction network
plot_phylo_comb

Plot shifts of diversifcation on a phylogeny
plot_BICompare

Display modalities on a phylogeny.
plot.fit_t.env

Plot the phenotypic evolutionary rate through time estimated by the fit_t_env function
plot_fit_env

Plot speciation, extinction & net diversification rate functions of a fitted environmental model
plot_ClaDS_phylo

Plot a phylogeny with branch-specific values
plot_net.BipartiteEvol

Plot the output of BipartiteEvol
remove.model

Removing a model from shift.estimates output
redalgae

Red algae diversity since the Jurassic
plot_fit_bd

Plot speciation, extinction & net diversification rate functions of a fitted model
plot_div.BipartiteEvol

Plot the output of BipartiteEvol
prob_dtt

Confidence intervals of diversity through time
radiolaria

Radiolaria diversity since the Jurassic
sim_MCBD

Simulation of macroevolutionary diversification under the integrated model described in Aristide & Morlon 2019
shifts_cetacea

Cetacean shift.estimates results
sim.BipartiteEvol

Simulation of the BipartiteEvol model
sim.convergence.geo

Simulation of trait data under the model of convergent character displacement described in Drury et al. 2017
silica

Silica data across the Cenozoic
sim.divergence.geo

Simulation of trait data under the model of divergent character displacement described in Drury et al. 2017
sim_ClaDS

Simulation of the ClaDS model
sealevel

Sea level data since the Jurassic
shift.estimates

Estimating clade-shifts of diversification
sim_env_bd

Simulate birth-death tree dependent on an environmental curve
simul.comb.shift

Simulating trees from shift.estimates() results to test model adequacy
sim_t_env

Recursive simulation (root-to-tip) of the environmental model
simulateTipData

Tip trait simulation under a model of phenotypic evolution.
sim_t_tworegime

Recursive simulation (root-to-tip) of two-regime models
simulateTipData-methods

~~ Methods for Function simulateTipData ~~
sim_t_env_ou

Recursive simulation (root-to-tip) of the OU environmental model
sim_t_comp

Recursive simulation (root-to-tip) of competition models
sim_sgd

Algorithm for simulating a phylogenetic tree under the SGD model
theta_estimator

Compute Watterson genetic diversity (Theta estimator)
taxo_cetacea

Cetacean taxonomy
spectR_t

Spectral density plot of phylogenetic trait data
spectR

Spectral density plot of a phylogeny
ClaDS0_example

An example run of ClaDS0.
BICompare

Identify modalities in a phylogeny
Calomys

Calomys phylogeny
Caprimulgidae_ClaDS2

An example run of ClaDS2.
BGB.examples

BioGeoBEARS stochastic maps
Anolis.data

Anolis dataset
Caprimulgidae

The Caprimulgidae phylogeny.
Balaenopteridae

Balaenopteridae phylogeny
Cetacea

Cetacean phylogeny
Cetacea_clades

Stochastic map of clade membership in Cetacean phylogeny