powered by
Estimate initial k from heavy isotope concentration data
pdr_estimate_k0(time, n, frac_k, quiet = FALSE)
Initial estimate of k0 (consumption rate constant)
Vector of numeric time values (e.g. days); first should be zero
Observed heavy isotope (as a volume), same length as time
Fractionation: 13C consumption as a fraction of 12C consumption
Suppress output message, logical
pdr_estimate_k0(1:5, c(1, 0.9, 0.7, 0.65, 0.4), frac_k = 0.98)
Run the code above in your browser using DataLab