2

単純な動物園のオブジェクトは

z <- zoo(11:15, as.Date(31:45))
as.data.frame(z) 
            z
1970-02-01 11
1970-02-02 12
1970-02-03 13
1970-02-04 14
1970-02-05 15
1970-02-06 11
1970-02-07 12
1970-02-08 13
1970-02-09 14
1970-02-10 15
1970-02-11 11
1970-02-12 12
1970-02-13 13
1970-02-14 14
1970-02-15 15
as.data.frame(rollapply(z, width=3, mean, align="right"))
           rollapply(z, width = 3, mean, align = "right")
1970-02-03                                       12.00000
1970-02-04                                       13.00000
1970-02-05                                       14.00000
1970-02-06                                       13.33333
1970-02-07                                       12.66667
1970-02-08                                       12.00000
1970-02-09                                       13.00000
1970-02-10                                       14.00000
1970-02-11                                       13.33333
1970-02-12                                       12.66667
1970-02-13                                       12.00000
1970-02-14                                       13.00000
1970-02-15                                       14.00000

開始観測を固定することは可能ですか?つまり、日付「1970-02-04」では、平均が「12, 13, 14」ではなく「11, 12, 13, 14」などになりますか? もちろん、ループを介してそれを行うこともできますが、問題がもう少し複雑な場合、ループは非常に遅く、apply、またはと同じくらい速いものを探していましたrollapply。

4

1 に答える 1

3

シリーズの長さに等しい幅を設定し、次を使用しますpartial = TRUE。

rollapplyr(z, length(z), mean, partial = TRUE)

使用することもできますcumsum(z) / seq_along(z)

于 2015-03-23T15:33:49.583 に答える