0

我有一些棘手的问题要解决，但不知何故我想不出任何答案。

我有以毫米为单位的阵雨日期和值。所以我需要做的是找出，如果连续 3 天内超过了 150 毫米的降雨量。我拥有的 data.frame ("rain") 看起来像这样。

   date amount
1. 1965-07-15     38.85315
2. 1968-08-02     33.22356
3. 1972-08-15     30.90391
4. 1978-08-07     33.34613
5. 1987-08-18     37.85918
6. 1989-07-24     38.23495
7. 1995-07-22     52.46519
8. 1995-07-23     70.29935
9. 1995-07-24     45.11350

最后，我只想得到满足条件的日期。在这种情况下，结果应该是。

7. 1995-07-22     52.46519
8. 1995-07-23     70.29935
9. 1995-07-24     45.11350

现在我试图想一些 if.condition 来解决这个问题。不幸的是，我想不出 R 中有任何可以帮助我的函数。如果有人对如何解决此案有任何想法，我将不胜感激？

提前致谢

4

1 回答 1

1

您的样本数据：

x <- read.table(header=T,
                colClasses=c(date='Date','numeric'),
text="date amount
1965-07-15     38.85315
1968-08-02     33.22356
1972-08-15     30.90391
1978-08-07     33.34613
1987-08-18     37.85918
1989-07-24     38.23495
1995-07-22     52.46519
1995-07-23     70.29935
1995-07-24     45.11350"
                )

求每个三元组的日期分布和总和：

library(zoo)
delta <- rollapply(x$date, 3, FUN=function(x) range(x)[2]-range(x)[1])
sums  <- rollapply(x$amount, 3, FUN=sum)

这些是您想要的行的（起始位置）：

delta <= 2 & sums >= 150
## [1] FALSE FALSE FALSE FALSE FALSE FALSE  TRUE

这是请求的数据子集：

w <- which(delta <= 2 & sums >= 150)
x[c(w,w+1,w+2),]
##         date   amount
## 7 1995-07-22 52.46519
## 8 1995-07-23 70.29935
## 9 1995-07-24 45.11350
于 2013-03-03T22:52:53.680 回答