Skip to contents

Simulate a lnK ~ lnQ relationship in the Kb format

Source: R/metab_sim.R
sim_Kb.Rd

Uses linear interpolation among "nodes" (lnQ, lnK points) to describe the lnK ~ lnQ relationship

Usage

sim_Kb(
  K600_lnQ_nodes_centers,
  K600_lnQ_cnode_meanlog,
  K600_lnQ_cnode_sdlog,
  K600_lnQ_nodediffs_meanlog,
  K600_lnQ_nodediffs_sdlog,
  ...
)

Arguments

K600_lnQ_nodes_centers

Data configuration argument for pool_K600='binned'. numeric vector giving the natural-log-space centers of the discharge bins. See also calc_bins().

K600_lnQ_cnode_meanlog

For a sim model with pool_K600='binned'. The mean of a lognormal distribution describing the y=K600 value of the middle (or just past middle) node in the piecewise lnK ~ lnQ relationship.

K600_lnQ_cnode_sdlog

For a sim model with pool_K600='binned'. The sd of a lognormal distribution describing the y=K600 value of the middle (or just past middle) node in the piecewise lnK ~ lnQ relationship.

K600_lnQ_nodediffs_meanlog

For a sim model with pool_K600='binned'. The average (in log space) difference between ln(K) values of successive nodes. A non-zero value introduces a trend in K ~ Q.

K600_lnQ_nodediffs_sdlog

Hyperparameter for pool_K600='binned'. The standard deviations of the differences in estimated K600 between successive lnQ_nodes (bins), where the means of those differences are always zero.

...

Additional simulation parameters, which are ignored.

Value

A numeric vector of simulated ln(K600) values at the supplied nodes.

Examples

set.seed(42)
sim_Kb(
  K600_lnQ_nodes_centers = c(-1, 0, 1),
  K600_lnQ_cnode_meanlog = log(10),
  K600_lnQ_cnode_sdlog = 0.1,
  K600_lnQ_nodediffs_meanlog = 0.2,
  K600_lnQ_nodediffs_sdlog = 0.05
)
#> [1] 2.211446 2.439681 2.657837