05.06.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!