Z <- normalmatrix(5151, 8, 24)
zeit <- 0:23
niveau <- seq(40, 110, length.out = 8)
reihen <- niveau + Z * 5
ziel <- 4 # die Serie, um die es geht
reihen[ziel, ] <- reihen[ziel, ] + seq(0, 25, length.out = 24)
bunt <- c("#1f77b4", "#ff7f0e", "#2ca02c", "#d62728",
"#9467bd", "#8c564b", "#e377c2", "#17becf")
mit_akzent <- rep(neutral, 8); mit_akzent[ziel] <- akzent
dicke_bunt <- rep(1.5, 8)
dicke_akzent <- rep(1.0, 8); dicke_akzent[ziel] <- 2.5
pfadlaenge <- function(y) {
dx <- diff(zeit) / (max(zeit) - min(zeit))
dy <- diff(y) / (max(reihen) - min(reihen))
sum(sqrt(dx^2 + dy^2))
}
gewichte <- function(pal, dicke) {
sapply(1:8, function(i) pfadlaenge(reihen[i, ]) * dicke[i] *
kontrast(pal[i], "#FFFFFF"))
}
for (nm in c("alles bunt", "ein Akzent")) {
g <- if (nm == "alles bunt") gewichte(bunt, dicke_bunt) else
gewichte(mit_akzent, dicke_akzent)
cat(sprintf("%-12s Zielserie %6.2f, Mittel der übrigen %5.2f, Verhältnis %5.2f\n",
nm, g[ziel], mean(g[-ziel]), g[ziel] / mean(g[-ziel])))
}alles bunt Zielserie 14.66, Mittel der übrigen 10.12, Verhältnis 1.45
ein Akzent Zielserie 33.40, Mittel der übrigen 3.13, Verhältnis 10.66
round(c(pfadlaenge_ziel = pfadlaenge(reihen[ziel, ]),
pfadlaenge_uebrige = mean(sapply((1:8)[-ziel],
function(i) pfadlaenge(reihen[i, ])))), 3) pfadlaenge_ziel pfadlaenge_uebrige
1.946 1.854
par(mfrow = c(1, 2), mar = c(4, 4, 3, 1))
matplot(zeit, t(reihen), type = "l", lty = 1, lwd = dicke_bunt, col = bunt,
xlab = "Monat", ylab = "Wert", main = "alles bunt")
matplot(zeit, t(reihen), type = "l", lty = 1, lwd = dicke_akzent,
col = mit_akzent, xlab = "Monat", ylab = "Wert", main = "ein Akzent")
par(mfrow = c(1, 1))








