get_pars_s
Obtains parameters from growth and reprod data at multiple food levels
Contents
Syntax
[par, U, Ep] = get_pars_s(p, w, par0)
Description
Obtains parameters from growth and reprod data at multiple food levels Maturity and somatic maintenance rate coefficients might differ. The parameter kapR is fixed
Input
- p: (n,6)-matrix with quantities (n > 1). The columns are:
1 f % - % scaled functioal response
2 L_b % mm % length at birth
3 L_p % mm % length at puberty
4 L_i % mm % ultimate length
5 \dot{r}_B % d^-1 % von Bertalanffy growth rate
6 \dot{R}_i % % d^-1 % maximum reproduction rate- w: (n,5)-matrix with weight coefficient (optional, default is ones)
- par0: optional (8,2)-matrix with initial estimate for par
Output
- par: (8,2)-matrix with first column
1 kap % - % fraction allocated to som maint + growth
2 kapR % - % fraction of energy allocated to reprod fixed in embryo
this parameter is set in fnget_pars_s
3 g % - % energy investment ratio
4 kJ % d^-1 % maturity maintenance rate coefficient
5 kM % d^-1 % somatic maintenance rate coefficient
6 v % mm/d % energy conductance
7 Hb % d mm^2 % scaled maturity at birth M_H^b/{J_EAm}
8 Hp % d mm^2 % scaled maturity at puberty M_H^p/{J_EAm}
choice of sequence for consistency with get_par_r
second column: zeros (fix) or ones (iterate)- U: (n,3)-matrix with U^0 = M_E^0/\{\dot{J}_{EAm}\}, U^b = M_E^b/\{\dot{J}_{EAm}\}, U^p = M_E^p/\{\dot{J}_{EAm}\}
- Ep: (n,6)-matrix as p, but now based on DEB parameters
Remarks
iget_pars_s is inverse to get_pars_s
Example of use
See mydata_get_pars