R在我的研究中計算XTS（或data.frame）對象中的相關日期

Question

我正在分析氣象站的歷史風速（如果趨勢正在增加或減少等）。 現在我每天都有一個價值。 我想統計彼此相關的日子（例如1958-03-18,19,20 ......）。對於每年我想有一個價值連續3天以上發生多少（風暴）天。

例如：

1982-01-30 41.04 
1982-02-02 45.72 
1982-02-03 46.8 
1982-02-04 41.04 
1982-02-12 39.24 
1982-02-17 53.28 
1982-02-18 49.68 
1982-02-19 40.32 
1982-03-01 46.08 

在1982年2月的2倍的情況是（2,3,4）和（17,18,19）。

有沒有人知道如何計算這個，並把它放在一個新的表中進一步分析/繪圖？

1982 23 
1983 7 
1984 11 
. 
. 
. 

至少每年計算所有日子也會對我有幫助。

親切的問候 薩沙

Answer 1

下面是（根據您的評論風> 39）一個簡單的方法來識別暴雨天，有三個連續的暴雨天的日期。

#create example data 
fakedates <- as.Date("1982-01-01")+(1:730) 
fakewind <- sample(1:80, length(fakedates), replace = TRUE) 
dateyear <- format(fakedates, "%Y") 

df <- data.frame("dates"=fakedates, "maxwind"=fakewind, "year"=dateyear) 

#identify 'storm' days where max wind is >39 
df$storm <- df$maxwind > 39 #identify storm days 

# ensure original data is in chronological order, otherwise you would need to sort first 
# e.g. if(is.unsorted(original$dates), sort(original$dates, .... etc) 

#identify 3-day storms - undetermined for last two days in the record 
df$storm3 <- 
    c(df$storm[1:(length(df$storm) - 2)] & 
     df$storm[2:(length(df$storm) - 1)] & 
     df$storm[3:(length(df$storm))], "unknown", "unknown") 

#select all storm day records 
dfs <- df[df$storm == TRUE, c(1:3)] 

#select only dates of 3-day storms 
dfs3 <- df[df$storm3== TRUE,(1:3)] 

你的輸出就會是這個樣子：

> head(dfs3) 
     dates maxwind year 
12 1982-01-13  57 1982 
13 1982-01-14  47 1982 
33 1982-02-03  65 1982 
49 1982-02-19  79 1982 
61 1982-03-03  49 1982 
68 1982-03-10  55 1982 

您可以通過一年總結，做你需要的任何其他分析，無論是與所有困天（DFS）或三天的暴雨天只有（dfs3）。

當您更熟悉R時，您可以通過各種其他方式解決任務，例如與包裹dplyr