f2o
From scaled functional response trajectories to otolith opacity profiles
Contents
Syntax
[LOO, teLLOO] = f2o (tcf, eb, LOb, par)
Description
From scaled functional response f to otolith opacity O profiles. The routine o2f reconstructs time-trajectories of scaled function response, scaled reserve density, body length, otolith length, opacity from otolith length-opacity data, given DEB parameters and temperature correction factors over time. The reconstructed body length at otolith collection can be compared to the observed body length as an independent check on the reliability of the reconstruction. The assumption kJ = kM is made to ensure stage transitions at fixed length. Opacity is zero for contribution of dissipation and one of growth. Opacity, therefore, is always on the interval (0,1). Measured opacity might need some rescaling to fit these requirements. Otolith and body length are treated as volumetric length; measured lengths might need rescaling. The theory behind o2f and f2o is discussed in PecqKooy2008
Input
- tcf: (n,3)-matrix with time t, tmep cor factor c_T, func resp f
- eb: scalar with e at t=0
- LOb: scalar with otolith length LO at t=0 (birth)
- par: 10-vector with parameters Lb, Lp, v, vOD, vOG, kM, g, kap, kapR, delS
Output
- LOO: (n,2)-matrix with L_O, O
- teLLOO: (n,5)-matrix with time, e, L, LO, O
Remarks
Routine f2o is reverse to o2f and can be used for checking.
Example of use
See mydata_of
code
unpack par