isotracer (version 1.1.8)

set_half_life: Set the half-life for radioactive tracers

Description

Indicating a non-zero value for half-life will add a decay to the marked portion of the tracer element. The decay constant is calculated from the half-life value as:

Usage

set_half_life(nm, hl, quiet = FALSE)

Value

A networkModel object.

Arguments

nm

A networkModel object.

hl

Half-life value, in the same time unit as the observations are (or will be) given. Setting half-life to zero is equivalent to using a stable isotope (no decay used in the model).

quiet

Boolean for verbosity.

Details

lambda_decay = log(2) / half_life

Note that for correct calculations the half-life value should be given in the same time unit (e.g. hour, day) that the time unit used for observations.

Examples

library(magrittr)
x <- new_networkModel() %>%
    set_topo("32P -> root -> leaf") %>%
    set_half_life(hl = 14.268)
x