The goal of paleoTS is to allow the user to simulate and fit time-series models commonly used to understand trait evolution in paleontology. Models include random walks, stasis, directional trends, OU, covariate-tracking, punctuations and more. Model fitting is done via maximum likelihood.
This is a simple example in which a time-series is generated, plotted, and then fit with three common models in paleobiology. The generating model is a general (also called biased) random walk, with a pretty strong trend parameter. Usually, this model receives just about all of the available model support with these generating parameters.
library(paleoTS)
y <- sim.GRW(ns = 40, ms = 0.3)
plot(y)
fit3models(y)
#>
#> Comparing 3 models [n = 40, method = Joint]
#>
#> logL K AICc dAICc Akaike.wt
#> GRW -26.86719 3 60.40106 0.00000 1
#> URW -37.85943 2 80.04318 19.64213 0
#> Stasis -113.33758 2 230.99949 170.59844 0
Take a look at the vignette “paleoTS_basics” for more of an introduction to this package.
paleoTS should be installed from CRAN.
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.