我想在r ts()
对象中找到最小值和最大值的日期。我尝试了which.min
和which.max
功能,但它们仅返回"行号"。我想输出实际日期。谢谢。
data <- ts(round(rnorm(60), 2), frequency = 12, end = c(2016, 12))
data
which.min(data)
which.max(data)
我最喜欢的时间序列工具是 xts
, ts
对象干净地翻译:
library(xts)
x = as.xts(data)
输出:
> min(index(x))
[1] "Jan 2012"
> max(index(x))
[1] "Dec 2016"
这是您在数据中找到最小和最大值的年份的方式:
> 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
最低日期的相应日期为
> 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"
您获得的最大值的日期类似
> 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"
以下是关于上述示例中使用的有用功能的摘要:
> 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"
这是我处理的方式,但我对ts
不熟悉,我敢肯定有一个更好的选择。
要从最大/min位置检索日期,您可以在ts
上索引time
创建的对象。例如:time(data)[which.max(data)]
;which.min
。
然后将其转换为适当的年度(简单)和月份(棘手)索引,我通常会创建此功能:
numyear2monthyear <- function(x){
c(trunc(x), # entire part = year
round((x-floor(x))*12 + 1)) # decimal part * 12 + 1 (Jan=0) = Month
}
这是一个示例:
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
通常我将其变成另一个方便的功能,例如:
extrema_dates <- function(ts){
ts_min_date <- numyear2monthyear(time(ts)[which.min(ts)])
ts_max_date <- numyear2monthyear(time(ts)[which.max(ts)])
list(min=min(ts),
min_year=ts_min_date[1],
min_month=ts_min_date[2],
max=max(ts),
max_year=ts_max_date[1],
max_month=ts_max_date[2])
}
> extrema_dates(data)
$min
[1] -1.97
$min_year
[1] 2013
$min_month
[1] 6
$max
[1] 2.17
$max_year
[1] 2015
$max_month
[1] 8
我希望它能解决您的问题(并且很高兴看到这样做的更好选择)。
这给出了最小的年度类对象。与which.max
相同的最大值。
library(zoo)
imin <- which.min(data)
tmin <- time(as.zoo(data))[imin]
tmin
## [1] "Dec 2016"
as.integer(tmin)
,cycle(tmin)
和as.Date(tmin)
分别使用本月的第一个给出年度,月份(以1到12之间的数字)和日期类对象。
交替,time(data)[imin]
,as.integer(time(data))[imin]
和 cycle(data)[imin]
将把时间作为年份 分数的时间一年为整数和月号。这三个不使用任何软件包。
如果您有类似日期列和另一个值(例如降雨)列的东西,则可以随时以这种方式返回值。如果使用动物园,您可能需要转换为数据框。
# Get your max and min Rainfall values. The brackets return the value as well as set it in a variable
(maxR <- max(df$Rainfall))
(minR <- min(df$Rainfall))
# Get the Dates of your max and min Rainfall values.
(maxD <- df$Date[df$Rainfall==maxR])
(minD <- df$Date[df$Rainfall==minR])