# Advent Calendar of Football Trivia Analyses

**rstats on Robert Hickman**, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

One of the most consistent fonts of posts on this blog is The Guardian’s football trivia page The Knowledge. A particular reason for this is that the small contained questions lend themselves to small blogposts that I can turn around in an hour or two, as opposed to being endlessly redrafted until I lose interest.

However, I still sometimes don’t quite get round to finishing some of these posts, or have trouble justifying a blog post on a very small and ‘trivial’ answer to a question. Therefore, as a sort of end-of-year round up, and a Christmas present to myself, I wanted to push out answers to questions I found particularly interesting over the last year and hadn’t quite got round to ^{1}. I’ll probably add them all to this post as I finish them up.

# 2nd December – Everything in its right place

I wonder if any of any sporting leagues have ever ended in alphabetical order? pic.twitter.com/you6u8Uzwz

— P A Hunt (@TeachFMaths) June 15, 2019

## Answer – yes, kind of. But also no.

This question has actually been answered (as many of these will have been). For a league of 20 teams (like the English Premier League), we might imagine if would have happened over the last ~150 years, but it’s almost certain from some basic maths that it won’t have, and moreover, will never happen.

Let’s load some data and see why.

#as per usual, going to heavily rely on tidyverse #and engsoccerdata throughout these posts library(tidyverse) library(engsoccerdata) #load English league data league_data <- engsoccerdata::england %>% #select and gather match results select(season = Season, division, home, visitor, hgoal, vgoal) %>% gather("location", "team", -season, -division, -hgoal, -vgoal) %>% mutate( g_for = case_when( location == "home" ~ hgoal, location == "visitor" ~ vgoal ), g_ag = case_when( location == "home" ~ vgoal, location == "visitor" ~ hgoal )) %>% #get correct point for a win/loss mutate( points = case_when( g_for > g_ag & season < 1981 ~ 2, g_for > g_ag & season > 1980 ~ 3, g_for == g_ag ~ 1, g_for < g_ag ~ 0 ), gd = g_for - g_ag ) %>% #group by season and league and get final tables group_by(season, division, team) %>% summarise(points = sum(points), gd = sum(gd), g_for = sum(g_for)) %>% arrange(-points, -gd, -g_for) %>% #rank league order and alphabetical order mutate(league_pos = rank(-points, ties.method = "first"), alph_order = rank(team, ties.method = "first")) %>% select(season, division, team, league_pos, alph_order) %>% #split by league and season split(., f = list(.$season, .$division)) %>% keep(function(x) nrow(x) > 0) #print the top of the first league table head(league_data[[1]]) ## # A tibble: 6 x 5 ## # Groups: season, division [1] ## season division team league_pos alph_order #### 1 1888 1 Preston North End 1 9 ## 2 1888 1 Aston Villa 2 2 ## 3 1888 1 Wolverhampton Wanderers 3 12 ## 4 1888 1 Blackburn Rovers 4 3 ## 5 1888 1 Bolton Wanderers 5 4 ## 6 1888 1 West Bromwich Albion 6 11

We can then run a load of Spearman’s rank correlation tests on the data to see which ones are perfectly correlated or anti-correlated in both league and alphabetical ranking. We’ll use the very handy broom package to tidy the results of our many tests into one data.frame (remove the filter at the end of the pipe chain to see what gets output).

#use broom to tidily do stats library(broom) #correlate league and alphabetical order by year exact_correlations <- league_data %>% map_df(., function(data) { cor.test( data$league_pos, data$alph_order, method = "spearman" ) %>% tidy() %>% mutate(season = unique(data$season), division = unique(data$division)) }) %>% #take only significantly filter(abs(statistic) == 1) print(exact_correlations) ## # A tibble: 0 x 7 ## # ... with 7 variables: estimate, statistic , p.value , ## # method , alternative , season , division

And so we find no exact correlations. There are no instances in 363 separate seasons of English league football where teams line up in either alphabetical, or anti-alphabetical order.

