HTTEM: Hydro-thermal-time model with log-logistic distribution of...

In OnofriAndreaPG/drcSeedGerm: Utilities for data analyses in seed germination/emergence assays

Hydro-thermal-time model with log-logistic distribution of germination time (Onofri et al., 2018)

Description

This model relates the time-course of the proportion of germinated seeds to the water potential and temperature in the substrate. It is based on a truncated log-logistic distribution of germination time:

P(t) = \frac{d}{1 + exp≤ft[ b (\log(x) - \log(e)\right]}

where two of the three usual parameters ('d' and 'e') are expressed as functions of water potential (Ψ) and temperature (T). In the function 'HTTEM()', we implemented the following submodels: (1) for the parameter 'd', we implemented a shifted exponential function:

d = G \, ≤ft[ 1 - \exp ≤ft( \frac{ Ψ - Ψ_b - k(T - T_b )}{σ_{Ψ_b}} \right) \right]

while, (2) for the parameter 'e' we considered that its inverse corresponds to the median Germination Rate within the population (i.e. 1/e = GR_{50}) and modelled this latter parameter as:

GR_{50} = \frac{T - T_b }{θ_{HT}} ≤ft[Ψ - Ψ_b - k(T - T_b )\right]

The 'HTTEM.fun()' is a generic function, which can be used for plotting or other applications, while the 'HTTEM()' function is meant to be used for model fitting with the 'drmte()' function in the 'drcte()' package.

Usage

HTTEM()
HTTEM.fun(time, Psi, Temp, G, Psib, kt, Tb, sigmaPsib, ThetaHT, b)

Arguments

The 'HTTEM' function has no arguments. The 'HTTEM.fun()' has the following arguments:

time	time
Psi	water potential in the substrate
Temp	temperature
G	maximum germination capability of the seed lot
Psib	base water potential
kt	parameter measuring the effect of temperature on base water potential
Tb	base temperature
sigmaPsib	standard deviation for the base water potential within the seed lot
ThetaHT	hydro-thermal-time parameter
b	shape parameter for the cumulative distribution function of germination time

Details

The detail of this time-to-event model and the meaning of parameters are described in Onofri et al. (2018). See Table 2, where 'HTTEM()' is abbreviated as HTTE.

Value

The 'HTTEM.fun()' function returns the proportion of germinated seeds, for any given values of time and water potential in the substrate. The 'HTTEM()' function returns a list containing the nonlinear function, the parameter names and other items which are internally used by the 'drmte()' function. At the moment, there is no self-starting function and starting parameters for fitting must be provided within the 'drcte' function.

Author(s)

Andrea Onofri

References

Onofri, A., Benincasa, P., Mesgaran, M.B., Ritz, C., 2018. Hydrothermal-time-to-event models for seed germination. European Journal of Agronomy 101, 129–139. https://www.statforbiology.com/2023/stat_drcte_10-examplehtte/

Examples

data(hordeum)
modHTTEM <- drmte(nSeeds ~ timeBef + timeAf + water + temp,
                 data=hordeum,
                 fct = HTTEM(),
  start=c(0.8,-2, 0.05, 3, 0.2, 2000, 0.5))
summary(modHTTEM)

OnofriAndreaPG/drcSeedGerm documentation built on March 14, 2023, 5:45 p.m.