profile likelihood

profile likelihood

prof.lik(
  Res,
  a = Res$pars$a,
  b = Res$pars$b,
  sd = Res$pars$sd,
  rho = ifelse(Res$input$out.AR, 0, Res$pars$rho)
)

Arguments

Res

再生産関係（fit.SR()）の結果オブジェクト