底辺私立医大を卒業した医者って頭悪いよね？ Part17

library(gtools)
library(fitdistrplus)
pm0=permutations(6,6,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) ; sd(d)
fit=fitdistrplus::fitdist(d,'norm')
plot(fit)
table(d)
plot(table(d))
plot(table(d)/S,ann=F)

#
x=1:6
n=length(x)
(m=mean(x))
(SS=sum((x-m)^2))
(var=var(1:6)) ; SS/(n-1) ; SS/n
(sd=sd(x))
N=100
(M=m*N)
(Var=var*N)
(Sd=sqrt(Var))
curve(dnorm(x,M,Sd),N,6*N)
pnorm(3*N,M,Sd,lower=T)