I am working with 16-day MODIS EVI (satellite) data and I want to fit a Piecewise Logistic Function through my 23 EVI data values. The following formula is for the Piecewise Logistic Function:
y(t)= [c / (1 + e^(a+bt))] + EVIb
t is the day of year
a is vegetation growth time
b is rate of leaf development
c is the amplitude
bEVI is the background EVI value [constant]
For some reason there is very little clarification about each component of the equation in most papers. I am assuming that the bEVI value is a constant but am supposing that one or two other components are constants too?
I have little experience with this function, but from the literature I have read, is applied separately to the left and then subsequently to the right of this bell-shaped curve - split at the maximum y value. I am hoping someone here with a stronger mathematical background and experience using the piecewise logistic function would be able to advise me on which component of the equation above should be a constant and which should be a variable.
I tried embedding the original journal paper that I have on my OneDrive but can't yet seem to figure out how to do it. I have, however, included the link here "Monitoring Vegetation Phenology using MODIS (Zhang et al., 2003)". You may require access to Science Direct but please do let me know if you have any problems with that. The piecewise logistic regression is mentioned on page 472, section 2.
I contacted the original author regarding the equation and he clarified that the function if applied separately for the left and right side of the curve, split at the maximum point, which could most likely be somewhere in the middle of the curve.
The series of days of the year are as follows:
Here is an example EVI dataset I currently have for one year of 23 data points: