# 2009 
umsatz09 <- data.frame(
  jahr=c(1999:2009), 
  singles=c(56.9, 54.0, 51.8, 38.9, 26.8, 23.3, 17.8, 15.9, 10.7, 7.6, 5.7), 
  CD.alben=c(209.7, 206.1, 184.6, 178.7, 146.8, 145.5, 147.6, 149.5, 148.6, 145.1, 147.3), 
  mobile.audio=c(NA, NA, NA, NA, NA, 7.5, 19.7, 29.2, 39.9, 41.3, 49.2), 
  euro=c(2648, 2630, 2365, 2201, 1816, 1753, 1748, 1706, 1652, 1582, 1530) # gesamtumsatz phonomarkt deutschland in mio
)

raub09 <- data.frame(
  jahr=c(2003:2009),
  p2p=c(602, 383, 412, 374, 312, 316, 258), 
  zivilverfahren=c(NA, 20, 395, 2811, 15699, 13562, 18000)
)

# funktion zum z-transformieren:
# von jedem wert mittelwert abziehen und durch standardabweichung teilen
z.mean.sd <- function(data){
 return.values <- (data - mean(data, na.rm=TRUE)) / (sd(data, na.rm=TRUE)) 
 return(return.values)
}

## plots erstellen
plot(z.mean.sd(raub09$p2p),
 type="b",
 main=paste("P2P-Nutzung, CD-Verkäufe, Umsatz und Zivilverfahren\n(z-Standardisiert)"),
 sub="Zahlen: Jahreswirtschaftsberichte 2008/2009 Bundesverband Musikindustrie",
 xlab="",
 ylab="z-Wert",
 xaxt="n",
 lwd=2,
 col="red")
# x-achse mit jahreszahlen beschriften
axis(1, at=1:7, labels=raub09$jahr)
lines(z.mean.sd(raub09$zivilverfahren), type="b", col="brown", lwd=2)
lines(z.mean.sd(umsatz09$CD.alben[5:11]), type="b", col="black", lwd=2)
lines(z.mean.sd(umsatz09$euro[5:11]), type="b", col="blue", lwd=2)
legend("topright",
 legend=c("illegale Downloads","CD-Verkäufe","Umsatz","Zivilverfahren"),
 text.col=c("red","black","blue","brown"),
 col=c("red","black","blue","brown"),
 lty=c(1,1,1),
 lwd=3)

# absolute zahlen
plot(umsatz09$CD.alben[5:11],
 type="b",
 ylim=c(100,600),
 main="Absatzentwicklung 2003-2009",
 sub="Zahlen: Jahreswirtschaftsberichte 2008/2009 Bundesverband Musikindustrie",
 xlab="",
 ylab="Absatz/Downloads in Mio.",
 xaxt="n",
 lwd=2)
lines(raub09$p2p, type="b", col="red", lwd=2)
# x-achse mit jahreszahlen beschriften
axis(1, at=1:7, labels=umsatz09$jahr[5:11])
legend("topright",
 legend=c("illegale Downloads","CD-Verkäufe"),
 text.col=c("red","black"),
 col=c("red","black"),
 lty=c(1,1,1),
 lwd=3)

# alben a 10 tracks statt songs
plot(umsatz09$CD.alben[5:11],
 type="b",
 ylim=c(0,180),
 main=paste("Korrigierte Absatzentwicklung 2003-2009\n(Näherung auf 10 Songs als 1 Album)"),
 sub="Zahlen: Jahreswirtschaftsberichte 2008/2009 Bundesverband Musikindustrie",
 xlab="",
 ylab="Absatz/Downloads in Mio.",
 xaxt="n",
 lwd=2)
lines(raub09$p2p/10, type="b", col="red", lwd=2)
# x-achse mit jahreszahlen beschriften
axis(1, at=1:7, labels=umsatz09$jahr[5:11])
legend("topright",
 legend=c("illegale Downloads (Alben)","CD-Verkäufe (Alben)"),
 text.col=c("red","black"),
 col=c("red","black"),
 lty=c(1,1,1),
 lwd=3)

## zahlen vor 2003
umsatz04.schaetz <- data.frame(
  jahr=c(1995:2009),
  CD.alben=c(190.0, 198.2, 209.2, 206.8, 209.7, 206.1, 184.6, 178.7, 146.8, 145.5, 147.6, 149.5, 148.6, 145.1, 147.3)
)

