Usage
pars.limits(phi = c(lower = 0, upper = +Inf),
sigmasq = c(lower = 0, upper = +Inf),
nugget.rel = c(lower = 0, upper = +Inf),
kappa = c(lower = 0, upper = +Inf),
kappa2 = c(lower = 0, upper = +Inf),
lambda = c(lower = -3, upper = 3),
psiR = c(lower = 1, upper = +Inf),
psiA = c(lower = 0, upper = 2 * pi),
tausq.rel = nugget.rel)