Code necessary to generate plots illustrating the evolution of playing time. Writeups for this plot can be found at http://www.massivereportdata.com/blog/2015/10/20/updating-the-evolution-of-playing-time/ and http://www.massivereportdata.com/blog/2015/09/26/plotting-the-evolution-of-playing-time/

evolution_plot.R

# Read in teams
team11 <- read.csv("minutes_evolution_11_2015.csv")
team12 <- read.csv("minutes_evolution_12_2015.csv")
team13 <- read.csv("minutes_evolution_13_2015.csv")
team14 <- read.csv("minutes_evolution_14_2015.csv")
team15 <- read.csv("minutes_evolution_15_2015.csv")
team16 <- read.csv("minutes_evolution_16_2015.csv")
team17 <- read.csv("minutes_evolution_17_2015.csv")
team18 <- read.csv("minutes_evolution_18_2015.csv")
team19 <- read.csv("minutes_evolution_19_2015.csv")
team20 <- read.csv("minutes_evolution_20_2015.csv")
team42 <- read.csv("minutes_evolution_42_2015.csv")
team43 <- read.csv("minutes_evolution_43_2015.csv")
team44 <- read.csv("minutes_evolution_44_2015.csv")
team45 <- read.csv("minutes_evolution_45_2015.csv")
team340 <- read.csv("minutes_evolution_340_2015.csv")
team427 <- read.csv("minutes_evolution_427_2015.csv")
team463 <- read.csv("minutes_evolution_463_2015.csv")
team479 <- read.csv("minutes_evolution_479_2015.csv")
team506 <- read.csv("minutes_evolution_506_2015.csv")
team547 <- read.csv("minutes_evolution_547_2015.csv")
colors=c("#e41a1c","#377eb8","#4daf4a","#984ea3","#ff7f00","#a65628","#f781bf")

# Define plot
plotPTEvolution <- function(team,name,representative,gameCount) {
  par(new=F)
  separator = 1
  j = 1
  for (i in names(team)) {
    if (j==representative) {
      plot(team[[i]],ylim=c(0,1),type="l",xlab="Game number",ylab="Percentage of playing time",main=name,col=colors[j%%6+1],pch=19,las=1,mai=c(1,2,1,8),xlim=c(1,36),cex.axis=0.7,cex.lab=0.85,yaxt="n")
      axis(2,at=pretty(team[[i]]), lab=paste0(pretty(team[[i]]) * 100,"%"), las=TRUE, cex.axis=0.7)
    } else {
      plot(team[[i]],ylim=c(0,1),type="l",xlab="",ylab="",main="",col=colors[j%%6+1],pch=19,mai=c(0,2,1,8),axes=FALSE,mar=c(1,2,1,8),xlim=c(1,36))
    }
    text(gameCount+separator,team[[gameCount,i]],names(team[i]),adj=c(0,0.5),cex=0.5)
    par(new=T)
    j = j + 1
  }
  return(TRUE)
}

# Example plots for Wil Trapp and Emanuel Pogatetz
par(mfrow=c(1,2))
plot(team11$Trapp,ylim=c(0,1),type="l",xlab="Game number",ylab="Percentage of playing time",main="Wil Trapp",col="red",pch=19,las=1,mai=c(1,2,1,8),xlim=c(1,36),cex.axis=0.7,cex.lab=0.85,yaxt="n")
plot(team11$Pogatetz,ylim=c(0,1),type="l",xlab="Game number",ylab="Percentage of playing time",main="Emanuel Pogatetz",col="red",pch=19,las=1,mai=c(1,2,1,8),xlim=c(1,36),cex.axis=0.7,cex.lab=0.85,yaxt="n")

# Choose for single-team plots
par(mfrow=c(1,1))

# Choose for an overall league plot with all teams
par(mfrow=c(4,5))

# League overview
plotPTEvolution(team13,"Chicago Fire", 2, 33)
plotPTEvolution(team14,"Colorado Rapids", 6, 32)
plotPTEvolution(team11,"Columbus Crew", 24, 33)
plotPTEvolution(team12,"DC United", 4, 33)
plotPTEvolution(team16,"FC Dallas", 19, 33)
plotPTEvolution(team427,"Houston Dynamo", 12, 33)
plotPTEvolution(team19,"Los Angeles Galaxy", 24, 33)
plotPTEvolution(team45,"Montreal Impact", 5, 33)
plotPTEvolution(team15,"New England Revolution", 16, 33)
plotPTEvolution(team547,"New York City", 29, 33)
plotPTEvolution(team20,"New York Red Bulls", 10, 33)
plotPTEvolution(team506,"Orlando City", 5, 33)
plotPTEvolution(team479,"Philadelphia Union", 18, 33)
plotPTEvolution(team42,"Portland Timbers", 6, 33)
plotPTEvolution(team340,"Real Salt Lake", 24, 33)
plotPTEvolution(team17,"San Jose Earthquakes", 11, 33)
plotPTEvolution(team43,"Seattle Sounders", 10, 33)
plotPTEvolution(team18,"Sporting Kansas City", 2, 32)
plotPTEvolution(team463,"Toronto FC", 4, 33)
plotPTEvolution(team44,"Vancouver Whitecaps", 1, 33)

# Game-specific preview
par(mfrow=c(1,2))
plotPTEvolution(team11,"Columbus Crew", 24, 33)
plotPTEvolution(team12,"DC United", 4, 33)

minutes_evolution_11_2015.csv