Fit all univariate models to an evolutionary sequence (time-series).
Source:R/fit.all.univariate.R
fit.all.univariate.RdWrapper function to find maximum likelihood solutions for all univariate models (excluding models with mode shifts) to an evolutionary sequence (time-series).
Value
The function returns a list of all investigated models and their highest log-likelihood (and their corresponding AICc and AICc weight).
References
Hunt, G. 2006. Fitting and comparing models of phyletic evolution: random walks and beyond. Paleobiology 32:578–601
Hunt, G., Bell, M. A. & Travis, M. P. Evolution towards a new adaptive optimum: Phenotypic evolution in a fossil stickleback lineage. Evolution 62:700–710 (2008)
Examples
## ##Generate a paleoTS object.
x <- paleoTS::sim.GRW(30)
## Fit univariate models to the data.
fit.all.univariate(x, pool = TRUE)
#>
#> Comparing 9 models [n = 30, method = Joint]
#>
#> logL K AICc
#> GRW -18.88096 3 44.68500
#> URW -19.30267 2 43.04979
#> Stasis -24.48833 2 53.42110
#> StrictStasis -72.51084 1 147.16454
#> Decel -19.30267 3 45.52842
#> Accel -18.31547 3 43.55402
#> OU -16.19082 4 41.98164
#> OU model with moving optimum (ancestral state at optimum) -16.73012 4 43.06024
#> OU model with moving optimum -16.19082 5 44.88164
#> dAICc Akaike.wt
#> GRW 2.703357 0.079
#> URW 1.068146 0.178
#> Stasis 11.439460 0.001
#> StrictStasis 105.182897 0.000
#> Decel 3.546778 0.052
#> Accel 1.572373 0.138
#> OU 0.000000 0.304
#> OU model with moving optimum (ancestral state at optimum) 1.078597 0.177
#> OU model with moving optimum 2.900002 0.071