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.
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