zoo <- matrix(
c(
1969,1119000,1560
,1971,1250000,1623
,1972,1250000,1624
,1973,1300000,1810
,1974,1300000,1700
,1975,1332194,1654
,1976,1264726,1948
,1977,1218617,1407
,1978,1380000,1700
,1979,1238000,1324
,1980,1240000,1388
,1981,1130000,NA
,1982,1330000,1367
,1983,1246500,1316
,1984,1330000,1183
,1985,1243000,1007
,1986,1230000,998
,1987,1268000,944
,1988,1330000,924
,1989,1560000,1075
,1990,1021293,1580
,1991,789272,3254
,1992,901496,3790
,1993,831426,2270
,1994,761563,1979
,1995,782951,1954
,1996,686731,1501
,1997,687398,1382
,1998,744749,2147
,1999,755736,2829
,2000,751619,4227
,2001,1172158,6392
,2002,1164392,8546
,2003,1181535,11122
,2004,1319767,14000
,2005,1457938,25616
,2006,1759963,35554
,2007,1774587,42854
,2008,1673608,39515
,2009,1723349,45209
,2010,1560791,27553
,2011,2098075,38732
,2012,2123349, 34701
,2013,1842000, 39024
),
ncol=3, byrow = TRUE)
jahr = zoo[,1]
besucher = zoo[,2]
dauerticket = zoo[,3]
par(mar = c(3.2, 3, 3, 1))
fontscale=0.7
#Leerer Plot zwecks Einstellungen alle am selben Ort etc
plot( jahr, besucher, ylim=c(0,2200000), type="l", col=NA,
xlab=" ", ylab=" ", main="Zoo Leipzig, Besucherzahlen in Millionen nach Jahr", yaxt="n", cex.main=.8, cex.axis=fontscale)
y<-c(0.5, 1, 1.5, 2)*1e6
axis(2, at=y, labels=formatC(y/1e6, big.mark = ".", format = "fg", decimal.mark=","), cex.axis=fontscale)
#Als Flaechen
polyX = c(jahr, max(jahr), min(jahr))
polyBesucher = c(besucher, 0, 0)
polygon( polyX, polyBesucher, col="royalblue", border=NA)
#Ohne Fehlende Werte:
zooDauer = matrix( zoo[!is.na(zoo[,3])] ,ncol=3)
jahrDauer = zooDauer[,1]
dauerticket = zooDauer[,3]
polyXDauer = c( jahrDauer, max(jahrDauer), min(jahrDauer))
polyDauerticket = c(dauerticket, 0, 0)
polygon( polyXDauer, polyDauerticket, col="gold", border=NA)
lines( jahr, besucher, type="l", col="navy")
lines( jahr, zoo[,3], type="l", col="orange")
grid()