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.

are inside this interval. Additionally, a τ = 0.5 median curve visualizes<br />

possible asymmetry and its change with time. When a representative<br />

sample <strong>of</strong> <strong>the</strong> whole size distribution at each time is available, typically<br />

with n > 50, points representing unconditional quantiles can be superposed<br />

on <strong>the</strong> graph to show how well quantile regressions match <strong>the</strong>m, as we did<br />

in Figs. 28B, 31 and 34A. There is no residuals analysis for this kind <strong>of</strong><br />

quantile regression, at least in <strong>the</strong> way it is conceived for least-square<br />

regression.<br />

Nonlinear quantile regression using interior point algorithm proposed<br />

by Koenker & Park (1996) is compatible with many growth <strong>model</strong>s. For<br />

<strong>the</strong> cases where representative samples are measured at each time interval,<br />

we have proposed a modified method to use with so-called 'envelope<br />

<strong>model</strong>s', that is, <strong>model</strong>s <strong>of</strong> <strong>the</strong> form Y = f(t, τ) as eq. 35. These <strong>model</strong>s<br />

calculate a 3D-surface enveloping data (Fig. 33). Curves for all quantiles<br />

0 < τ < 1 are calculated at once. They contain most <strong>of</strong> <strong>the</strong> information in<br />

<strong>the</strong> initial dataset, including individual variability. They are particularly<br />

useful when individual variability in itself expresses some aspects <strong>of</strong> <strong>the</strong><br />

phenomenon studied, like growth in <strong>the</strong> presence <strong>of</strong> an intraspecific<br />

competition. We have developed basic tools to fit and diagnose <strong>the</strong>m<br />

(namely, basic graphical residuals analysis, Fig. 34). Tests <strong>of</strong> significance<br />

<strong>of</strong> <strong>the</strong> fitting can be elaborated as extensions <strong>of</strong> some tests for<br />

unconditional size-distributions, like χ 2 or Kolmogorov-Smirnov tests or<br />

by o<strong>the</strong>r means (Koenker & Machado, 1999). Confidence intervals for <strong>the</strong><br />

parameters do not yet exist. Tests for comparing various <strong>model</strong>s fitted on<br />

<strong>the</strong> same dataset, or to compare a single <strong>model</strong> fitted on various datasets<br />

remain also to be developed. However, one difficulty in conceptualizing<br />

such tests for envelope <strong>model</strong>s is that one <strong>of</strong> <strong>the</strong> "independent variables",<br />

τ, is estimated according to <strong>the</strong> dependent variable y (eq. 36) and is thus<br />

not really independent.<br />

Flexible, unconstrained envelope <strong>model</strong>s can be incredibly complex,<br />

with dozens <strong>of</strong> parameters. They are impossible to fit in practice.<br />

Part IV: A growth <strong>model</strong> with intraspecific competition<br />

165

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

Saved successfully!

Ooh no, something went wrong!