rm(list=ls(all=TRUE))

# absolute Häufigkeiten von Aufgabe 2, Blatt 4
h <- matrix(c(665, 161, 1160, 103, 2361, 108), nrow = 3, byrow = TRUE)
n <- sum(h)

# Chi-Quadrat-Koeffizient berechnen
T <- 0
for (i in 1:2) {
	for (j in 1:3) {
		T <- T + (h[j, i] - sum(h[, i]) * sum(h[j, ]) / n)^2 / (sum(h[, i]) * sum(h[j, ]) / n)
	}
}

# Kontingenzkoeffizient berechnen
T2 <- sqrt(T / (n + T))

# korrigierten Kontingenzkoeffizient berechnen
T3 <- T2 / sqrt((2 - 1) / 2)