std.paleoTS: Convert time-series to standard deviation units

std.paleoTSR Documentation

Convert time-series to standard deviation units


Convert time-series to standard deviation units


std.paleoTS(y, center = c("mean", "start"))



a paleoTS object


optional translation of time-series according to "mean" or "start"; see Details


The standardization expresses each sample mean as the deviation from the overall mean, divided by the pooled within-sample variance. Sample variances are also divided by the pooled sample variance.

Essentially, this converts paleontological time-series data into standard deviation units, similar to the computation of evolutionary rates in haldanes. This operation does not change the relative fit of models, but it does facilitate the comparison of parameter estimates across time-series of traits measured in different units.

If argument center = "start" the time-series is translated such that the trait mean of the first sample is zero.


the converted paleoTS object


x <- sim.Stasis(ns = 8, theta = 1, omega = 4, vp = 2)
xs <- std.paleoTS(x, center = "start")
plot(x, ylim = range(c(x$mm, xs$mm)))
plot(xs, col = "red", add = TRUE)
legend(x = "topright", c("unstandardized", "standardized"), lty=1, col=c("black", "red"), cex=0.7)

paleoTS documentation built on Aug. 9, 2022, 1:06 a.m.