get_init_values: Extract parameter initial values from parameter information

In SticsRPacks/CroptimizR: A Package to Estimate Parameters of Crop Models

Extract parameter initial values from parameter information

Description

Extract parameter initial values from parameter information

Usage

get_init_values(param_info)

Arguments

param_info

Information on the parameters to estimate. Either a list containing: ub and lb, named vectors of upper and lower bounds (-Inf and Inf can be used if init_values is provided), init_values, a data.frame containing initial values to test for the parameters (optional, if not provided, or if less values than number of repetitions of the minimization are provided, the, or part of the, initial values will be randomly generated using LHS sampling within parameter bounds). or a named list containing for each parameter: sit_list, list the groups of situations for which the current estimated parameter must take different values (see here for an example), ub and lb, vectors of upper and lower bounds (one value per group), init_values, the list of initial values per group (data.frame, one column per group, optional).

Value

A dataframe containing initial values for the different parameters to estimated (one column per parameter)

Examples

# Simple cases
param_info <- list(
  init_values = c(dlaimax = 0.001, durvieF = 200),
  lb = c(dlaimax = 0.0001, durvieF = 50),
  ub = c(dlaimax = 0.01, durvieF = 400)
)
CroptimizR:::get_init_values(param_info)

param_info <- list(
  init_values = data.frame(dlaimax = c(0.001, 0.002), durvieF = c(50, 200)),
  lb = c(dlaimax = 0.0001, durvieF = 50),
  ub = c(dlaimax = 0.01, durvieF = 400)
)
CroptimizR:::get_init_values(param_info)

# Cases with groups of situations per parameter
param_info <- list()
param_info$dlaimax <- list(
  sit_list = list(c(
    "bou99t3", "bou00t3", "bou99t1", "bou00t1",
    "bo96iN+", "lu96iN+", "lu96iN6", "lu97iN+"
  )),
  init_values = 0.001, lb = 0.0001, ub = 0.1
)
param_info$durvieF <- list(
  sit_list = list(
    c("bo96iN+", "lu96iN+", "lu96iN6", "lu97iN+"),
    c("bou99t3", "bou00t3", "bou99t1", "bou00t1")
  ),
  init_values = c(200, 300), lb = 50, ub = 400
)
CroptimizR:::get_init_values(param_info)

param_info <- list()
param_info$dlaimax <- list(
  sit_list = list(c(
    "bou99t3", "bou00t3", "bou99t1", "bou00t1",
    "bo96iN+", "lu96iN+", "lu96iN6", "lu97iN+"
  )),
  init_values = c(0.001, 0.002), lb = 0.0001, ub = 0.1
)
param_info$durvieF <- list(
  sit_list = list(
    c("bo96iN+", "lu96iN+", "lu96iN6", "lu97iN+"),
    c("bou99t3", "bou00t3", "bou99t1", "bou00t1")
  ),
  init_values = data.frame(c(200, 300), c(250, 350)), lb = 50, ub = 400
)
CroptimizR:::get_init_values(param_info)

SticsRPacks/CroptimizR documentation built on Dec. 16, 2024, 11:54 a.m.