evoTS (version 1.0.3)

loglik.surface.decel: Calculate the log-likelihood surface for a part of parameter space

Description

Function to calculate the log-likelihood surface for a part of parameter space for an Unbiased Random Walk with an decelerated rate of evolution.

Usage

loglik.surface.decel(y, vstep.vec, r.vec, pool = TRUE)

Value

the function returns the range of parameter values that are within two log-likelihood units from the best (maximum) parameter estimate and a log-likelihood surface.

Arguments

y

an univariate paleoTS object.

vstep.vec

vector containing the parameter values of the variance parameter to be evaluated

r.vec

vector containing the parameter values of the r parameter to be evaluated

pool

indicating whether to pool variances across samples

Author

Kjetil Lysne Voje

Examples

## Generate a paleoTS objects
x <- sim.accel.decel(30, r=-0.5)

## Fit the model to the data.
x1<-opt.joint.decel(x)

# \donttest{
## Create log-likelihood surface (the example may take > 5 seconds to run)
loglik.surface.decel(x, vstep.vec = seq(0, 5, 0.1), r.vec = seq(-1, 0, 0.01))
# }