#0:
# Experiment zum Uebergang von der geometrischen zur Exponentialverteilung
source(file="http://www.wiwi.uni-bielefeld.de/~wolf/software/R-wtools/slider/slider.R")

geo.to.exp<-function(){
  refresh.code<-function(...){
  # Vorbereitung
    n<-slider(no=1); p<-slider(no=2)/100; x.0<-slider(no=3)
  # Berechnung
     pn<-p/n; x<-1:(1/pn*3);xn<-x/n; lambda<-pn
     f.geo<-dgeom(x-1,pn); f.exp<-dexp(x-1,lambda)
     F.geo<-pgeom(x-1,pn); F.exp<-pexp(x-1,lambda)
     if(x.0>max(xn)) x.0<-slider(set.no.value=c(3,ceiling(max(xn))))
# Plot
    par(mfrow=c(2,1))
    plot(xn,f.geo,type="h",bty="l",xlim=c(0,max(xn)),
        xlab="Untersuchungszeitpunkte",ylab="f.geo | f.exp",
         main=paste("Anzahl Versuche, p=",round(p,2),
         ",\nTeilungen n=",n,", lambda/n=",round(lambda*n,2),sep=""))
    lines(xn-xn[1],f.exp,lty=2); abline(v=x.0,col="red",lty=2)

    plot(c(0,xn),c(0,F.geo),type="s",bty="l",
        xlab="Untersuchungszeitpunkte",ylab="F.geo | F.exp",
        xlim=c(0,max(xn)),ylim=0:1)
    lines(xn-xn[1],F.exp,lty=2); abline(v=x.0,col="red",lty=2)
    xt<-0.7*par()$usr[2]
    text(xt,0.4,pos=4,paste("P(X.exp<=",x.0,")=",
                              round(pexp(x.0,lambda),4),sep=""))
  par(mfrow=c(1,1))
  }
  slider(refresh.code,
       c("n: Anzahl Teilungen",
         "p: Event-WS in %",
          "x.0: x-Wert"),
       c(1,1,.02),c(25,99,300),c(1,1,.2),c(1,50,2)  )
}
#:0
#3:
geo.to.exp()
#:3