# Aufgabe 3 t <- read.table("temperatur.txt", col.names="temp") quant <- 1/61*1:60 # k/(n+1) par(mfrow=c(2,2)) # Exponentialverteilung plot(-1/.05 * log(1-quant), sort(t$temp), main = "QQ-Plot, Exponentialverteilug", xlab = "Quantile der Exponentialverteilung", ylab = "empirische Quantile") abline(0,1) # Gleichverteilung plot(20+2*quant, sort(t$temp), main = "QQ-Plot, Gleichverteilug", xlab = "Quantile der Gleichverteilug", ylab = "empirische Quantile") abline(0,1) # Normalverteilung plot(qnorm(quant, 21, 1), sort(t$temp), main = "QQ-Plot, Normalverteilung", xlab = "Quantile der Normalverteilung", ylab = "empirische Quantile") abline(0,1) # Weibull-Verteilung plot((-21 * log(1-quant)), sort(t$temp), main = "QQ-Plot, Weibull-Verteilug", xlab = "Quantile der Weibull-Verteilug", ylab = "empirische Quantile") abline(0,1)