
Last chance! 50% off unlimited learning
Sale ends in
simulate_error(
ph_out,
chamber_pars,
n = 1L,
use_tealeaves = ("T_air" %in% colnames(ph_out))
)
A data frame with n * nrow(ph_out)
rows. It contains all the
original output in ph_out
as well as a column .rep
indicating replicate
number from 1 to n
. Other new columns are assumed or measured chamber
parameters and 'measured' values estimated from synthetic data with
measurement error:
column name | assumed or derived? | description |
flow | assumed | chamber flow rate |
leaf_area | assumed | leaf area in chamber |
sigma_CO2_r | assumed | standard deviation of measurement error in CO2_r |
sigma_CO2_s | assumed | standard deviation of measurement error in CO2_s |
sigma_H2O_r | assumed | standard deviation of measurement error in H2O_r |
sigma_H2O_s | assumed | standard deviation of measurement error in H2O_s |
c_0 | derived | CO |
w_i | derived | Water vapor concentration within leaf [mmol / mol] |
w_a | derived | Water vapor concentration in chamber [mmol / mol] |
w_0 | derived | Water vapor concentration before entering chamber [mmol / mol] |
g_tw | derived | Leaf conductance to water vapor [mol/m |
E_area | derived | Evaporation rate per area [mmol/m |
E | derived | Total evaporation rate [mmol/s] |
CO2_r | derived | CO |
CO2_s | derived | CO |
H2O_s | derived | Water vapor concentration in chamber with measurement error [mmol / mol] |
H2O_r | derived | Water vapor concentration before entering chamber with measurement error [mmol / mol] |
E_meas | derived | Total evaporation rate (measured) [mmol/s] |
E_area_meas | derived | Evaporation rate per area (measured) [mmol/m |
g_tw_meas | derived | Leaf conductance to water vapor (measured) [mol/m |
g_sc_meas | derived | Stomatal conductance to CO |
g_tc_meas | derived | Leaf conductance to CO |
A_meas | derived | Net photosynthetic CO |
C_i | derived | Intercellular CO |
A data frame of output from photo()
or photosynthesis()
with units.
A data frame with a single row of chamber parameters. See Note below for table of required parameters.
Integer. Number of replicated simulations per row of ph_out
.
Flag. The tealeaves package uses a slightly
different equation to calculate the saturating water content of air as a
function temperature and pressure than LI-COR. If FALSE, the function uses
LI-COR's equation in the LI6800 manual. If TRUE, it uses the tealeaves
function for internal consistency. The function attempts to guess whether
ph_out
was run with tealeaves, but this can be manually overridden by
providing a value for the argument.
library(photosynthesis)
# Use photosynthesis() to simulate 'real' values
# `replace = ...` sets parameters to meet assumptions of `simulate_error()`
lp = make_leafpar(replace = list(
g_sc = set_units(0.1, mol/m^2/s),
g_uc = set_units(0, mol/m^2/s),
k_mc = set_units(0, 1),
k_sc = set_units(0, 1),
k_uc = set_units(0, 1)
),
use_tealeaves = FALSE)
ep = make_enviropar(replace = list(
wind = set_units(Inf, m/s)
), use_tealeaves = FALSE)
bp = make_bakepar()
cs = make_constants(use_tealeaves = FALSE)
chamber_pars = data.frame(
flow = set_units(600, umol / s),
leaf_area = set_units(6, cm ^ 2),
sigma_CO2_s = 0.1,
sigma_CO2_r = 0.1,
sigma_H2O_s = 0.1,
sigma_H2O_r = 0.1
)
ph = photosynthesis(lp, ep, bp, cs, use_tealeaves = FALSE, quiet = TRUE) |>
simulate_error(chamber_pars, n = 1L)
Run the code above in your browser using DataLab