# Teil (a) empkvf <- function(x){ y <- y <- sort(x, index.return = TRUE) F <- rep(0, length(x)) for(i in 1:length(x)){ F[y$ix[i]] <- i } for(i in 1:length(x)){ F[i] = max(F[x == x[i]]) } F <- F / length(x) M <- matrix(c(x,F), ncol=2) return(M) } # Teil (b) plot_empkvf <- function(M, s){ L <- max(t) + 1 t <- sort(M[,1]) F <- sort(M[,2]) t2 <- c(t[2:length(t)], L) segments(t, F, t2, F, col=s, lwd=2) } range <- max(X)+1 t <- seq(0,range,0.01) X <- rexp(1000, 1) X1 <- X[1:10] X2 <- X[1:100] F0 <- 1 - exp(-t) plot(t, F0, col='red', type='l', xlim = c(0,range), lwd=2, main='Verteilungsfunktionen') plot_empkvf(empkvf(X1),'blue') plot_empkvf(empkvf(X2),'green') plot_empkvf(empkvf(X),'black')