[R] Fisherman simulation (Poisson process)

fisherman-simulation.r

number.of.fish <- function() {
  k <- rpois(1, 1.2) # Use lambda = 1.2 because we care about the first 2 hours.
  ifelse(k > 0, k, 1)
}

simulate <- function(simulation.count) {
  if (simulation.count > 0)
    mean(sapply(1:simulation.count, function(i) number.of.fish()))
}

set.seed(13)
simulate(100000) # Outputs 1.50563