Let’s see why this is. To make things simpler, I’m going to imagine a cutdown league of only 6 teams using teams starting with each of the first 6 letter of the alphabet:

first_letter_names <- league_data %>% bind_rows() %>% ungroup() %>% #get first letter of a team name mutate(first_letter = gsub("(^.)(.*)", "\\1", team)) %>% filter(season > 1992 & division == 1 & first_letter %in% toupper(letters[1:6]) ) %>% #get one team beginning with A, B, C... filter(!duplicated(first_letter)) %>% select(team) %>% arrange(team) %>% print() ## # A tibble: 6 x 1 ## team #### 1 Arsenal ## 2 Blackburn Rovers ## 3 Coventry City ## 4 Derby County ## 5 Everton ## 6 Fulham

For the league to finish in alphabetical order, we first need the team that is first alphabetically (Arsenal) to finish in first position. Assuming all teams have an equal chance of winning the league, the chance of this is obviously

\[ p(Arsenal = 1) = \frac{1}{n}\]

Then we need the second team (Blackburn Rovers), to finish in second. This is predicated on Arsenal already finishing in first position, so the chance becomes

\[ p(Blackburn = 2 | Arsenal = 1) = \frac{1}{n-1} \]

and so on until the last team (Fulham) just have to slot into the only position left (n, 6th in our example)

Thus the total chance becomes

\[ \frac{1}{n} \cdot \frac{1}{n-1} … \cdot \frac{1}{1} \]

which can also be written

\[ p(ordered) = \prod_{n = 1}^{N} \frac{1}{n}\]

which multiplies out to

\[ p(ordered) = \frac{1}{n!} \]

so for our very small league the chance of n (assumed equally strong teams)

factorial(nrow(first_letter_names)) ## [1] 720

so we have a 1/720 chance that this league ends perfectly in alphabetical order. For bigger leagues (for reference most large European leagues contain 18-24 teams) this denominator grows *super-exponentially* and becomes tiny.

For the English Premier League (20 teams) for instance the chance becomes

league_data %>% bind_rows() %>% ungroup() %>% filter(season == max(season) & division == 1) %>% nrow() %>% factorial() ## [1] 2.432902e+18

or 1 in 2.4 quintillion. In short, if it’s assumed that there’s no relation between order of names and team strength, we might expect the universe to end before all 20 teams finish in perfect order.

We can test if our predictions bear out by looking at tiny leagues with small numbers of teams, e.g. the group stages of the Champions/Europa Leagues.

First we need to scrape the final tables for the last 8 years of data from both competitions:

library(rvest) #website to scrape group stage data from fb_data <- "https://footballdatabase.com" ucl_links <- sprintf( "/league-scores-tables/uefa-champions-league-20%s-%s", 10:18, 11:19 ) europa_links <- sprintf( "/league-scores-tables/uefa-europa-league-20%s-%s", 10:18, 11:19 ) #function to scrape the data from these links get_competition_data <- function(competition, links) { data <- links %>% paste0(fb_data, .) %>% map_df(., function(year) { page_read <- read_html(year) groups <- letters[1:8] %>% map_df(., function(group) { page_read %>% html_nodes(sprintf("#total-group-%s > div > table", group)) %>% html_table(fill = TRUE) %>% as.data.frame() %>% mutate(group) }) %>% mutate(year = gsub("(.*-)([0-9]{4}-[0-9]{2})", "\\2", year)) }) %>% mutate(competition) } #scrape and bind the data uefa_data <- bind_rows( get_competition_data("champions", ucl_links), get_competition_data("europa", europa_links) ) #print a cutdown version of the scraped data head(uefa_data %>% select(club = Club, points = P, year, competition)) ## club points year competition ## 1 Tottenham Hotspur 11 2010-11 champions ## 2 Inter Milan 10 2010-11 champions ## 3 FC Twente 6 2010-11 champions ## 4 Werder Bremen 5 2010-11 champions ## 5 Schalke 04 13 2010-11 champions ## 6 Lyon 10 2010-11 champions

