4

我正在尝试在 R 中复制以下公式:

Xt = Xt-1 * b + Zt * (1-b)

我正在使用以下代码

t %>% 
  mutate(x= ifelse(week == 1, z, NaN)) %>% # Initial value for the first lag
  mutate(x= ifelse(week != 1, lag(x,1 ,default = 0) * b + z, z)

但是除了第二个元素,我得到了所有的 NaN。

     z      b           x
     (dbl)  (dbl)    (dbl)
1  168.895  0.9      168.8950
2   20.304  0.9      131.7472
3   14.943  0.9         NA
4   11.028  0.9         NA
5    8.295  0.9         NA
6    8.024  0.9         NA
7    6.872  0.9         NA
8    7.035  0.9         NA
9    4.399  0.9         NA
10   4.158  0.9         NA

这在 excel 中相当简单,但我必须在 R 中完成,你有什么方法吗?

可重现的例子:

set.seed(2)
t = data.frame(week = seq(1:52),
               z = runif(52, 0, 100),
               b = 0.2)
4

1 回答 1

1

感谢@Frank 和@docendo discimus,我找到了运行以下循环的解决方案

for (row in 2:dim(t)[1]) {
    t[row,] <- mutate(t[1:row,], x= lag(x,1) * b + z * (1 - b))[row,]
  }
于 2016-06-29T09:18:40.227 回答