Growth model of the reared sea urchin Paracentrotus ... - SciViews
Growth model of the reared sea urchin Paracentrotus ... - SciViews
Growth model of the reared sea urchin Paracentrotus ... - SciViews
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
d. Results<br />
Theoretical considerations: use <strong>of</strong> quantile regression instead<br />
<strong>of</strong> least-square regression<br />
The most widespread method to fit a curve is <strong>the</strong> use <strong>of</strong> a least-square<br />
method with one <strong>of</strong> <strong>the</strong> many minimization algorithms available [simplex,<br />
(quasi-)Newton, etc.; Sen & Srivastava, 1990; Draper & Smith, 1998;<br />
Nocedal & Wright, 1999]. The algorithm finds <strong>the</strong> combination <strong>of</strong> values<br />
for <strong>the</strong> various parameters in <strong>the</strong> <strong>model</strong> (<strong>the</strong> solution) that leads to a<br />
minimal value for <strong>the</strong> objective function, which is here <strong>the</strong> sum <strong>of</strong> <strong>the</strong><br />
square <strong>of</strong> <strong>the</strong> residuals (that is, <strong>the</strong> sum <strong>of</strong> squared distances between<br />
observed values for <strong>the</strong> dependent variable and values predicted by <strong>the</strong><br />
<strong>model</strong> at <strong>the</strong> same levels for <strong>the</strong> independent variables).<br />
Least-square regression has many advantages over o<strong>the</strong>r methods. In<br />
particular, when partial first (gradient matrix) and second (Hessian matrix)<br />
derivatives <strong>of</strong> <strong>the</strong> function are calculable for each parameter, convergence<br />
through a solution is accelerated and can be verified (at least for a local<br />
solution, Nocedal & Wright, 1999). In <strong>the</strong> counterpart, that regression<br />
supposes that <strong>the</strong> fluctuation around <strong>the</strong> <strong>model</strong> (called <strong>the</strong> error term) is<br />
additive, independent, normally distributed and with a constant standard<br />
deviation (heteroscedasticity). It is also very sensitive to outliers because it<br />
uses <strong>the</strong> squared residuals. Those constraints, even if not strictly met every<br />
time, particularly in many nonlinear phenomena like growth, appear to be<br />
<strong>of</strong> minor importance for many authors. Indeed, outliers are eliminated, or<br />
weighing methods are applied to limit <strong>the</strong>ir impact.<br />
Yet, <strong>the</strong>re are two arguments against <strong>the</strong> least-square regression used in<br />
<strong>the</strong> framework <strong>of</strong> growth <strong>model</strong>s, particularly when individuals' growth is<br />
influenced by <strong>the</strong> presence <strong>of</strong> conspecifics or <strong>of</strong> o<strong>the</strong>r species (indeed a<br />
general case to be verified in each study, except when a single individual is<br />
grown alone in a cage or an aquarium!): first it is sensitive to <strong>the</strong><br />
Part IV: A growth <strong>model</strong> with intraspecific competition<br />
144