R Zeitreihenobjekt ts() Minimal- und Maximaldatum

Question

Ich möchte das Datum der Minimal- und Maximalwerte in einem R ts() Objekt finden. Ich versuchte die which.min und which.max Funktionen, aber sie geben nur die "Zeilennummer" zurück. Ich möchte das aktuelle Datum ausgeben. Vielen Dank.

data <- ts(round(rnorm(60), 2), frequency = 12, end = c(2016, 12)) 
data 
which.min(data) 
which.max(data) 

Answer 1

Dies ist, wie Sie das Jahr und Monat des Min- und Max-Wert im Daten finden:

> data <- ts(round(rnorm(60), 2), frequency = 12, end = c(2016, 12)) 
> data 
     Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
2012 0.18 -0.07 -0.77 1.23 -0.97 1.20 -1.41 1.39 -0.72 -0.94 0.28 0.97 
2013 -0.86 -0.57 -0.16 -1.24 -0.35 -0.06 0.78 1.32 1.80 -0.51 -1.91 1.14 
2014 -0.51 1.21 0.14 0.30 1.18 -0.32 -0.92 -0.46 -0.97 -0.94 -1.56 -0.63 
2015 0.13 0.93 -1.45 1.97 0.04 0.55 0.45 0.13 1.14 0.27 0.15 -1.39 
2016 0.68 2.16 -1.56 -0.44 1.07 1.27 1.01 -2.93 -0.19 -0.70 1.44 0.09 

Der entsprechende Zeitpunkt des Minimums Datum ist

> data_min_value <- data[which.min(data)] 
> data_min_value 
[1] -2.93 
> data_min_value_time <- time(data)[which.min(data)] 
> data_min_value_time 
[1] 2016.583 
> data_min_value_year <- floor(time(data)[which.min(data)]) 
> data_min_value_year 
[1] 2016 
> data_min_value_month <- (time(data)[which.min(data)] %% 1)*12 
> data_min_value_month 
[1] 7 
> data_min_value_month_abb <- month.abb[(time(data)[which.min(data)] %% 1)*12+1] 
> data_min_value_month_abb 
[1] "Aug" 

Das Datum des maximalen Wert erhalten Sie in ähnlicher Weise

> data[which.max(data)] 
[1] 2.16 
> floor(time(data)[which.max(data)]) # Year 
[1] 2016 
> month.abb[(time(data)[which.max(data)] %% 1)*12+1] # Abbreviation of month 
[1] "Feb" 

Unten ist eine Zusammenfassung der hilfreichen Funktionen, die verwendet wurden in den obigen Beispielen:

> floor(2016.563) # find out the integer part on the number 
[1] 2016 
> 2016.563 %% 1 # find out the fractional part on the number 
[1] 0.563 
> month.abb[0.563*12+1] # find out the abbreviation of the month name 
[1] "Jul" 

Answer 2

Hier ist, wie ich damit umgehen, aber ich bin nicht vertraut mit ts und ich bin sicher, es ist eine bessere Option.

Um das Datum von der Max/Min-Position abzurufen, können Sie das von time erstellte Objekt unter ts indexieren. ZB: time(data)[which.max(data)]; gleich für which.min.

Dann ist diese in ein richtiges Jahr zu konvertieren (leicht) und Monat (tricky) Index, schaff' ich in der Regel diese Funktion:

set.seed(123) # for the sake of reproducibility 
data <- ts(round(rnorm(60), 2), frequency = 12, end = c(2016, 12)) 
data 
     Jan Feb Mar Apr May Jun Jul Aug 
2012 -0.56 -0.23 1.56 0.07 0.13 1.72 0.46 -1.27 
2013 0.40 0.11 -0.56 1.79 0.50 -1.97 0.70 -0.47 
2014 -0.63 -1.69 0.84 0.15 -1.14 1.25 0.43 -0.30 
2015 0.55 -0.06 -0.31 -0.38 -0.69 -0.21 -1.27 2.17 
2016 0.78 -0.08 0.25 -0.03 -0.04 1.37 -0.23 1.52 
Sep Oct Nov Dec 
2012 -0.69 -0.45 1.22 0.36 
2013 -1.07 -0.22 -1.03 -0.73 
2014 0.90 0.88 0.82 0.69 
2015 1.21 -1.12 -0.40 -0.47 
2016 -1.55 0.58 0.12 0.22 

which.min(data) 
[1] 18 
which.max(data) 
[1] 44 

numyear2monthyear(time(data)[which.max(data)]) 
[1] 2015 8 

numyear2monthyear(time(data)[which.min(data)]) 
[1] 2013 6 

Und in der Regel I:

numyear2monthyear <- function(x){ 
    c(trunc(x),     # entire part = year 
    round((x-floor(x))*12 + 1)) # decimal part * 12 + 1 (Jan=0) = Month 
} 

Hier ist ein Beispiel verwandle das in eine andere handliche Funktion, wie zum Beispiel:

Ich hoffe, es löst dein Problem (und wäre happ py, um eine bessere Option zu sehen, es zu tun).

Answer 3

Dies gibt ein Jahrmonat Klasse Objekt für das Minimum. Das gleiche mit which.max würde für das Maximum funktionieren.

library(zoo) 

imin <- which.min(data) 
tmin <- time(as.zoo(data))[imin] 
tmin 
## [1] "Dec 2016" 

as.integer(tmin), cycle(tmin) und as.Date(tmin) würde das Jahr, den Monat (als Zahl zwischen 1 und 12) geben und ein Objekt Date-Klasse jeweils den 1. des Monats verwenden.

Alternativ time(data)[imin], as.integer(time(data))[imin] und cycle(data)[imin] wird die Zeit als Jahr + Fraktion, Jahr als integer und Monatszahl geben. Diese drei verwenden keine Pakete.

Answer 4

Mein Lieblings-Tool für Zeitreihen ist xts und ts Objekte übersetzen sauber:

library(xts) 
x = as.xts(data) 

Ausgänge:

> min(index(x)) 
[1] "Jan 2012" 
> max(index(x)) 
[1] "Dec 2016"