getSens: Estimate temperature sensitivty value and time varying $R_b$...
In dpabon/ecofunr: Deriving Ecosystem Functional Properties

Function to determine the temperature sensitivity sensitivity and time varying

General Function to determine the temperature sensitivity (S value) and time varying basal efflux (R_b) from a given temperature and efflux (usually respiration) time series. The following general model is assumed:

Resp(i) = Rb exp(S/tau),

where i is the time index. It has been shown, however, that this model is misleading when Rb is varying over time which can be expected in many real world examples (e.g. Sampson et al. 2008).

If Rb varies slowly, i.e. with some low frequency then the 'scale dependent parameter estimation, SCAPE' allows us to identify this oscillatory pattern. As a consequence, the estimation of Q10 can be substantially stabilized (Mahecha et al. 2010). The model becomes

Resp(i) = Rb(i) * Q10^((T(i)-Tref)/(gamma),

where Rb(i) is the time varying 'basal respiration', i.e. the respiration expected at Tref. The convenience function getQ10 allows to extract the Q10 value minimizing the confounding factor of the time varying Rb. Four different spectral methods can be used and compared. A surrogate technique (function by curtsey of Dr. Henning Rust, written in the context of Venema et al. 2006) is applied to propagate the uncertainty due to the decomposition.

The user is strongly encouraged to use the function with caution, i.e. see critique by Graf et al. (2011).

getSens(
  temperature,
  respiration,
  sf,
  gettau,
  fborder = 30,
  M = -1,
  nss = 0,
  method = "Fourier",
  weights = NULL,
  lag = NULL,
  gapFilling = TRUE,
  doPlot = FALSE
)

`temperature`	numeric vector: temperature time series
`respiration`	numeric vector: respiration time series
`sf`	numeric: sampling rate, number of measurements (per day)
`gettau`	numeric: function to transform the exponent in the sensitivity model
`fborder`	numeric: boundary for dividing high- and low-frequency parts (in days)
`M`	numeric vector: size of SSA window (in days)
`nss`	numeric vector: number of surrogate samples
`method`	String: method to be applied for signal decomposition (choose from 'Fourier','SSA','MA','EMD','Spline')
`weights`	numeric vector: optional vector of weights to be used for linear regression, points can be set to 0 for bad data points
`lag`	numeric vector: optional vector of time lags between respiration and temprature signal
`gapFilling`	Logical: Choose whether Gap-Filling should be applied
`doPlot`	Logical: Choose whether Surrogates should be plotted