if (FALSE) {
data(bartlett2009.fitted)
derivatives.phenophases <- PhenoDeriv(
x=bartlett2009.fitted$fit$predicted, fit=bartlett2009.fitted$fit
)
plot(bartlett2009.fitted$fit$predicted)
abline(v=derivatives.phenophases[c(1:2,4)], col=palette())
mtext(names(derivatives.phenophases[c(1:2,4)]),
at=derivatives.phenophases[c(1:2,4)],
line=-2,
col=palette()[1:3])
}
Run the code above in your browser using DataLab