Calculating College Basketball rankings using functional programming in R

March Madness is officially upon us as College Basketball teams across the US try to get a few more signature wins on their resume during conference tournaments before Selection Sunday. With only 30 or so games in a given season and around 350 teams in Division I, being able to measure the relative strength of teams who have not played one another becomes an important challenge.

Unsurprisingly, there are literally hundreds of different ranking systems for College Basketball and other sports, some created by large news organizations and others by passionate fans. The focus here today is to describe two popular methods, the Rating Percentage Index (RPI) and the Simple Rating System (SRS), that you can easily calculate and track yourself using functional programming in R. Then, in future posts to be released throughout the tournament season over the next few weeks, I will be able to use the rating systems described here to answer some interesting questions about tournament teams in the past and what is happening in 2018.

For this, we will need the following R packages, and a small dataset outlining an example season.

library(readr) # for reading in data
library(dplyr) # for data cleaning
library(purrr) # for functional programming
library(forcats) # for recoding factors
library(limSolve) # for solving linear equations

# Read in data example

ex_data <- read_csv("../data/SmallExampleResults.csv")
ex_data

## # A tibble: 6 x 5
##   WTeamID WScore   LTeamID LScore  WLoc
##     <chr>  <int>     <chr>  <int> <chr>
## 1   UConn     64    Kansas     57     H
## 2   UConn     82      Duke     68     H
## 3   UConn     72 Wisconsin     71     A
## 4  Kansas     69     UConn     62     H
## 5    Duke     81 Wisconsin     70     H
## 6  Kansas     62 Wisconsin     52     A

Ratings Percentage Index (RPI)

The RPI is one of the most widely known (and widely criticized) rating systems used in College Basketball. It is a weighted formula, where 25% comes from a team’s own winning percentage (WP), 50% from its opponents’ winning percentage (OWP), and 25% from its opponents’ opponents’ winning percentage (OOWP). In other words, it is trying to use a team’s win/loss record and their strength of schedule to infer their rating. The main criticisms of the RPI is that it fails to take into account margin of victory and it places too much weight on facing strong opponents (to the point where it clearly rewards teams for losing badly to many strong opponents and penalizes teams for playing well against weaker ones). That being said, it is still one of the biggest factors used by the committee for tournament selection and seeding, so it’s worth looking at and understanding.

Though calculating the RPI seems simple on the surface, there are two additional caveats that make it a bit more challenging (and a great example for functional programming).

1. An additional weighting factor (1.4 for away wins/losses, 0.6 for home wins/losses ) is included for the reference team’s WP calculation, but not for OWP and OOWP

2. For OWP and OOWP, the reference team is not included in the win/loss record, making both of these calculations for the reference team unique to that team. In other words (and using the example data above), Kansas as an opponent will have a different WP for when calculating the OWP for UConn than it would for Wisconsin. If you are still confused, I highly recommend checking out the example here, which walks you through the calculation step-by-step.

So putting it all together, our R implementation requires the creation of three main functions, one each for WP, OWP, and OOWP. Below is the WP function, which is able to calculation a team’s winning percentage in two different ways. The first way occurs if a team should be excluded from the calculation, which will come in handy for the OWP and OOWP functions. The second is the weighted formula for home and away wins/losses, which is relevant for the WP part of the RPI formula for a given team. Here, we see that while UConn’s WP is 75% (3/4), its weighted WP is actually 81% (2.6/3.2).

# Team winning percentage

calc_wp <- function(game_data, team_id, exclusion_id = NULL){
  
  games_played <- game_data[game_data$WTeamID == team_id | game_data$LTeamID == team_id, ]

  if(!is.null(exclusion_id)){
    
    games_played <- 
      games_played[games_played$WTeamID != exclusion_id & games_played$LTeamID != exclusion_id, ]
    
    wp <- sum(games_played$WTeamID == team_id)/length(games_played$WTeamID)
    
  } else{
    
    wwins <- 1.4 * sum(games_played$WTeamID == team_id & games_played$WLoc == "A") +
      0.6 * sum(games_played$WTeamID == team_id & games_played$WLoc == "H") +
      sum(games_played$WTeamID == team_id & games_played$WLoc == "N")
    
    wlosses <- 1.4 * sum(games_played$LTeamID == team_id & games_played$WLoc == "A") +
      0.6 * sum(games_played$LTeamID == team_id & games_played$WLoc == "H") +
      sum(games_played$LTeamID == team_id & games_played$WLoc == "N")
    
    wp <- wwins/(wwins + wlosses)
    
  }

  return(wp)
  
}

