Datei Diskussion:Parlamentswahl in Polen 2019-Umfragen.svg

Code für R

library("Cairo")
CairoSVG("parlamentswahlPL-2019.svg",width=8.5,height=5.5)
daten3<-data.frame(zeit=c(25,1461),
"inst"=c("",""),
"ENDE"=c(50,50))
daten<-data.frame(zeit=c(), 
"inst"=c(), 
"PIS"=c(), 
"PO"=c(),
"NOWO"=c(), 
"ZL"=c(), 
"PSL"=c(),
"KORWIN"=c(),  
"RAZEM"=c()) 

blau<-"#073A76" #PIS
orange<-"#FCA241" #PO
schwarz<-"#000000" #KUKIZ
tuerkis<-"#005CA9" #NOWO
rot<-"#E2001A" #ZL
gruen<-"#1bb100" #PSL
dunkelblau<-"#E8B909" #KORWIN
violett<-"#660033" #RAZEM 

cexwert=0.5

D<-dim(daten)[1]

plot(daten3$zeit,daten3$'ENDE',pch=as.character(daten3$inst), col="#000000",ylim=c(0,max(daten3$'ENDE')),
cex=cexwert, axes=FALSE,xlab="",ylab="")
for (i in 1:50)lines(c(0,3000),i*c(1,1),col="lightgrey",lwd=.3)
for (i in 1:10)lines(c(0,3000),i*c(5,5),col="lightgrey",lwd=.6)
lines(c(0,3000),rep(0,2),col="black",lty=1,lwd=2)

axis(2,at=c(0,5,10,20,30,40,50),labels=c("0%","5%","10%","20%","30%","40%","50%"), tick=FALSE, las=2)

verschieben=2

axis(1,at=seq(25,122,365.25/3),labels=paste(rep(c("Okt."),1),"\n",rep(2015:2015,each=1),sep=""),tick=FALSE)
axis(1,at=seq(153,2600,365.25/3),labels=paste(rep(c("Feb.","Jun.","Okt."),7),"\n",rep(2016:2019,each=3),sep=""),tick=FALSE)
points(daten$zeit,daten$PIS,pch=as.character(daten$inst),col=blau,cex=cexwert)
points(daten$zeit,daten$PO,pch=as.character(daten$inst),col=orange,cex=cexwert)
points(daten$zeit,daten$KUKIZ,pch=as.character(daten$inst),col=schwarz,cex=cexwert)
points(daten$zeit,daten$NOWO,pch=as.character(daten$inst),col=tuerkis,cex=cexwert)
points(daten$zeit,daten$ZL,pch=as.character(daten$inst),col=rot,cex=cexwert)
points(daten$zeit,daten$PSL,pch=as.character(daten$inst),col=gruen,cex=cexwert)
points(daten$zeit,daten$KORWIN,pch=as.character(daten$inst),col=dunkelblau,cex=cexwert)
points(daten$zeit,daten$RAZEM,pch=as.character(daten$inst),col=violett,cex=cexwert)

lines(c(0,3000),rep(5,2),col="darkgrey",lty=2,lwd=1.5)
title(main="Parlamentswahl in Polen 2019",font=2)

points(min(daten$zeit),daten$"PIS"[1],col="#073A76",pch=19,cex=1.25)
points(min(daten$zeit),daten$PO[1],col=orange,pch=19,cex=1.25)
points(min(daten$zeit),daten$KUKIZ[1],col=schwarz,pch=19,cex=1.25)
points(min(daten$zeit),daten$NOWO[1],col=tuerkis,pch=19,cex=1.25)
points(min(daten$zeit),daten$ZL[1],col=rot,pch=19,cex=0.75)
points(min(daten$zeit),daten$PSL[1],col=gruen,pch=19,cex=1.25)
points(min(daten$zeit),daten$KORWIN[1],col=dunkelblau,pch=19,cex=0.75)
points(min(daten$zeit),daten$RAZEM[1],col=violett,pch=19,cex=0.75)

daten<-daten[order(daten$zeit),]

movingaverage<-function(x,n) {
  l <- length(x)-n+1
  y <- numeric(l)
  for (i in 1:l)   y[i]<-mean(x[(1:n)+i-1], na.rm=TRUE) 
  return(y) 
}

for (i in 3:10) {
 daten[7:D,i]<-movingaverage(daten[,i],7)
}

daten<-daten[-(1:6),]

lines(daten$zeit[-D],daten$PIS[-D],pch=as.character(daten$inst),col="#073A76",lwd=2) 
lines(daten$zeit[-D],daten$PO[-D],pch=as.character(daten$inst),col=orange,lwd=2) 
lines(daten$zeit[-D],daten$KUKIZ[-D],pch=as.character(daten$inst),col=schwarz,lwd=2) 
lines(daten$zeit[-D],daten$NOWO[-D],pch=as.character(daten$inst),col=tuerkis,lwd=2) 
lines(daten$zeit[-D],daten$ZL[-D],pch=as.character(daten$inst),col=rot,lwd=2)
lines(daten$zeit[-D],daten$PSL[-D],pch=as.character(daten$inst),col=gruen,lwd=2) 
lines(daten$zeit[-D],daten$KORWIN[-D],pch=as.character(daten$inst),col=dunkelblau,lwd=2) 
lines(daten$zeit[-D],daten$RAZEM[-D],pch=as.character(daten$inst),col=violett,lwd=2) 

dev.off()