Learn R Programming

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

paleoTS

paleoTS allows the user to analyze paleontological time-series implementing many common models that are considered by paleontologists when looking at trait changes in a species over time, including random walks, directional change, Ornstein-Uhlenbeck models of adaptation, and stasis. In addition, more complex models are available that allow for punctuated change, and for shifts in dynamics within a time-series.

Example

As a simple example, we simulate directional change, and then fit three common models used in the literature to these data:

 library(paleoTS)
 x <- sim.GRW(ns = 20, ms = 1, vs = 0.3)
 fit3models(x)
#> 
#> Comparing 3 models [n = 20, method = Joint]
#> 
#>             logL K      AICc    dAICc Akaike.wt
#> GRW    -15.10736 3  37.71472  0.00000         1
#> URW    -28.45691 2  61.61970 23.90498         0
#> Stasis -61.93319 2 128.57227 90.85754         0

Copy Link

Version

Install

install.packages('paleoTS')

Monthly Downloads

412

Version

0.5.3

License

GPL-3

Maintainer

Gene Hunt

Last Published

August 8th, 2022

Functions in paleoTS (0.5.3)

cantius_L

Time-series of the length of lower first molar for the Cantius lineage
compareModels

Compare model fits for a paleontological time-series
fit.sgs

Fit a model of trait evolution with a protracted punctuation.
dorsal.spines

Time-series of dorsal spine data from a fossil stickleback lineage
ln.paleoTS

Approximate log-transformation of time-series data
lynchD

Compute Lynch's Delta rate metric
bootSimpleComplex

Bootstrap test to see if a complex model is significantly better than a simple model.
as.paleoTSfit

Create a paleoTSfit object
opt.punc

Fit a model of trait evolution with specified punctuation(s)
plot.paleoTS

Plot a paleoTS object
read.paleoTS

Read a text-file with data from a paleontological time-series
fit9models

Fit large set of models to a time-series
pool.var

Compute a pooled variance
fit3models

Fit a set of standard evolutionary models
opt.GRW.shift

Fit random walk model with shift(s) in generating parameters
fitSimple

Fit simple models of trait evolution
opt.covTrack

Fit a model in which a trait tracks a covariate
fitMult

Fit the same simple model across multiple time-series
LRI

Log-rate, Log-interval (LRI) method of Gingerich
as.paleoTS

Make a Paleontological Time-series object
opt.joint.GRW

Fit evolutionary models using the "Joint" parameterization
sim.Stasis.RW

Simulate trait evolution with a mode shift
sim.OU

Simulate an Ornstein-Uhlenbeck time-series
sim.punc

Simulate a punctuated time-series
fitGpunc

Fit trait evolution model with punctuations estimated from the data
sim.sgs

Simulate protracted punctuation
fitModeShift

Fit model in which the mode of trait evolution shifts once
sub.paleoTS

Subsample a paleontological time-series
std.paleoTS

Convert time-series to standard deviation units
mle.GRW

Analytical ML estimator for random walk and stasis models
sim.GRW

Simulate random walk or directional time-series for trait evolution
opt.GRW

Fit evolutionary model using "AD" parameterization
sim.GRW.shift

Simulate (general) random walk with shift(s) in generating parameters
opt.joint.OU

Fit Ornstein-Uhlenbeck model using the "Joint" parameterization
sim.covTrack

Simulate trait evolution that tracks a covariate
sim.Stasis

Simulate Stasis time-series for trait evolution
test.var.het

Test for heterogeneity of variances among samples in a time-series
ESD

Compute Expected Squared Divergence (ESD) for Evolutionary Models
IC

Compute Information Criteria