calc_wp(ex_data, team_id = "UConn")

## [1] 0.8125

Here we have the function to calculate the second part of the RPI formula, OWP. In order to perform the calculation, we need to calculate the winning percentage of each opponent on a given team’s schedule. This is where functional programming comes in, which allows us to apply the function calc_wp element-wise to a vector of opponents. Here, I am using the map family of functions (specifically map_dbl, as I am returning a double) in the purrr package, which I find to be more user-friendly than the apply family equivalents. If you are unfamiliar, I highly recommend Jenny Bryan’s tutorials on the subject. At any rate, once I get the WP for each opponent, I simply return the mean. In this case, UConn’s OWP is 75%.

# Opponents winning percentage (exclusive of the reference team)

calc_owp <- function(game_data, team_id){
  
  opp_games <- game_data[game_data$WTeamID == team_id | game_data$LTeamID == team_id, ]
  opps <- if_else(opp_games$WTeamID == team_id, opp_games$LTeamID, opp_games$WTeamID)
  
  owp <- opps %>%
    map_dbl(~ calc_wp(game_data, team_id = .x, exclusion_id = team_id))
  
  return(mean(owp))
  
}

calc_owp(ex_data, team_id = "UConn")

## [1] 0.75

Finally, the OOWP calculation is actually quite simple to compute now that the calc_owp function has been created. It is identical in structure to that function, except instead of calling calc_wp on the vector of opponents, we call calc_owp on that vector. This returns the OWP for the vector of opponents, which is the OOWP. For the case of UConn, their OOWP is 51%.

# Opponents opponents winning percentage

calc_oowp <- function(game_data, team_id){
  
  opp_games <- game_data[game_data$WTeamID == team_id | game_data$LTeamID == team_id, ]
  opps <- if_else(opp_games$WTeamID == team_id, opp_games$LTeamID, opp_games$WTeamID)
  
  oowp <- opps %>%
    map_dbl(~ calc_owp(game_data, team_id = .x))
  
  return(mean(oowp))
  
}

calc_oowp(ex_data, team_id = "UConn")

## [1] 0.5138889

Putting it all together, we just need a final function that creates the linear combination of WP, OWP, and OOWP defined above. For UConn, this value is 0.7066.

# RPI using weighted formula

calc_rpi <- function(game_data, team_id){
  
  rpi <- 0.25 * calc_wp(game_data, team_id) +
    0.5 * calc_owp(game_data, team_id) +
    0.25 * calc_oowp(game_data, team_id)
  
  return(round(rpi, 4))
 
}

calc_rpi(ex_data, team_id = "UConn")

## [1] 0.7066

Again using the map_dbl function, we can easily apply our RPI function to all teams we have in our dataset and get an idea of how teams relate to one another. In this simple example, we see that UConn and Kansas are close in rating (which makes sense given they went 1-1 against the other, with both wins for each team coming at home). There is also a clear difference between the top two teams and the bottom two.

ex_teams <- unique(c(ex_data$WTeamID, ex_data$LTeamID))

data_frame(Team = ex_teams,
           RPI = map_dbl(ex_teams, ~ calc_rpi(ex_data, team_id = .x)))

## # A tibble: 4 x 2
##        Team    RPI
##       <chr>  <dbl>
## 1     UConn 0.7066
## 2    Kansas 0.6830
## 3      Duke 0.4340
## 4 Wisconsin 0.3403

The Simple Rating System (SRS)

While the RPI is fairly straightforward to understand, the relative scores between teams have very little meaning. Because it does not take margin of victory into account, it is hard to say how teams would actually match up to one another, and how that match up might differ on a home, away, or neutral court. This is where a system that explicitly takes margin of victory into account comes in handy, and the SRS is one of these systems.

