在R中的循环系列中找到3个连续的最高值
Find 3 consecutive highest values in circular series in R
我有以下形式的数据（2 个示例）：
p1 <- structure(c(1.38172177074188, 1.18601365390563, 1.25131938561825,
1.07175353794277, 0.887770295772917, 0.806599968169486, 0.843543355495394,
0.889051695167723, 0.764131945540256, 0.699309441111923, 0.945165791967098,
1.31310409471336), .Dim = 12L)
p2 <- structure(c(1.24801075135611, 1.06280347993594, 1.21410288703334,
1.36797720634294, 1.07291218307332, 0.936924063490867, 0.819723966406961,
0.854960740335283, 0.718565087630857, 0.649827141012991, 0.785853771875901,
1.04368795443605), .Dim = 12L)
这些是水文时间序列的标准化月平均值；所谓的 Pardé 制度，它给出了一些关于年度季节性的迹象。要进行进一步分析，我需要从这些 Pardé 系列中得出最高和最低的 3 个月份。因为季节性可能是双峰的，所以我需要确定 3 个 highest/lowest 连续 数据点（通常不是三个绝对 highest/lowest 数据点，请参见示例）得出最潮湿和最干燥时期的时间。到目前为止，由于时间序列的循环特性，我失败了，这提出了一个特殊的挑战。
有什么建议吗？
您可以使用 filter。它对连续值求和，可以处理循环时间序列。
f1 <- stats::filter(p1, c(1, 1, 1), circular = TRUE, sides = 1)
#Time Series:
# Start = 1
#End = 12
#Frequency = 1
#[1] 3.639992 3.880840 3.819055 3.509087 3.210843 2.766124 2.537914 2.539195 2.496727 2.352493 2.408607 2.957579
((which.max(f1) - (3:1)) %% 12) + 1
#[1] 12 1 2
我有以下形式的数据（2 个示例）：
p1 <- structure(c(1.38172177074188, 1.18601365390563, 1.25131938561825,
1.07175353794277, 0.887770295772917, 0.806599968169486, 0.843543355495394,
0.889051695167723, 0.764131945540256, 0.699309441111923, 0.945165791967098,
1.31310409471336), .Dim = 12L)
p2 <- structure(c(1.24801075135611, 1.06280347993594, 1.21410288703334,
1.36797720634294, 1.07291218307332, 0.936924063490867, 0.819723966406961,
0.854960740335283, 0.718565087630857, 0.649827141012991, 0.785853771875901,
1.04368795443605), .Dim = 12L)
这些是水文时间序列的标准化月平均值；所谓的 Pardé 制度，它给出了一些关于年度季节性的迹象。要进行进一步分析，我需要从这些 Pardé 系列中得出最高和最低的 3 个月份。因为季节性可能是双峰的，所以我需要确定 3 个 highest/lowest 连续 数据点（通常不是三个绝对 highest/lowest 数据点，请参见示例）得出最潮湿和最干燥时期的时间。到目前为止，由于时间序列的循环特性，我失败了，这提出了一个特殊的挑战。 有什么建议吗？
您可以使用 filter。它对连续值求和，可以处理循环时间序列。
f1 <- stats::filter(p1, c(1, 1, 1), circular = TRUE, sides = 1)
#Time Series:
# Start = 1
#End = 12
#Frequency = 1
#[1] 3.639992 3.880840 3.819055 3.509087 3.210843 2.766124 2.537914 2.539195 2.496727 2.352493 2.408607 2.957579
((which.max(f1) - (3:1)) %% 12) + 1
#[1] 12 1 2