setwd('C:/Users/matthias/Lehre/RaeumlicheStatistik/Blatt2') data <- read.table('Blatt2dat.txt') X <- data$X Y <- data$Y plot(X,Y) c <- sum(log(X)) theta1 <- -N / c - 1 print(theta1)