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.

plot(edat0q[, "age"], edat0q[, "0.975"], ylim=c(0,65), xlab="Time elapsed from<br />

metamorphosis in days", ylab="Diameter increase in mm", main="Quantile regressions<br />

with SSfuzremOrig1")<br />

points(edat0q[, "age"], edat0q[, "0.5"], col=2)<br />

points(edat0q[, "age"], edat0q[, "0.025"], col=4)<br />

lines(edat0q[, 1], predict(edat0f.975, newdata=list(age=edat0q[, 1])), col=1)<br />

lines(edat0q[, 1], predict(edat0f.5, newdata=list(age=edat0q[, 1])), col=2)<br />

lines(edat0q[, 1], predict(edat0f.025, newdata=list(age=edat0q[, 1])), col=4)<br />

legend(1500, 20, c("quantile 0.975", "quantile 0.5 (median)", "quantile 0.025"),<br />

col=c(1,2,4), lty=1, pch=1)<br />

#Rem: for <strong>the</strong> graph <strong>of</strong> Fig.4 in <strong>the</strong> paper, it is:<br />

#windows(8, 6)<br />

#plot(edat0q[, "age"], edat0q[, "0.975"], ylim=c(0,65), xlab=expression(paste("Time ",<br />

italic("t'"), " in days")), ylab=expression(paste("Diameter increase ", italic("D'"),<br />

" in mm")), pch=2)<br />

#points(edat0q[, "age"], edat0q[, "0.5"], pch=1)<br />

#points(edat0q[, "age"], edat0q[, "0.025"], pch=6)<br />

#lines(edat0q[, 1], predict(edat0f.975, newdata=list(age=edat0q[, 1])), lty=2)<br />

#lines(edat0q[, 1], predict(edat0f.5, newdata=list(age=edat0q[, 1])), lty=1)<br />

#lines(edat0q[, 1], predict(edat0f.025, newdata=list(age=edat0q[, 1])), lty=4)<br />

#legend(1700, 18, c("quantile 0.975", "quantile 0.5 (median)", "quantile 0.025"),<br />

lty=c(2,1,4), pch=c(2,1,6))<br />

Pause()<br />

#Rem: <strong>the</strong> following code takes very long to run, so it is deactivated<br />

#just remove comment marks to run it...<br />

#cat("\n ---- Calculating parameters estimation for quantiles every 5% step ----\n\n")<br />

#cat("... it takes some time, please, be patient...\n")<br />

#edat0f.pars

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

Saved successfully!

Ooh no, something went wrong!