質問

みんな、通常は次のようなことをするとき:

tmp = zoo(rnorm(100), 1:100)
rollapply(tmp, 10, function(x) quantile(x, 0.05), align="right")

まったく正しく rollapply 瞬間から値の計算を開始します10要素が利用可能です。

残念ながら、拳10の観測にできるだけ多くのデータを使用するものが必要です。スライドウィンドウを使用するのに十分なデータがあるまで、基本的に成長しているデータのウィンドウ、例:1、1:2、1:3、1:4など。少なくとも10個の要素があるまで、いつものようにウィンドウをスライドさせるまで。

ループの醜いよりも、これを行うより良い方法はありますか?

役に立ちましたか?

解決

最初に9 NASでシリーズを埋めてみませんか? 「ループの醜い」よりも間違いなく優れています:

tmp = zoo(c(rep(NA,9), rnorm(100)), 1:109)
zoo(rollapply(tmp, 10, function(x) quantile(x, 0.05, na.rm = TRUE), 
              align="right"), 1:100)

他のヒント

rollapply 動物園では、指定することでそれを行うことができます partial=TRUE, 、例:

> library(zoo)
> 
> rollapplyr(zoo(1:20), 3, sum, partial=TRUE)
 1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 
 1  3  6  9 12 15 18 21 24 27 30 33 36 39 42 45 48 51 54 57
ライセンス: CC-BY-SA帰属
所属していません StackOverflow
scroll top