Growth model of the reared sea urchin Paracentrotus ... - SciViews

lines(edatq[, 1], predict(edat.025, newdata=list(age=edatq[, 1])), col=4)
legend(1500, 20, c("quantile 0.975", "quantile 0.5 (median)", "quantile 0.025"),
col=c(1,2,4), lty=1, pch=1)
#Rem: for the graph of Fig. 1B in the paper, it is:
#windows(8, 6)
#plot(edatq[, "age"], edatq[, "0.975"], ylim=c(0,65), xlab=expression(paste("Time ",
italic("t"), " in days")), ylab=expression(paste("Diameter ", italic("D"), " in mm")),
pch=2)
#points(edatq[, "age"], edatq[, "0.5"], pch=1)
#points(edatq[, "age"], edatq[, "0.025"], pch=6)
#lines(edatq[, 1], predict(edat.975, newdata=list(age=edatq[, 1])), lty=2)
#lines(edatq[, 1], predict(edat.5, newdata=list(age=edatq[, 1])), lty=1)
#lines(edatq[, 1], predict(edat.025, newdata=list(age=edatq[, 1])), lty=4)
#legend(1700, 18, c("quantile 0.975", "quantile 0.5 (median)", "quantile 0.025"),
lty=c(2,1,4), pch=c(2,1,6))
Pause()
cat("\n ---- Quantile regression with the fuzzy-remanent growth function ----\n\n")
cat("... it takes some time, please, be patient...\n")
# Rem: not able to calculate self-starting conditions with this data set... give
initial plausible values!
# starting values are very important here! We often get stuck in a local minimum or
the algorithm fails!
# Rem: time-scale starts at metamorphosis, and is thus shifted by 30 days (see paper)
edatf.975