So now we have 128 (8 groups x 8 years x 2 competitions) ‘mini-leagues’ each of 4 teams.

We can then munge this data to find all the groups where the teams finish in alphabetical order. We’d expect 128/4! leagues to finish in alphabetical order (or 5.33 to be exact).

ordered_groups <- uefa_data %>% #select relevant informatiob select(team = Club, league_pos = X., group, year, competition) %>% #by group find where teams finish in alphabetical order group_by(year, group, competition) %>% mutate(alph_order = rank(team, ties.method = "first")) %>% filter(league_pos == alph_order) %>% #keep only group where all (4) teams finish in order summarise(n = n()) %>% filter(n == 4) %>% #join and filter back data left_join(uefa_data, ., by = c("group", "year", "competition")) %>% filter(!is.na(n)) %>% #select useful information select(team = Club, points = P, gd = X..., league_pos = X., group, year, competition) %>% #split groups up split(., list(.$year, .$group, .$competition)) %>% keep(function(x) nrow(x) > 0)

which leaves us with 5 leagues that have finished in order! almost exactly what we’d predict by chance if the first letter of a teams name had no effect on the outcome.

ordered_groups ## $`2011-12.c.champions` ## team points gd league_pos group year competition ## 5 Benfica 12 4 1 c 2011-12 champions ## 6 FC Basel 11 1 2 c 2011-12 champions ## 7 Manchester United 9 3 3 c 2011-12 champions ## 8 Otelul Galati 0 -8 4 c 2011-12 champions ## ## $`2015-16.c.champions` ## team points gd league_pos group year competition ## 9 Atlético Madrid 13 8 1 c 2015-16 champions ## 10 Benfica 10 2 2 c 2015-16 champions ## 11 Galatasaray 5 -4 3 c 2015-16 champions ## 12 Lokomotiv Astana 4 -6 4 c 2015-16 champions ## ## $`2010-11.f.champions` ## team points gd league_pos group year competition ## 1 Chelsea FC 15 10 1 f 2010-11 champions ## 2 Marseille 12 9 2 f 2010-11 champions ## 3 Spartak Moskva 9 -3 3 f 2010-11 champions ## 4 Žilina 0 -16 4 f 2010-11 champions ## ## $`2015-16.g.champions` ## team points gd league_pos group year competition ## 13 Chelsea FC 13 10 1 g 2015-16 champions ## 14 Dynamo Kyiv 11 4 2 g 2015-16 champions ## 15 FC Porto 10 1 3 g 2015-16 champions ## 16 Maccabi Tel Aviv FC 0 -15 4 g 2015-16 champions ## ## $`2018-19.h.champions` ## team points gd league_pos group year competition ## 17 Juventus 12 5 1 h 2018-19 champions ## 18 Manchester United 10 3 2 h 2018-19 champions ## 19 Valencia 8 0 3 h 2018-19 champions ## 20 Young Boys 4 -8 4 h 2018-19 champions

We can also do a larger test by randomly selecting teams out of the English league data we looked at earlier. To do this I need two quick functions: one to sample randomly from the data, and another to carry out the correlation test.

The first takes a number of samples (how many tests to run) and then selects a number of teams from each league sample. For instance, if I chose 3 teams, it might select Liverpool, Manchester United, and Watford, from the last season of the Premier League. These teams finished 2nd, 6th, and 11th respectively, so this ‘sampled league’ would fulfill the criteria of finishing in alphabetical order.

