#0:
erfolge.bei.bernoulli.experimenten<-function(){
  refresh.code<-function(...){
  # Modellierung
     prozent.ok<-slider(no=1); f.ok<-prozent.ok/100; h.ok<-21000*f.ok
     gg<-c(rep(1,h.ok),rep(0,21000-h.ok))
  # Experimentparameter
     n.stpr<-slider(no=2); wd<-slider(no=3); set.seed(zz<-slider(no=4))
  # Umsetzung
     result<-unlist(lapply(1:wd,function(x) { stpr<-sample(gg,size=n.stpr); sum(stpr) }))
  # Ergebnisdarstellung
     f.anz<-table(result)/wd
     plot(as.character(names(f.anz)),f.anz,bty="n",
        xlab="Erfolgsanzahl",ylab="relative H\344aufigkeit",
        xlim=c(0,n.stpr),ylim=c(0,1.5*dbinom(round(n.stpr*.5),n.stpr,.5)),
        main=paste("H\344ufigkeiten verschiedener Erfolgsanzahlen\n","% Erfolge in GG=",
                           prozent.ok,", n.stpr=",n.stpr,", WD=",wd,", Zufall=",zz,sep=""),type="h")
    result<<-result
  }
  slider(refresh.code,
       c("GG: Erfolgsprozentsatz","Stichprobenumfang","Experiment-Wiederholungen","Zufallsstart"),
       c(1,10,100,1),c(99,100,3000,30),c(1,5,200,1),c(58,50,500,13)  )
}

#:0
#1:
#Experiment zur Ermittlung empirischer Haeufigkeiten von Erfolgsanzahlen.
source(file="http://www.wiwi.uni-bielefeld.de/~wolf/software/R-wtools/slider/slider.R")

#:1

#3:
erfolge.bei.bernoulli.experimenten()
#:3