Binomujakauma
Todennäköisyysfunktio
======================
n<-20; x<-0:n; p<- 1/2
dbinom(x, n, p)
plot(dbinom(x, n, p))
plot(dbinom(x, n, p),type="h")
points(dbinom(x, n, p))

# Kertymäfunktiofunktio
y<-pbinom(x,n,p)
plot(x,y,type="s")

cy<-c(0,y)
Kf<-stepfun(x, cy, f = 0)
plot(Kf,verticals= FALSE)


# Ryhmätestaus
p<-0.0001
curve((x+1)/x -(1-p)^x,5,100,ylim=c(0,0.6))
p<-0.001
curve((x+1)/x -(1-p)^x,5,100,ylim=c(0,0.6),add=T)
p<-0.005
curve((x+1)/x -(1-p)^x,5,100,ylim=c(0,0.6),add=T)
p<-0.01
curve((x+1)/x -(1-p)^x,5,100,ylim=c(0,0.6),add=T)