regionen <- c("Nord", "Sued", "Ost", "West")
Z <- normalmatrix(4711, 4, 24)
zeit <- 0:23
niveau <- c(120, 60, 95, 80)
trend <- c(0.0, 1.2, -0.15, 0.1)
daten <- niveau + trend %o% zeit + Z * 6 # 4 Zeilen, 24 Spalten
kennzahlen <- t(sapply(1:4, function(i) {
y <- daten[i, ]
m <- lm(y ~ zeit)
c(mittel = mean(y), sd = sd(y), steigung = unname(coef(m)[2]))
}))
rownames(kennzahlen) <- regionen
round(kennzahlen, 3) mittel sd steigung
Nord 121.126 5.197 0.084
Sued 73.634 9.747 1.258
Ost 95.110 6.436 -0.197
West 81.192 6.682 -0.208
round(c(groesste_region_mittel = max(kennzahlen[, "mittel"]),
anteil_sued_monat_1 = daten[2, 1] / sum(daten[, 1]) * 100,
anteil_sued_monat_24 = daten[2, 24] / sum(daten[, 24]) * 100), 1)groesste_region_mittel anteil_sued_monat_1 anteil_sued_monat_24
121.1 16.6 24.3
par(mfrow = c(1, 3), mar = c(4, 4, 3, 1))
# Aussage 1: Wer ist am groessten?
barplot(sort(kennzahlen[, "mittel"]), horiz = TRUE, las = 1, col = neutral,
border = NA, main = "1: Wer ist am grössten?", xlab = "Umsatz")
# Aussage 2: Wer waechst?
matplot(zeit, t(daten), type = "l", lty = 1, lwd = 1.6,
col = c(neutral, akzent, neutral, neutral),
xlab = "Monat", ylab = "Umsatz", main = "2: Wer wächst?")
# Aussage 3: Wie verschiebt sich der Anteil?
anteile <- prop.table(daten, margin = 2) * 100
plot(zeit, anteile[2, ], type = "l", lwd = 2, col = akzent, ylim = c(0, 40),
xlab = "Monat", ylab = "Anteil Süd in Prozent",
main = "3: Wie verschiebt sich der Anteil?")
par(mfrow = c(1, 1))






