dice

サイコロ３個ふって出目をカウント
６進法を使ってみた。

library(oro.dicom)

## oro.dicom: Rigorous - DICOM Input / Output (version = 0.3.6)

dice <- function(dec = 2, size = 2) (sum(mapply(as.integer, strsplit(dec2base(dec, 
    6, size), "")) + 1))
a <- rep(0, 18)
names(a) <- 1:18

for (i in 0:215) {
    b <- dice(i, size = 3)
    a[b] <- a[b] + 1
}
a

##  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 
##  0  0  1  3  6 10 15 21 25 27 27 25 21 15 10  6  3  1 

round(a[3:18]/sum(a), 3)

##     3     4     5     6     7     8     9    10    11    12    13    14 
## 0.005 0.014 0.028 0.046 0.069 0.097 0.116 0.125 0.125 0.116 0.097 0.069 
##    15    16    17    18 
## 0.046 0.028 0.014 0.005 

barplot(a)

サイコロ10個振ると？
裏RjpWiki ダメ出し：方針が間違えている
裏RjpWiki ダメ出し：方針が間違えている　その2