Умножая все элементы вектора в R

Я хочу, чтобы функция возвращала произведение всех значений в векторе, например sum, но с умножением вместо добавления. Я ожидал, что это уже существует, но если это произойдет, я не смогу его найти. Здесь мое решение:

product <- function(vec){
    out <- 1
    for(i in 1:length(vec)){
         out <- out*vec[i]
    }
    out
}

Это ведет себя так, как я хочу. Например:

> product(1:3)
[1] 6

Есть ли лучший способ сделать это, либо с существующей функцией, либо с помощью улучшения этого пользовательского?