setwd('C:/Users/matthias/Lehre/StochastikI/Blatt11') data <- read.table('Bleiakkus.txt', header=TRUE) x1 <- data$Typ_A x2 <- data$Typ_B n1 <- length(x1) n2 <- length(x2) s2 <- sqrt((sum((x1 - mean(x1))^2) + sum((x2 - mean(x2))^2)) / (n1 + n2 - 2)) T1 <- sqrt((n1 * n2) / (n1 + n2)) * (mean(x1) - mean(x2)) / s2 alpha <- 0.05 c <- qt(1 - alpha/2, df = n1 + n2 - 2) print(T1) print(c) if(abs(T1) > c){ print('H0 wird verworfen') }else{ print('H0 wird nicht verworfen') } s21 <- var(x1) s22 <- var(x2) T2 <- s22 / s21 c1 <- qf(alpha/2, df1 = n2-1, df2 = n1 - 1) c2 <- qf(1 - alpha/2, df1 = n2-1, df2 = n1 - 1) if(T2 < c1 || T2 > c2){ print('H0 wird verworfen') }else{ print('H0 wird nicht verworfen') }