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