# plotten
plot(umsatz04.schaetz$CD.alben,
 type="b",
 ylim=c(100,250),
 main=paste("CD-Verkäufe 1995-2009\n(1995/1996 geschätzt, da anders berechnet)"),
 sub="Zahlen: Jahreswirtschaftsberichte 2004-2009 Bundesverband Musikindustrie",
 xlab="",
 ylab="Absatz in Mio.",
 xaxt="n",
 lwd=2)
# x-achse mit jahreszahlen beschriften
axis(1, at=1:15, labels=umsatz04.schaetz[["jahr"]])
abline(v=4, lty=4, col="red", lwd=3)
text(x=4, y=240, labels="Napster", pos=4, col="red")
abline(v=6, lty=4, col="violet", lwd=3)
text(x=6, y=220, labels="\"Kopierschutz\"", pos=4, col="violet")

# CD-singles und downloads
plot(umsatz09$mobile.audio[5:11],
 type="b",
 col="orange",
 ylim=c(0,80),
 main=paste("Absatzentwicklung 2003-2009\n(CD-Singles und bezahlte Downloads)"),
 sub="Zahlen: Jahreswirtschaftsbericht 2009 Bundesverband Musikindustrie",
 xlab="",
 ylab="Absatz/Downloads in Mio.",
 xaxt="n",
 lwd=2)
axis(1, at=1:7, labels=umsatz09$jahr[5:11])
legend("topright", legend=c("CD-Singles","Downloads (Einzeltrack)"), text.col=c("black","orange"), col=c("black","orange"), lty=c(1,1,1), lwd=3)
lines(umsatz09$singles[5:11], type="b", lwd=2)


## statistik -- korrelationen
cor.test(raub09$p2p, raub09$zivilverfahren)
cor.test(raub09$p2p, umsatz09$CD.alben[5:11])
cor.test(raub09$p2p, umsatz09$euro[5:11])

# zahlen von 2008
umsatz08 <- data.frame(
jahr=c(1999:2008),
single=c(56.9, 54.0, 51.8, 38.9, 26.8, 23.3, 17.8, 15.9, 10.7, 7.6), 
CD.alben=c(209.7, 206.1, 184.6, 178.7, 146.8, 145.5, 147.6, 149.5, 148.6, 145.1), 
MC=c(21.5, 20.5, 22.3, 14.3, 15.5, 13.2, 8.7, 5.8, 4.6, 3.2), 
LP.vinyl=c(0.6, 0.9, 1.1, 1.1, 1.1, 0.8, 0.7, 0.6, 0.7, 0.9), 
DVD.audio=c(0.0, 0.0, 0.1, 0.2, 0.6, 0.3, 0.5, 0.2, 0.2, 0.4), 
sum.longplay=c(231.8, 227.5, 208.1, 194.3, 164.0, 159.8, 157.5, 156.1, 154.1, 149.6), 
DVD.musikvideo=c(0.0, 0.6, 1.5, 3.8, 9.0, 11.5, 12.9, 14.2, 13.6, 11.7), 
VHS.musikvideo=c(1.0, 0.6, 0.6, 3.2, 1.9, 1.2, 0.8, 0.1, 0.0, 0.0), 
sum.physisch=c(289.7, 282.7, 262.0, 240.2, 201.7, 195.8, 189.0, 186.3, 178.4, 168.9), 
mobile.audio=c(NA, NA, NA, NA, NA, 7.5, 19.7, 29.2, 39.9, 44.6), 
bundles=c(NA, NA, NA, NA, NA, 0.4, 1.4, 1.9, 2.6, 3.9), 
sum.download=c(NA, NA, NA, NA, NA, 7.9, 21.1, 31.1, 42.5, 48.5), 
klingelton=c(NA, NA, NA, NA, NA, NA, NA, 11.9, 8.3, 4.1), 
ring.back=c(NA, NA, NA, NA, NA, NA, NA, 1.3, 1.2, 1.4), 
sum.mobile=c(NA, NA, NA, NA, NA, NA, NA, 13.2, 9.5, 5.5), 
sum.digital=c(NA, NA, NA, NA, NA, 7.9, 21.1, 44.3, 52.0, 54.0), 
sum=c(289.7,282.7,262.0,240.2,201.7,210.1,230.6,230.4,203.7,222.9)
)