die1 <- 1:4
die1
die2 <- 1:4
die2
dieM <- die1*die2
dieRoll <- expand.grid(die1, die2)
dieRoll
dieRoll$Value <- dieRoll$Var1 * dieRoll$Var2
head(dieRoll,16)
roll <- function(){
die <- 1:4
dice <- sample(die, size = 2, replace = TRUE)
sum(dice)
}
roll()
roll()
roll()
Rolls <- replicate(10000, roll())
library("ggplot2")
qplot(Rolls, bindwidth = 1)
output: html_document: toc:yes fig_caption:yes number_sections:yes keep_md:yes