rm(list=ls())
source('toolmini.R')
library(gtools)
library(fitdistrplus)
N=6
(pm0=permutations(6,N,rep=T))
f<-function(x) sum(x) > 26
s=sum(apply(pm0,1,f))
S=nrow(pm0)
cat(s,'/',S,'\n')
s/S

d=apply(pm0,1,sum)
hist(d)
mean(d) ; var(d)
N*mean(1:6) ; N*sum((1:6-mean(1:6))^2)/N
fit=fitdistrplus::fitdist(d,'norm')
plot(fit)
table(d)
plot(table(d))
plot(table(d)/S,ann=F)