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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Finally, ∆D∞ should follow a normal distribution, as size distributions<br />

when approaching asymptotic maximum size are normal or close to<br />

normal (see Fig. 28A, t > 1500 days, and also Grosjean et al, 1996, see<br />

Part III):<br />

with µ D∞<br />

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

∆D ( τ) ∼ N ( µ , σ )<br />

(34)<br />

∞ ∆D ∆D<br />

∆ being <strong>the</strong> mean and σ ∆D∞<br />

normal distribution <strong>of</strong> ∆D∞(τ).<br />

∞ ∞<br />

Replacing eqs. 32-34 into eq. 31, we get:<br />

−k1⋅t' 1+ e<br />

D'( t', τ) =∆D<br />

( τ)<br />

1 −s⋅(1 −τ) ⋅e<br />

being <strong>the</strong> standard deviation <strong>of</strong> <strong>the</strong><br />

∞ −k2⋅t' (35)<br />

which links curves for all quantiles 0 < τ < 1 and has 5 parameters to be<br />

estimated: k1, k2, s, µ D∞<br />

∆ and σ ∆ D∞<br />

. It includes individual variations into<br />

<strong>the</strong> <strong>model</strong> and is a kind <strong>of</strong> 3D-surface that envelops data (see Fig. 33). For<br />

this reason, it will be called an 'envelope <strong>model</strong>'.<br />

The quantile regression method is modified as follows. Considering<br />

that every individual present in <strong>the</strong> batch is measured at each sampling<br />

time, unconditional quantiles at each size distribution can be regarded as<br />

estimators <strong>of</strong> conditional quantiles τ at corresponding time t' in eq. 35<br />

(note that this is fundamentally different than <strong>the</strong> previous quantile<br />

regression method in eqs. 21-22 where unconditional quantiles were not<br />

used at all in <strong>the</strong> regression). Estimators <strong>of</strong> τ, noted ˆ τ , are <strong>the</strong>n calculated<br />

as:<br />

ˆ τ =<br />

i<br />

nt' ( i )<br />

∑<br />

j=<br />

1<br />

( D'j t'i < D'i)<br />

I ( )<br />

nt' ( )<br />

i<br />

(36)<br />

where t'i is t' corresponding to <strong>the</strong> i th observation, n(t'i) is <strong>the</strong> total number<br />

<strong>of</strong> individuals measured at time t'i, D'j(t'i) is <strong>the</strong> j th observation among all<br />

160

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

Saved successfully!

Ooh no, something went wrong!