\donttest{
ricker() |>
spect_objfun(
est=c("r","sigma","N_0"),
partrans=parameter_trans(log=c("r","sigma","N_0")),
paramnames=c("r","sigma","N_0"),
kernel.width=3,
nsim=100,
seed=5069977
) -> f
f(log(c(20,0.3,10)))
f |> spect() |> plot()
if (require(subplex)) {
subplex(fn=f,par=log(c(20,0.3,10)),control=list(reltol=1e-5)) -> out
} else {
optim(fn=f,par=log(c(20,0.3,10)),control=list(reltol=1e-5)) -> out
}
f(out$par)
f |> summary()
f |> spect() |> plot()
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.