set_half_life: Set the half-life for radioactive tracers

In isotracer: Isotopic Tracer Analysis Using MCMC

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)

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.

Value

A networkModel object.

Examples

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

isotracer documentation built on Sept. 22, 2023, 1:07 a.m.