evoTS (version 1.0.3)

loglik.surface.GRW: 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 a General Random Walk.

Usage

loglik.surface.GRW(y, mstep.vec, vstep.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.

mstep.vec

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

vstep.vec

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

pool

indicating whether to pool variances across samples

Author

Kjetil Lysne Voje

Examples

## Generate a paleoTS objects
x <- paleoTS::sim.GRW(30)

## Fit the the model to the data.
x1<-paleoTS::opt.joint.GRW(x)

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