R cumsum with value-based multiplication

Question

I have a vector of positive numbers such that - if the number is greater than or equal to 1, it is added - if the number is less than 1, it multiplies the cumulative sum.

for instance

> set.seed(0)
> x <- sample(c(0.9, 0.9, 0.9, 1:3), 10, replace=T)
> x
 [1] 3.0 0.9 0.9 1.0 3.0 0.9 3.0 3.0 1.0 1.0

Received Vector

3  2.7  2.43  3.43  6.43  5.787  8.787  11.787  12.787  13.787

Is there a way to do this without using a for loop?

Answer 1

Here are a few alternatives using only the R base:

1) Reduce

sumprod <- function(x, y) if (y < 1) x * y else x + y
Reduce(sumprod, x, acc = TRUE)
## [1]  3.000  2.700  2.430  3.430  6.430  5.787  8.787 11.787 12.787 13.787

2) Recursion

cumsumprod_ <- function(x, init) {
      if (length(x)) c(init, Recall(x[-1], sumprod(tail(init, 1), x[1])))
      else init
}
cumsumprod <- function(x) {
      if (length(x) < 2) x
      else cumsumprod_(x[-1], x[1])
}

cumsumprod(x)
## [1]  3.000  2.700  2.430  3.430  6.430  5.787  8.787 11.787 12.787 13.787

Update: Fixed edge in (2).

Answer 2

And here the tidyversealternative @ G.Grothendieck answers with purrr::accumulate:

library(tidyverse)
accumulate(x, function(x, y) if_else(y < 1, x * y, x + y))
#>  [1]  3.000  2.700  2.430  3.430  6.430  5.787  8.787 11.787 12.787 13.787