I am working with climate data which I want to split up into seasons (DJF, MAM, JJA and SON) and apply a function to each individual season. xts offers apply.quarterly but that splits it up into quarters (JFM, AMJ, JAS, OND).
How can I do this seasonally with period.apply? endpoints also only offers quarters and not seasons.
I am using xts and zoo.