set.seed(3459) #take a random sample of leagues and teams withing those leagues sample_cutdown_leagues <- function(nteams, nsamples, data) { samples <- sample(length(data), nsamples, replace = TRUE) sampled_league_data <- data[samples] league_team_serials <- sampled_league_data %>% lapply(., nrow) %>% lapply(., sample, size = nteams) #carry out the correlation test league_cor_test <- map2_df( .x = sampled_league_data, .y = league_team_serials, .f = cor_test_data ) } #function for correlation test cor_test_data <- function(full_league_data, sampled_teams) { sampled_league <- full_league_data[sampled_teams,] %>% arrange(league_pos) cor_test <- cor.test( sampled_league$league_pos, sampled_league$alph_order, method = "spearman" ) %>% tidy() %>% #mutate on information about that season and teams chosen mutate(teams = paste(sampled_league$team, collapse = ", "), season = unique(sampled_league$season), division = unique(sampled_league$division)) }

So for instance if I just run it once, randomly selecting 4 teams:

test <- sample_cutdown_leagues(4, 1, league_data) #print the teams selected test$teams ## [1] "Brentford, Bristol Rovers, Brighton & Hove Albion, Chester" test ## # A tibble: 1 x 8 ## estimate statistic p.value method alternative teams season division #### 1 0.8 2.00 0.333 Spearma~ two.sided Brentfor~ 1994 3

It gives me 4 teams from the 1994 division 3 who didn’t finish in alphabetical order (though, amusingly, all have a very similar starting letter).

We can then carry this out with 10000 samples for n_team numbers of 2:6 to see if we get roughly the expected numbers of exactly correlated league finish positions (this will take 1-2mins) by finding out how many tests give an estimate of 1 (finished exactly correlated with alphabetical order) or -1 (finished exactly anti-correlated with alphabetical order).

Both these numbers should be roughly equal to the number of samples (10000) divided by the factorial of the number of teams selected.

test_n_numbers <- function(nteams) { #run sampling function n times #10k should do sampling <- sample_cutdown_leagues(nteams, 10000, league_data) #find exactly correlated and anti-correlated examples #where teams are in exact alphabetical order ascending or descending correlated <- length(which(sampling$estimate == max(sampling$estimate))) anti_correlated <- length(which(sampling$estimate == min(sampling$estimate))) expected <- nrow(sampling) / factorial(nteams) df <- data.frame(n = nteams, sample_cor = correlated, sample_anticor = anti_correlated, sample_expected = expected) } #run the function testing <- map_df(2:6, test_n_numbers) #print results print(testing) ## n sample_cor sample_anticor sample_expected ## 1 2 5010 4990 5000.00000 ## 2 3 1676 1665 1666.66667 ## 3 4 367 398 416.66667 ## 4 5 101 81 83.33333 ## 5 6 14 15 13.88889

And the numbers line up, as we would expect if there is no effect of the first letter of a team’s name upon final league position.

Finally, we can do a Kendall’s correlation test to really see if there is any relationship between alphabetical team name order and final league finish for all out our English league data. We use Kendall instead of a Spearman test here because we grouping all the data together we’re going to have a lot of ties (one team has to finish 1st in every league each year).

all_data <- league_data %>% bind_rows() #do a big correlation test kendall_test <- cor.test(all_data$alph_order, all_data$league_pos, alternative = "two.sided", method = "kendall") %>% tidy() %>% print() ## # A tibble: 1 x 5 ## estimate statistic p.value method alternative #### 1 0.0135 1.74 0.0826 Kendall's rank correlation tau two.sided

And we can see that, even though our p-value is ‘approaching significance’, it’s not significant at our fairly liberal threshold of 0.05. Even then, the effect size (0.013) is tiny, so there’s no need for Watford to start worrying just yet.

- SMALL DIGRESSION: I love blogging on this site and it also has been a great help to me in numerous ways (practice coding/writing, feeling like a “programmer”, for job interviews), but quite a lot of the time feel posts are not quite where I want them (I’m sure this feeling isn’t restricted to me) and so won’t put them up and so that time (sometimes quite a few hours!) I put into them in my spare time feels wasted and makes me feel worse about myself. I’m hoping that pushing out fairly rushed/half formed ideas like this will help with this.

**leave a comment**for the author, please follow the link and comment on their blog:

**rstats on Robert Hickman**.

R-bloggers.com offers

**daily e-mail updates**about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.