n <- 10; p <- 0.6; X = rbinom(1000,n,p); hist(X/n)

n <- 10

p <- 0.6; X <- rbinom(1000,n,p); hist(X/n,xlim=c(0,1),breaks=30)




library(MASS)

n<-2; for(i in (1:16) ){
                 X <- rbinom(1000,n,p); 
                 truehist(X/n,xlim=c(0,1),nbins="FD",prob=TRUE,ylim=c(0,20),main=paste("n = ",n)); 
                 locator(1) ; n<-2*n }

p <- 0.6; X <- rbinom(1000,1,p); M <- cumsum(X)/(1:1000); plot(M)





K <- 5
p <- 0.6; X <- rbinom(1000,10,p); verwerp <- sum(X>K) ;
cat("We verwerpen de nulhypothese in ",verwerp," uit 1000 keer")

K <- K+1
K <- K-1

p <- 0.5; X <- rbinom(1000,10,p); verwerp <- sum(X>K) ;