假設我有以下系列資料框。
library(zoo)
n=10
date = as.Date(1:n);date
y = rnorm(10);y
dat = data.frame(date,y)
dat
date y
1 1970-01-02 -0.02052313
2 1970-01-03 0.28255304
3 1970-01-04 -0.10718621
4 1970-01-05 -1.19299366
5 1970-01-06 1.17072468
6 1970-01-07 0.55849119
7 1970-01-08 0.30474050
8 1970-01-09 -0.30777180
9 1970-01-10 -0.01874367
10 1970-01-11 -0.74233556
計算寬度（滾動視窗）為 3 天的滾動標準偏差，我這樣做：
a = zoo::rollapply(dat$y,3,sd)
結果：
2 3 4 5 6
1.7924102 0.4189522 0.4599979 0.4164408 0.3786601
7 8 9
0.9481849 0.7048868 0.2494578
并找到最大標準偏差
max(a)
1.79241
現在我想知道這個最大值指的是哪個時間 3 天的時間間隔。我該怎么做？假設我的系列是 20 年，所以我想找到這 20 年的最大標準差，并提取具體的 3 天時間間隔。
uj5u.com熱心網友回復：
``` r
set.seed(123)
n=10
date = as.Date(1:n, origin = "1970-01-01")
y = rnorm(10)
dat = data.frame(date ,y)
roll_win <- 3
dat$a = c(rep(NA_real_, roll_win - 1), zoo::rollapply(dat$y, roll_win ,sd))
dat <- subset(dat, !is.na(a))
dat_max <- dat[dat$a == max(dat$a, na.rm = TRUE), ]
dat_max$date_start <- dat_max$date - (roll_win - 1)
dat_max
#> date y a date_start
#> 8 1970-01-09 -1.265061 1.496275 1970-01-07
由reprex 包于 2022-02-18 創建(v2.0.1)
uj5u.com熱心網友回復：
在居中視窗的情況下，您可以執行以下操作：
step = (3-1)/2
wmax = which.max(a)
int = (wmax - step):(wmax step)
dat[int[int > 0],] #Avoids error when wmax < step.
# date y
# 1 1970-01-02 -0.6264538
# 2 1970-01-03 0.1836433
# 3 1970-01-04 -0.8356286
#Data
library(zoo)
set.seed(1)
n=10
date = as.Date(1:n)
y = rnorm(10)
dat = data.frame(date,y)
dat
a = rollapply(dat$y,3,sd)
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/429344.html