R/sum_dice.R

In tleblevecIMP/ICgeostats: Basic geostatistics

# Monte carlo simulation of the sum of dice

sum_dice<-function(n_dice,n_trials){

  # Preparation
  dice<-c(1,2,3,4,5,6)
  monte_carlo <-integer(length=n_trials)

  # Simulation
  for ( i in 1:n_trials){
    realization = sample(x=dice,size=n_dice,replace=TRUE)
    monte_carlo[i] = sum(realization)
  }
  hist(monte_carlo)
}