call_method: Call a reconstruction method

In ldsr: Linear Dynamical System Reconstruction

Description

Call a reconstruction method subroutine according to the method required

Usage

call_method(
  y,
  u,
  v,
  method,
  init,
  num.restarts,
  return.init,
  ub,
  lb,
  num.islands,
  pop.per.island,
  niter,
  tol
)

Arguments

y	Catchment output, preprocessed from data
u	Input matrix for a single-model reconstruction, or a list of input matrices for an ensemble reconstruction.
v	Same as u.
method	By default this is "EM". There are experimental methods but you should not try.
init	A list, each element is a vector of initial values for the parameters. If NULL, will be created by make_init(). See make_init for details.
num.restarts	The number of initial conditions to start the EM search; ignored if init is provided.
return.init	If TRUE, the list of initial values (init) will be returned. This can be useful if you want to reproduce the model from this one set of initial values.
ub	Upper bounds, a vector whose length is the number of parameters
lb	Lower bounds
num.islands	Number of islands (if method is GA; experimental)
pop.per.island	Initial population per island (if method is GA; experimental)
niter	Maximum number of iterations, default 1000
tol	Tolerance for likelihood convergence, default 1e-5. Note that the log-likelihood is normalized by dividing by the number of observations.

Value

The results as produced by LDS_EM_restart() when the default method (EM) is called. Other methods are experimental.

ldsr documentation built on May 4, 2020, 5:06 p.m.