0

解決しなければならないトリッキーな質問がありますが、どういうわけか答えがわかりません。

にわか雨の日付と値は mm です。そこで私がしなければならないことは、3 日連続で 150mm の雨量を超えたかどうかを調べることです。私が持っているdata.frame（「雨」）は次のようになります。

   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 に答える