R-код для расчета рекламного запаса / распада с пределом эффекта переноса
Учитывая вектор х сказать,
x <- c(10,0,0,20,0,0)
и рекламный запас 0,8 - как рассчитать эффект рекламного запаса так, чтобы, если эффект переноса был меньше некоторого предела k, скажем, половины первоначального значения, он больше не переносился.
Без этого условия я могу просто использовать это:
filter(x, filter=0.8, method="recursive")
Который дает
10.0000 8.0000 6.4000 25.1200 20.0960 16.0768
Но учитывая это условие, если мы применяем запас от 0,2 до 10, который дает,
10.0000 8.0000 6.4000 5.1200 4.096 3.2768 ...
эффект переноса должен остановиться на 5,12, поскольку 4.096 меньше половины 10. Таким образом, ожидаемый результат будет:
10.0000 8.0000 6.4000 25.1200 16.0000 12.8000