定義運算骰子總合之function

dice.sum <- function(n.dice){
  dice <- sample(1:6, size = n.dice, replace = T)
  return(sum(dice))
}

接著應以三個骰子可能出現的組合6x6x6=216個狀況重複擲骰,然機率問題可能會有遺漏一兩個數值,因此此處重複1000次

dice_3 <- replicate(1000, dice.sum(3))
plot(table(dice_3),ylab = 'Frequency', main = 'Three dice rolling')

END