MAT8186 -- Cours 3

fonctions.R

## Debug
fact <- function(x) {
    x < 0 && return(x)

    x * fact(x - 1.0)
}

## Profile
cutVector <- function(vec, m) {
    n <- length(vec) / m
    lapply(0:(m-1),
           \(x) vec[seq(n * x + 1, n * (x + 1))])
}

dist0 <- function(m, n = 1e4) {
    dat <- sample(0:9,
                  size = n * m,
                  replace = TRUE)
    dat_chopped <- cutVector(dat, m)
    sapply(dat_chopped, \(v) mean(v == 0))
}

## Benchmark
fact_for <- function(x) {
    acc = 1
    for (k in 1:x)
        acc <- acc * k
    acc
}