[Rcpp-devel] R vectorisation vs. C++ vectorisation

[Hadley Wickham]

> To be fair with the R language, I would have compare with the use of real
> primitives :
>
> vaccBubu <- function(age, female, ily) {
>     gender <- female * 1.25
>     gender[!female] <- 0.75
>     p <- (0.25 + 0.3 * 1 / (1 - exp(0.04 * age)) + 0.1 * ily) * gender
>     p[p < 0] <- 0
>     p[p > 1] <- 1
>     p
> }
>

I think you could do it a bit more elegantly with

gender <- 0.75 + female * 0.5

but pmin and pmax are "real" primitives, so I think it's fine to use them.

Hadley


-- 
RStudio / Rice University
http://had.co.nz/