AFL results since 1998 from http://probabilistic-footy.monash.edu/
Load all results.
res <- read.delim('results.tsv')
#colnames(res)
# Add a score margin relative to home team
res <- res %>% mutate(margin=h_total-a_total)
Mean and std.dev of margin for home-team by year
res %>% group_by(year) %>% summarise(mean=mean(margin), std.dev=sd(margin)) %>% arrange(desc(year))
Check number of games by year
res %>% group_by(year) %>% summarise(n=n()) %>% arrange(desc(year))
Check home-team margin per year
res %>% ggplot(aes(factor(year),margin)) + geom_boxplot() + ggtitle("Margin by year")
Margins as a density for last 5 years
res %>% filter(year>2011) %>% ggplot(aes(color=factor(year),margin)) + geom_density() + ggtitle("Margin by year")
Check for normality of margins.
shapiro.test(res$margin)
Shapiro-Wilk normality test
data: res$margin
W = 0.99951, p-value = 0.509
qq-plot to visualise normality
qqnorm(res$margin)
qqline(res$margin, col=2)
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKQUZMIHJlc3VsdHMgc2luY2UgMTk5OCBmcm9tIFtodHRwOi8vcHJvYmFiaWxpc3RpYy1mb290eS5tb25hc2guZWR1L10oaHR0cDovL3Byb2JhYmlsaXN0aWMtZm9vdHkubW9uYXNoLmVkdS8pCgoKYGBge3IgaW5pdCxlY2hvPUZBTFNFLHJlc3VsdHM9J2hpZGUnLG1lc3NhZ2U9RkFMU0V9CmxpYnJhcnkoZHBseXIpCmxpYnJhcnkodGlkeXIpCmxpYnJhcnkoZ2dwbG90MikKYGBgCkxvYWQgYWxsIHJlc3VsdHMuCmBgYHtyfQpyZXMgPC0gcmVhZC5kZWxpbSgncmVzdWx0cy50c3YnKQojY29sbmFtZXMocmVzKQojIEFkZCBhIHNjb3JlIG1hcmdpbiByZWxhdGl2ZSB0byBob21lIHRlYW0KcmVzIDwtIHJlcyAlPiUgbXV0YXRlKG1hcmdpbj1oX3RvdGFsLWFfdG90YWwpCmBgYAoKIyMgTWVhbiBhbmQgc3RkLmRldiBvZiBtYXJnaW4gZm9yIGhvbWUtdGVhbSBieSB5ZWFyCmBgYHtyfQpyZXMgJT4lIGdyb3VwX2J5KHllYXIpICU+JSBzdW1tYXJpc2UobWVhbj1tZWFuKG1hcmdpbiksIHN0ZC5kZXY9c2QobWFyZ2luKSkgJT4lIGFycmFuZ2UoZGVzYyh5ZWFyKSkKYGBgCgpDaGVjayBudW1iZXIgb2YgZ2FtZXMgYnkgeWVhcgpgYGB7cn0KcmVzICU+JSBncm91cF9ieSh5ZWFyKSAlPiUgc3VtbWFyaXNlKG49bigpKSAlPiUgIGFycmFuZ2UoZGVzYyh5ZWFyKSkKYGBgCgojIyBDaGVjayBob21lLXRlYW0gbWFyZ2luIHBlciB5ZWFyCmBgYHtyfQpyZXMgJT4lIGdncGxvdChhZXMoZmFjdG9yKHllYXIpLG1hcmdpbikpICsgZ2VvbV9ib3hwbG90KCkgKyBnZ3RpdGxlKCJNYXJnaW4gYnkgeWVhciIpCmBgYAoKIyMgTWFyZ2lucyBhcyBhIGRlbnNpdHkgZm9yIGxhc3QgNSB5ZWFycwpgYGB7cn0KcmVzICU+JSBmaWx0ZXIoeWVhcj4yMDExKSAlPiUgZ2dwbG90KGFlcyhjb2xvcj1mYWN0b3IoeWVhciksbWFyZ2luKSkgKyBnZW9tX2RlbnNpdHkoKSArIGdndGl0bGUoIk1hcmdpbiBieSB5ZWFyIikKYGBgCgojIyBDaGVjayBmb3Igbm9ybWFsaXR5IG9mIG1hcmdpbnMuCmBgYHtyfQpzaGFwaXJvLnRlc3QocmVzJG1hcmdpbikKYGBgCgpxcS1wbG90IHRvIHZpc3VhbGlzZSBub3JtYWxpdHkKYGBge3J9CnFxbm9ybShyZXMkbWFyZ2luKQpxcWxpbmUocmVzJG1hcmdpbiwgY29sPTIpCmBgYAoKCg==