The main idea is that a team’s strength is comprised of its margin of victory in games and strength of schedule derived from the margin of victory of teams in other games. It also allows for the estimation of a global parameter that yields the effect of where the game was played (i.e., home vs. away). The key benefit of this system is that it is both interpretable and easy to estimate. Because of this, it is a common tool used in sports reporting, like in some of the articles by 538.

That being said, there is no such thing as a free lunch, and the SRS’s simplicity is not without drawbacks. Namely, it ignores wins/losses, is slightly biased to offensive-oriented and faster pace teams where larger victories are more common, and it weights all games in a schedule equally.

The calculation of the SRS relies on three matrices. The first, G(ames), is an MxN+1 matrix where M is each game in a season, N are the teams that played in the season (with an additional column for home/away designation), and each element has a 1 to indicate the winner, a -1 to indicate the loser, and a 0 if the teams were not involved in that game. The additional column uses 1/-1 to indicate a game was home/away. The second, R(atings), is a N+1×1 matrix of rating scores to be estimated. The third, S(scores), is an Mx1 matrix with the margin of victory for each game. Then, we just set up the equation GR = S and solve for R using a matrix solver (for more information on the math behind the SRS, I would recommend this resource).

With this in mind, we can quickly use functional programming again to iterate over all the teams in our dataset and use a custom function, transform_wl, to create a series of column vectors with the relevant win/loss information to comprise our G matrix.

all_teams <- unique(c(ex_data$WTeamID, ex_data$LTeamID))

# Function to create column vector of wins/losses for each team in every game

transform_wl <- function(game_data, team_id){
  
  col_w <- if_else(game_data$WTeamID == team_id, 1, 0) %>%
    na_if(0)
  
  col_l <- if_else(game_data$LTeamID == team_id, -1, 0) %>%
    na_if(0)
  
  col_all <- coalesce(col_w, col_l) %>%
    tbl_df()
  
  return(col_all)
  
}

# Replace NAs with 0 and cbind home/away column

srs_ex <- map(all_teams, ~ transform_wl(ex_data, team_id = .x)) %>%
  bind_cols() %>%
  setNames(all_teams) %>%
  replace(is.na(.), 0) %>%
  mutate(loc = fct_recode(ex_data$WLoc, "1" = "H", "-1" = "A", "0" = "N")) %>%
  mutate(loc = as.numeric(as.character(loc))) %>%
  select(loc, everything()) %>%
  as.matrix()

srs_ex

##      loc UConn Kansas Duke Wisconsin
## [1,]   1     1     -1    0         0
## [2,]   1     1      0   -1         0
## [3,]  -1     1      0    0        -1
## [4,]   1    -1      1    0         0
## [5,]   1     0      0    1        -1
## [6,]  -1     0      1    0        -1

The matrix S is much easier to create with a few simple steps, then we can use the lsei function in the limSolve package to estimate the parameters of interest.

In interpreting the results of the SRS, everything can be related to the expected performance relative to an average team on a neutral court. This time, we see that Kansas is actually rated higher than UConn, likely due to them having a much better margin of victory against Wisconsin, the weakest team in the schedule. Additionally, there seems to be about a 7 point advantage to playing at home. The great thing about this system is that it is straightforward to predict the outcomes of future matches. For example, if Duke played Kansas at home for their next game, we would expect Duke to be 2 point underdogs (-3 – 6 + 6.9).

scorediff_ex <- ex_data %>%
  mutate(scorediff = WScore - LScore) %>%
  select(scorediff) %>% 
  as.matrix()

results_ex <- lsei(srs_ex, scorediff_ex)

data_frame(Team = colnames(srs_ex),
           SRS = results_ex[[1]])

## # A tibble: 5 x 2
##        Team       SRS
##       <chr>     <dbl>
## 1       loc  6.873239
## 2     UConn  4.204225
## 3    Kansas  6.021127
## 4      Duke -3.007042
## 5 Wisconsin -7.218310

Hopefully this post provided some background into some rating systems in College Basketball and how these metrics can be easily calculated for any team (and any sport) using a few custom functions in R. In my next post, I will use these systems to look at how tournament teams have performed in the past to get some insight into what might happen in 2018.