Subsetting Winter (Dez, Jan, Feb) de la serie temporal diaria (zoológico)

Question

Tengo un zoológico diario (XTS) con algunas décadas de datos en el siguiente formato:

head(almorol)
1973-10-02 1973-10-03 1973-10-04 1973-10-05 1973-10-06 1973-10-07
     183.9      208.2      153.7       84.8       52.5       35.5

Y me gustaría trazar solo los datos de invierno (los meses completos de diciembre, enero y febrero). Encontré el subconjunto para XTS, así que pensé que podría extraer todos los decembers usando:

x<-apply.yearly(almorol, FUN=last(almorol, "1 month"))

Y luego haz algo similar para Jan y Feb, pero recibo el siguiente error:

Error in get(as.character(FUN), mode = "function", envir = envir) :
object 'FUN' of mode 'function' was not found

Puedo usar el apply.yearly y last(almorol, "1 month") por separado, pero cuando los combino no funciona. ¿Alguien conoce una forma de submirmar esos 3 meses de la serie temporal? ¡Gracias por ayudar!

Answer 1

Prueba esto:

z.winter <- z[months(time(z), TRUE) %in% c("Dec", "Jan", "Feb")]
plot(z.winter)