Displaying 1 result from an estimated 1 matches for "dateq50".
2012 Feb 14
1
cumsum function to determine plankton phenology
...m<-cumsum(datayear$sp)/max(cumsum(datayear$sp))
# take index of minimum value of cumsum
idq15<-which.min(mean(.15-datayear$cumsum))
idq50<-which.min(mean(.50-datayear$cumsum))
idq85<-which.min(mean(.85-datayear$cumsum))
#take corresponding date
dateq15<-datayear$date[idq15]
dateq50<-datayear$date[idq50]
dateq85<-datayear$date[idq85]
# output results in day in year form
print(format(c(dateq15,dateq50,dateq85),"%j")
# results should be similar to this:
[1] "2010-03-02 CET" "2010-06-29 CEST" "2010-10-27 CEST"
with dates formatted...