svg("Wanderungsbilanz Deutschland.svg", width = 15, height = 9)
jahr <- seq(1901, 2022)
wanderungsbilanz <- data.frame(jahr)
########################################
rownames( wanderungsbilanz) <- jahr
wanderungsbilanz[91:122, "zuwanderung"] <- c(1198978, 1502198, 1277408, 1082553, 1096048, 959691, 840633, 802456, 874023, 841158, 879217, 842543, 768975, 780175, 707352, 661855, 680766, 682146, 721014, 798282, 958299, 1080936, 1226493, 1464724, 2136954, 1865122, 1550721, 1585112, 1558612, 1186702, 1323466, 2672008)
comment( wanderungsbilanz$zuwanderung) <- "Sources: 1991- Statistisches Bundesamt ::Wanderungen zwischen Deutschland und dem Ausland, Zugezogene, Fortgezogene und Saldo::"
########################################
wanderungsbilanz[91:122, "abwanderung"] <- c(596455, 720127, 815312, 767555, 698113, 677494, 746969, 755358, 672048, 674038, 606494, 623255, 626330, 697632, 628399, 639064, 636854, 767889, 733796, 670605, 678969, 711991, 797886, 914241, 997552, 1365178, 1134641, 1185432, 1231552, 996451, 994303, 1217431)
comment( wanderungsbilanz$abwanderung) <- "Sources: 1991- Statistisches Bundesamt ::Wanderungen zwischen Deutschland und dem Ausland, Zugezogene, Fortgezogene und Saldo::"
########################################
wanderungsbilanz[91:122, "wanderungssaldo"] <- c(602523, 782071, 462096, 314998, 397935, 282197, 93664, 47098, 201975, 167120, 272723, 219288, 142645, 82543, 78953, 22791, 43912, -55743, -12782, 127677, 279330, 368945, 428607, 550483, 1139402, 499944, 416080, 399680, 327060, 220251, 329163, 1454577)
comment( wanderungsbilanz$wanderungssaldo) <- "Sources: 1991- Statistisches Bundesamt ::Wanderungen zwischen Deutschland und dem Ausland, Zugezogene, Fortgezogene und Saldo::"
########################################
jahre <- 1991:2022
par(mar = c(4.5, 6.5, 4.0, 0)) #bottom, left, top, right
plot(jahre, wanderungsbilanz$zuwanderung[91:122], ylim = c(-250000, 2700000), bty = "l", type = "b", yaxs = "i", xlab = "", ylab = "", pch = 16, cex = 2, lab = c(10, 10, 1), col = "orange1", main = "Wanderungsbilanz\n(1991-2022)", cex.main = 2, sub = "", lwd = 2, yaxt = "n", xaxt = "n", lty = 1)
lines(jahre, wanderungsbilanz$wanderungssaldo[91:122], lty = 1, type = "b", pch = 15, cex = 2, lwd = 2, col = "blue")
lines(jahre, wanderungsbilanz$abwanderung[91:122], lty = 1, type = "b", pch = 17, cex = 2, col = "limegreen", lwd = 2)
axis(1, at = seq(1991, 2022, by = 1), las = 2, tck = 1, lty = 2, col = "#c9c9c9", lwd = 1, cex.axis = 1.5)
axis(2, at = seq(-250000, 2700000, by = 250000), las = 2, tck = 1, lty = 2, col = "#c9c9c9", lwd = 1, cex.axis = 1.5)
axis(2, at = seq(0, 0, by = 1), tck = 1, las = 2, lty = 2, col = "#000000", cex.axis = 1.5)
legend(x = "topleft", # Position
inset = 0.02,
legend = c("Zuwanderungen", "Abwanderungen", "Wanderungssaldo"), # Legend texts
lty = c(1, 1, 1), # Line types
pch = c(16, 17, 15),
col = c("orange1", "limegreen", "blue"), # Line colors
lwd = 2,
cex = 1.7)
dev.off()