Skip to main content

Win-stay lose-shift strategy in formation changes in football

Abstract

Managerial decision making is likely to be a dominant determinant of performance of teams in team sports. Here we use Japanese and German football data to investigate correlates between temporal patterns of formation changes across matches and match results. We found that individual teams and managers both showed win-stay lose-shift behavior, a type of reinforcement learning. In other words, they tended to stick to the current formation after a win and switch to a different formation after a loss. In addition, formation changes did not statistically improve the results of succeeding matches. The results indicate that a swift implementation of a new formation in the win-stay lose-shift manner may not be a successful managerial rule of thumb.

1 Introduction

Exploring rules governing decision making has been fascinating various fields of research, and its domain of implication ranges from our daily lives to corporate and governmental scenes. In economic contexts in a widest sense, individuals often modify their behavior based on their past experiences, attempting to enhance the benefit received in the future. Such decision making strategies are generally called reinforcement learning. In reinforcement learning, behavior that has led to a large reward will be selected with a larger frequency, or the behavior will be incrementally modified toward the rewarded one. Reinforcement learning is common in humans [1, 2] and non-humans [3], is implemented with various algorithms [4], has theoretical underpinnings [1, 4], and has neural substrates [5, 6].

A simple version of reinforcement learning is the so-called win-stay lose-shift (WSLS) strategy [7, 8]. An agent adopting this strategy sticks to the current behavior if the agent is satisfied. The agent changes its behavior if unsatisfied. Experimental studies employing human participants have provided a line of evidence in favor of WSLS in situations such as repeated Prisoner’s Dilemma [9, 10], gambling tasks [11, 12], and tasks in which participants construct virtual stone tools [1315]. It has also been suggested in nonscientific contexts that decisions by athletes and gamblers are often consistent with WSLS patterns even if the outcome of games seems to be independent of the decision [16].

Association football (also known as soccer; hereafter refer to it as football) is one of the most popular sports in the world and provides huge business opportunities. The television rights of the English Premier League yield over two billion euros per year [17]. Transfer fees of top players can be tens of millions of euros [18]. Various aspects of football, not only watching but also betting [19] and the history of tactics [20], enjoy popularity. Football and other team sports also provide data for leadership studies because a large amount of sports data is available and the performance of teams and players can be unambiguously measured by match results [2123].

In the present study, using data obtained from football matches, we examine the possibility that managers of teams use the WSLS strategy. Managers can affect the performance of teams through selections of players, training of players, and implementation of tactics including formations [18]. In particular, a formation is a part of tactics to determine how players participate in offense and defense [24] and considered to affect match results [24, 25]. Managerial decision making in substituting players during a match may affect the probability of winning [25]. We hypothesize that a manager continues to use the same formation if he has won the previous match, whereas he experiments on another formation following a loss in the previous match.

The WSLS and more general reinforcement learning posit that unsuccessful individuals modify their behavior to increase the probability of winning. Therefore, we are interested in whether a formation change improves the performance of a team. To clarify this point, we also investigate effects of formation changes on the results of succeeding matches.

2 Materials and methods

2.1 Data set

We collected data on football matches from two websites, J-League Data Site (officially, “J. League”) [26] on J-League, and Kicker-online [27] on Bundesliga. J-League and Bundesliga are the most prestigious professional football leagues in Japan and Germany, respectively. We refer to the two data sets as the J-League and the Bundesliga data sets. The two data sets contain, for each team and match, the season, date, manager’s name, result (i.e., win, draw, or loss), and starting formation. Basic statistics of the data sets are summarized in Table 1. The distributions of the probability of winning for teams and managers are shown in Figure 1 for the two data sets. Because the strength of a team apart from the manager was considered to affect the probability of winning, in Figure 1, we treated a manager as different data points when he directed different teams. The same caveat applies to all the following analysis focusing on individual managers (Figures 3-6).

Figure 1
figure 1

Distributions of the probability of winning. (a) Distribution for the teams in the J-League data. (b) Distribution for the managers in the J-League data. (c) Distribution for the teams in the Bundesliga data. (d) Distribution for the managers in the Bundesliga data. The colored bars correspond to the teams or managers that have played at least 100 matches.

Table 1 Statistics of the J-League and Bundesliga data sets

Between 1993 and 2004, except for 1996, each season of J-League was divided into two half seasons. After the two half seasons had been completed, two champion teams, each representing a half season, played play-off matches. We regarded each half season as a season because intervals between two half seasons ranged from ten days to two months and therefore are longer than one week, which was a typical interval between two matches within a season. We also carried out the same analysis when we regarded one year, not one half season, as a season and verified that the main results were unaltered (Appendix A).

We also collected data on Bundesliga from another website, Fussballdaten [28]. We focused on the Kicker-online data rather than the Fullballdaten data because the definition of the position was coarser for the Fussballdaten data (i.e., a player was not assumed to change his position during a season) than the Kicker-online data. Nevertheless, to verify the robustness of the following results, we also analyzed the Fussballdaten data (Appendix B).

2.2 Definition of formation

The definition of formation was different between the two data sets. In the J-League data, each of the ten field players was assigned to either defender (DF), midfielder (MF), or forward (FW) in each match. We defined formation as a triplet of the numbers of DF, MF, and FW players, which sum up to ten. For example, a formation 4-4-2 implies four DFs, four MFs, and two FWs.

In the Bundesliga data, the starting positions of the players were given on a two-dimensional map of the pitch (Figure 2). For this data set, we defined formation as follows. First, we measured the distance between the goal line and the bottom edge of the image representing each player along the vertical axis (e.g., 113 shown in Figure 2). We referred to the HTML source code of Kicker-online to do this. The unit of the distance is pixel (px). The distance between the goal line and the half-way line is between 45 and 60 m in real fields. The same distance is approximately equal to 500 px in Kicker-online. Therefore, 1 px in Kicker-online roughly corresponds to 10 cm in real fields. Although the HTML source code also included the distance of players from the left touch line, we neglected this information because the primary determinant of the player’s position seems to be the distance from the goal line rather than that from the left or right touch line, as implied by the terms DF, MF, and FW. Second, we grouped players whose distances from the goal line were the same. Third, we ordered the groups of players in terms of the distance, resulting in an ordered set of the numbers of players at each distance value. The set of numbers defined a formation. For example, when the distances of the ten field players are equal to 113 px, 113 px, 113 px, 113 px, 236 px, 236 px, 359 px, 359 px, 359 px, and 441 px, the formation of the team is defined to be 4-2-3-1 (Figure 2).

Figure 2
figure 2

Definition of formation in the Bundesliga data. Kicker-online gives the starting positions of the players as two-dimensional coordinates on the pitch. Field players with the identical distance from the goal line are aggregated into the same position. The starting positions shown in the figure are coded as 4-2-3-1.

Among all matches in the Bundesliga data, the smallest nonzero distance between two players was equal to 31 px. Therefore, we did not have to worry about the possibility that players possessed almost the same distance values while being classified into distinct positions. For example, there was no case in which the distances of two field players from the goal line were equal to 113 px and 114 px.

For both data sets, a formation was defined as an ordered set of numbers, whereas the definition differs for the two data sets. For example, forward players possessing distance values 359 and 441 were classified into different positions in the Bundesliga data, whereas they belonged to the same position in the J-League data if they were both assigned to FW. In the following, we regarded that formation was changed when the ordered set of numbers differed between two consecutive matches.

Figures 3(a)-(d) show the distribution of the probability that a team or manager has changed the formation in the two leagues. To calculate the probability of formation changes for a team, we excluded the first match in each season and the matches immediately after a change in the manager. As in the case of formation changes, we regarded that a manager was changed when the manager directed a team in a given match but did not do so in the next match. With this definition, a short absence of a chief manager due to illness, for example, may induce formation changes. However, we adhered to this definition because of the lack of further information on behavior of managers. In addition, as explained in Section 2.1, we treated a manager as different data points when he led different teams.

Figure 3
figure 3

Distributions of the probability of formation changes. (a) Distribution for the teams in the J-League data. (b) Distribution for the managers in the J-League data. (c) Distribution for the teams in the Bundesliga data. (d) Distribution for the managers in the Bundesliga data. In (a)-(d), the colored bars correspond to the teams or managers that have played at least 100 matches. (e) Probability of formation change in each season in J-League, aggregated over the different teams and managers. A circle represents a season. Because the J-League data set did not have the information on managers between 1993 and 1998, we neglected changes of managers between 1993 and 1998. Between 1993 and 2004, except for 1996, intervals between two circles are dense because a season consists of two half seasons. (f) Probability of formation change in each season in Bundesliga, aggregated over the different teams and managers.

The frequency of formation changes as a function of time is shown in Figures 3(e) and 3(f) for J-League and Bundesliga, respectively. The figures suggest that the frequency of formation change is stable over years in J-League, but not in Bundesliga. Finally, we measured burstiness and memory coefficient [29] for interevent times of formation changes to quantify temporal patterns of formation changes. The results are shown in Appendix C.

2.3 GLMM

To statistically examine whether patterns of formation changes were consistent with WSLS behavior, we investigated effects of previous matches and other factors on the likelihood of formation change for each team. If managers used the WSLS, the effect of the win and loss in the previous match on the likelihood of formation change should be significantly negative and positive, respectively. We used a generalized linear mixed model (GLMM) with binomial errors and a logit-link function.

The dependent variable was the occurrence or lack thereof of formation changes, which was binary. As independent variables, we included the binary variable representing whether or not the stadium was the home of the team (i.e., home or away) and the ternary result of the previous match (i.e., win, draw, or loss). We designated the draw as the reference category for the match result. Because the likelihood of formation changes may be affected by a streak of wins or losses, we also included the result of the second last match as an independent variable. The difference between the focal team’s strength and the opponent’s strength was also an independent variable. The strength of a team was defined by the probability of winning in the season. We estimated the strength of a team separately for each season because it can vary across seasons. The name of the manager was included as a random effect (random intercept).

In this and the following analysis, we excluded the first match in each season for each team because we considered that the result of the last match in the preceding season would not directly affect the first match in a new season. In addition, we excluded matches immediately after a change of manager because we were not interested in formation changes induced by a change of manager. We further excluded the second match in each season for each team from the GLMM analysis when we employed the result of the second last match as an independent variable. Because the J-League data set did not have the information on managers between 1993 and 1998, we only used data between 1999 and 2014 in the GLMM analysis. We performed the statistical analysis using R 3.1.2 [30] with lme4 package [31].

2.4 Ordered probit model

We also investigated the effects of formation changes on match results. We used the ordered probit model because a match result was ternary. Because the strength was considered to heavily depend on teams, we controlled for the strength of teams. The same model was used for fitting match results in football in the Netherlands [32] and the UK [18].

The dependent variable of the model was a match result. We assumed that the occurrence of formation change (change or no change), the stadium (home or away), the strength of teams, and the result of the previous match (win, draw, or loss) can affect a match result. As a linear combination of these factors, we defined the unobserved potential variable for team i, denoted by \(\alpha_{i}\), by

$$ \alpha_{i} = \beta_{\mathrm{f}}f_{i} + \beta_{\mathrm{h}}h_{i}+\beta_{\mathrm{w}}w_{i} + \beta_{\ell}\ell_{i} +\beta_{\mathrm{ r}}r_{i}, $$
(1)

where \(f_{i}=1\) if team i changed the formation, and \(f_{i}=0\) otherwise; \(h_{i}=1\) if the stadium was the home of team i, and \(h_{i}=0\) otherwise; \(w_{i} = 1\) if team i won the previous match, and \(w_{i}=0\) otherwise; \(\ell_{i} = 1\) if team i lost the previous match, and \(\ell_{i}=0\) otherwise; the strength of team i denoted by \(r_{i}\) was defined as the fraction of matches that team i won in the given season. In Appendix D, we conducted the analysis by assuming that \(r_{i}\) was a latent variable obeying the normal distribution and then using the hierarchical Bayesian model [33].

Consider a match between home team i and away team j. We assumed that the match result, denoted by \(k_{ij}\), was determined by the difference between the potential values of the two teams, i.e.,

$$ y_{ij} \equiv\alpha_{i} - \alpha_{j}. $$
(2)

Variables \(y_{ij}\) and \(k_{ij}\) are related by

$$ k_{ij} = \textstyle\begin{cases} 2 \mbox{ (home team wins)} & \mbox{if } c_{1}< y_{ij} + \epsilon _{ij}, \\ 1 \mbox{ (draw)} & \mbox{if } c_{0}< y_{ij} + \epsilon_{ij} < c_{1}, \\ 0 \mbox{ (home team loses)} & \mbox{if } y_{ij} + \epsilon_{ij}< c_{0}, \end{cases} $$
(3)

where \(c_{0}\) and \(c_{1}\) are threshold parameters, and \(\epsilon_{ij}\) is an error term that obeys the normal distribution with mean 0 and standard deviation 1. Because \(h_{i}-h_{j}=1\), \(\beta_{\mathrm{h}}\) appears as a constant term on the right-hand side of Eq. (2). In fact, it is impossible to estimate \(\beta_{\mathrm{h}}\) because \(\beta _{\mathrm{h}}\) effectively shifts \(c_{0}\) and \(c_{1}\) by the same amount such that there are only two degrees of freedom in the parameter space spanned by \(c_{0}\), \(c_{1}\), and \(\beta_{\mathrm{h}}\). Therefore, we assumed \(c_{0}=-c_{1}\) and estimated \(c_{0}\) and \(\beta _{\mathrm{h}}\). This assumption did not alter the estimates of the other parameters. Equation (3) results in

$$\begin{aligned}& P(k_{ij}=2) = 1-\Phi(c_{1}-y_{ij}), \end{aligned}$$
(4)
$$\begin{aligned}& P(k_{ij}=1) = \Phi(c_{1}-y_{ij}) - \Phi(c_{0}-y_{ij}), \end{aligned}$$
(5)
$$\begin{aligned}& P(k_{ij}=0) = \Phi(c_{0}-y_{ij}), \end{aligned}$$
(6)

where P denotes the probability, and \(\Phi(\cdot)\) is the cumulative standard normal distribution function.

We excluded the matches that were the first game in a season at least for either team. We also excluded matches immediately after a change of manager in either team. Because the J-League data set did not contain the information on managers between 1993 and 1998, we only used data between 1999 and 2014 in this analysis. We performed the analysis using R 3.1.2 [30] and maxLik package [34].

2.5 Influence of individual manager’s behavior on match results

Different managers may show WSLS behavior to different extents to respectively affect match results. Therefore, we analyzed data separately for individual managers. For each manager i, we calculated the probability of winning under each of the following four conditions: (i) i’s team won the previous match, and i changed the formation, (ii) i’s team won the previous match, and i did not change the formation, (iii) i’s team lost the previous match, and i changed the formation, and (iv) i’s team lost the previous match, and i did not change the formation. We then compared the probability of winning between cases (i) and (ii), and between cases (iii) and (iv) using the paired t-test. In the t-test, we included the managers who directed at least ten pairs of consecutive matches in both of the two cases in comparison. In this and the next sections, we treated a manager as different data points when he directed different teams, as explained in Section 2.1. In addition, we excluded the pairs of consecutive matches when the managers changed the team between the two matches.

2.6 Degree of win-stay lose-shift

To further examine possible relationships between manager’s behavior and match results, we looked at the relationships between the tendency of the WSLS behavior for each manager (degree of WSLS for short) and the probability of winning. The degree of WSLS is defined by

$$\begin{aligned}& \mbox{degree of WSLS} \\& \quad= \bigl|P(\mbox{change}|\mbox{win}) - P_{\mathrm{WSLS}}(\mbox{change}|\mbox{win}) \bigr|+ \bigl|P(\mbox{change}|\mbox{loss})- P_{\mathrm{WSLS}}(\mbox{change}|\mbox{loss}) \bigr| \\& \quad= \bigl|P(\mbox{change}|\mbox{win}) - 0 \bigr|+ \bigl|P(\mbox{change}|\mbox{loss}) - 1 \bigr| \\& \quad= P(\mbox{change}|\mbox{win})+1-P(\mbox{change}|\mbox{loss}), \end{aligned}$$
(7)

where \(P_{\mathrm{WSLS}}(\mbox{change}|\mbox{win})\) (=0) is the conditional probability that a perfect WSLS manager changes the formation after winning, and likewise for \(P_{\mathrm{WSLS}}(\mbox{change}|\mbox{loss})\) (=1). The degree of WSLS ranges from 0 to 2.

3 Results

3.1 Win-stay lose-shift behavior in formation changes

We examined the extent to which managers possibly changed the formation of the team after losing a match and persist to the current formation after a win. The results of the GLMM analysis with the results of the previous matches being the only independent variables are shown in Table 2. For both data sets, winning in a match significantly decreased the probability of formation change in the next match, and losing in a match increased the probability of formation change. The results did not essentially change when we used the full set of independent variables (Table 3). Formation changes are consistent with WSLS patterns.

Table 2 Results of the GLMM analysis when the results of the previous match were used as the sole independent variables
Table 3 Results of the GLMM analysis when all the independent variables were considered

For the J-League data, the effects of all the additional independent variables were insignificant. We analyzed the J-League data by regarding a pair of half seasons (i.e., an yearly season) as a season to confirm that the results remained qualitatively the same except that winning in the second last match also significantly decreased the probability of formation change (Appendix A). We also confirmed that matches played in the further past affect the probability of formation change to progressively small extents (Appendix E).

For the Bundesliga data, winning in the second last match also significantly decreased the probability of formation change in the extended GLMM model (Table 3). These results are consistent with WSLS behavior. We also found for the Bundesliga data that stronger teams less frequently changed the formation and that a team would not change the formation to fight home games. We also investigated the Fussballdaten data for Bundesliga, in which the definition of formation was different, and confirmed that managers tended to use the WSLS strategy (Appendix B).

3.2 Determinants of match results

The results obtained from the ordered probit model are shown in Table 4. For both data sets, formation changes did not significantly affect a match result. The result remained qualitatively the same when each pair of half seasons was considered as a season in the J-League data (Appendix A), and when the strength of team was assumed to be a latent variable in the ordered probit model (Appendix D). However, when the Fussballdaten data were used, formation changes significantly decreased the probability of winning (Appendix B). Table 4 also tells us the following. Trivially, stronger teams were more likely to win in both data sets. The home advantage was significant in both data sets, consistent with previous literature [18, 35]. In Bundesliga, a win tended to yield a poor result in the next match. This is consistent with negative persistence effects reported in previous literature [18], i.e., the results of the current and previous matches tend to be the opposite.

Table 4 Effects of variables on match results as obtained from the ordered probit model

Figure 4(a) shows the probability of winning after individual managers changed or did not change the formation after a win in the J-League data. A large circle in Figure 4 represents a manager who presented both types of actions (i.e., formation change after winning and no formation change after winning) at least ten times. A small circle represents a manager who presented either type of action less than ten times. The formation change does not appear to affect the probability of winning. This is also apparently the case for the actions after a loss (Figure 4(b)) and the Bundesliga data (Figures 4(c) and 4(d)). The results also appear to be insensitive to the unconditional probability of winning, which roughly corresponds to the position along the diagonal in Figure 4. To be quantitative, we conducted the paired t-test on the managers who submitted the two types of actions at least ten times in each case (managers shown by the large circles in Figure 4). For the J-League data, there was no significant effect of formation change on the probability of winning both after winning (\(p = 0.441\), \(n=3\); corresponding to Figure 4(a)) and losing (\(p = 0.404\), \(n=4\); Figure 4(b)). For the Bundesliga data, formation changes after winning significantly decreased the probability of winning in the next match (\(p = 0.026\), \(n=46\); Figure 4(c)), whereas there was no significant effect after losing (\(p = 0.533\), \(n=42\); Figure 4(d)). These results suggest that formation changes did not at least increase the possibility of winning.

Figure 4
figure 4

Conditional probability of winning for individual managers after they changed or did not change the formation. A circle represents a manager associated with a team, who directed the team in at least 100 matches. A large circle represents a manager who showed both types of the actions (e.g., formation change after a win and no formation change after a win in (a) and (c)) at least ten times. A small circle represents a manager who presented either type of action less than ten times. (a) When the team won the previous match in the J-League data. (b) When the team lost the previous match in the J-League data. (c) When the team won the previous match in the Bundesliga data. (d) When the team lost the previous match in the Bundesliga data.

The analysis with the ordered probit model aggregated the data from all managers. Therefore, we examined the relationship between the degree of WSLS and the probability of winning for individual managers. The results are shown in Figure 5. A circle in Figure 5 represents a manager. We did not find a significant relationship between the usage of the WSLS and the probability of winning for both J-League (Pearson’s \(r = 0.213\), \(p = 0.411\), \(n=17\)) and Bundesliga (\(r = 0.058\), \(p = 0.668\), \(n=58\)) data.

Figure 5
figure 5

Relationship between the degree of WSLS and the probability of winning. (a) J-League. (b) Bundesliga. A circle represents a manager associated with a team, who directed the team in at least 100 matches.

4 Discussion

We have provided evidence that football managers tend to stick to the current formation until the team loses, consistent with the WSLS strategy previously shown in laboratory experiments with social dilemma games [9, 10] and gambling tasks [11, 12]. Formation changes did not significantly affect (at least did not improve) a match result in most cases. This result seems to be odd because managers change formation to lead the team to a success. Generally speaking, when the environment in which an agent is located is fixed or exogenously changing, reinforcement learning usually improves the performance of the agent [4]. However, computational studies have suggested that it is not always the case when agents employing reinforcement learning are competing with each other, because the competing agents try to supersede each other [3639]. The present finding that manager’s WSLS behavior does not improve team’s performance is consistent with these computational results. Empirical studies also suggest that humans obeying reinforcement learning does not improve the performance in complex environments. For example, players in the National Basketball Association were more likely to attempt 3 point shots after successful 3 point shots although their probability of success decreased for additional shots [40]. Also in nonscientific accounts, it has been suggested that humans engaged in sports and gambles often use the WSLS strategy even if outcome of games is determined merely at random [16]. We have provided quantitative evidence underlying these statements.

Many sports fans possess the hot hand belief in match results, i.e., belief that a win or good performance persists [41]. However, empirical evidence supports that streaks of wins and those of losses are less likely to occur than under the independence assumption [41]. By analyzing patterns of matches in the top division of football in England, Dobson and Goddard suggested the existence of negative persistence effects, i.e., a team with consecutive wins tended to perform poorly in the next match and vice versa [18]. Their results are consistent with the present results; we observed the negative persistence effects, i.e., anticorrelation between the results of the previous and present matches.

In the present study, we have neglected various factors that potentially affect the likelihood of formation change because our data sets did not contain the relevant information. For example, managers may change formations due to injuries, suspensions of players, and other strategic reasons including transfer of players. More detailed data will be able to provide further understanding of the relative importance of strategic versus accidental factors in formation changes.

An important limitation of the present study is that we have oversimplified the concept of formation. Effective formations dynamically change during a match owing to movements of players. Because of the availability of data and our interests in the manager’s long-term behavior rather than formation changes during a match [25], we used the formation data released in the beginning of the matches. Based on recent technological developments, formations can be extracted from tracking data on movement patterns of players [42, 43]. Investigations on manager’s decision making using such technologies warrant further research.

References

  1. Fudenberg D, Levine DK (1998) The theory of learning in games. MIT Press, Cambridge

    Google Scholar 

  2. Camerer C (2003) Behavioral game theory: experiments in strategic interaction. Princeton University Press, Princeton

    Google Scholar 

  3. Pearce JM (2013) Animal learning and cognition: an introduction. Psychology Press, Hove

    Google Scholar 

  4. Sutton RS, Barto AG (1998) Reinforcement learning: an introduction. MIT Press, Cambridge

    Google Scholar 

  5. Schultz W, Dayan P, Montague PR (1997) A neural substrate of prediction and reward. Science 275(5306):1593-1599

    Article  Google Scholar 

  6. Glimcher PW, Camerer C, Fehr E, Poldrack RA (2009) Neuroeconomics: decision making and the brain. Academic Press, New York

    Google Scholar 

  7. Kraines D, Kraines V (1989) Pavlov and the prisoner’s dilemma. Theory Decis 26(1):47-79

    Article  MathSciNet  Google Scholar 

  8. Nowak M, Sigmund K (1993) A strategy of win-stay, lose-shift that outperforms tit-for-tat in the Prisoner’s Dilemma Game. Nature 364(6432):56-58

    Article  Google Scholar 

  9. Wedekind C, Milinski M (1996) Human cooperation in the simultaneous and the alternating Prisoner’s Dilemma: Pavlov versus Generous Tit-for-Tat. Proc Natl Acad Sci USA 93(7):2686-2689

    Article  Google Scholar 

  10. Milinski M, Wedekind C (1998) Working memory constrains human cooperation in the Prisoner’s Dilemma. Proc Natl Acad Sci 95(23):13755-13758

    Article  Google Scholar 

  11. Hayden BY, Platt ML (2009) Gambling for Gatorade: risk-sensitive decision making for fluid rewards in humans. Anim Cogn 12(1):201-207

    Article  Google Scholar 

  12. Scheibehenne B, Wilke A, Todd PM (2011) Expectations of clumpy resources influence predictions of sequential events. Evol Hum Behav 32(5):326-333

    Article  Google Scholar 

  13. Mesoudi A, O’Brien MJ (2008) The cultural transmission of Great Basin projectile-point technology I: an experimental simulation. Am Antiq 73(1):3-28

    Google Scholar 

  14. Mesoudi A, O’Brien MJ (2008) The cultural transmission of Great Basin projectile-point technology II: an agent-based computer simulation. Am Antiq 73(4):627-644

    Google Scholar 

  15. Mesoudi A (2014) Experimental studies of modern human social and individual learning in an archaeological context: people behave adaptively, but within limits. In: Akazawa T, Ogihara N, Tanabe HC, Terashima H (eds) Dynamics of learning in Neanderthals and modern humans, vol 2. Springer, Heidelberg, pp 65-76

    Chapter  Google Scholar 

  16. Vyse SA (2013) Believing in magic: the psychology of superstition, updated edn. Oxford University Press, Oxford

    Google Scholar 

  17. Panja T (2013) Top soccer leagues get 25% rise in TV rights sales, report says. In: Bloomberg Business. http://www.bloomberg.com/news/articles/2013-11-11/top-soccer-leagues-get-25-rise-in-tv-rights-sales-report-says. Accessed 28 Feb 2015

  18. Dobson S, Goddard J (2001) The economics of football. Cambridge University Press, Cambridge

    Book  Google Scholar 

  19. Dixon MJ, Coles SG (1997) Modelling association football scores and inefficiencies in the football betting market. J R Stat Soc, Ser C, Appl Stat 46(2):265-280

    Article  Google Scholar 

  20. Wilson J (2013) Inverting the pyramid: the history of soccer tactics. Nation Books, New York

    Google Scholar 

  21. Audas R, Dobson S, Goddard J (1997) Team performance and managerial change in the English Football League. Econ Aff 17(6):30-36

    Article  Google Scholar 

  22. Dawson P, Dobson S, Gerrard B (2000) Estimating coaching efficiency in professional team sports: evidence from English association football. Scott J Polit Econ 47(4):399-421

    Article  Google Scholar 

  23. Audas R, Dobson S, Goddard J (2002) The impact of managerial change on team performance in professional sports. J Econ Bus 54(6):633-650

    Article  Google Scholar 

  24. Bangsbo J, Peitersen B (2000) Soccer systems & strategies. Human Kinetics, Champaign

    Google Scholar 

  25. Hirotsu N, Wright MB (2006) Modeling tactical changes of formation in association football as a zero-sum game. J Quant Anal Sports 2(2):4

    MathSciNet  Google Scholar 

  26. J-League Data site https://data.j-league.or.jp/SFTP01/. Accessed 4-6 June 2014

  27. Kicker-online http://www.kicker.de/. Accessed 14-16 Aug 2014

  28. Fussballdaten http://www.fussballdaten.de/. Accessed 28-31 July 2014

  29. Goh KI, Barabási AL (2008) Burstiness and memory in complex systems. EPL 81:48002

    Article  Google Scholar 

  30. R Core Team (2014) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. http://www.R-project.org/

    Google Scholar 

  31. Bates D, Maechler M, Bolker B, Walker S (2014) lme4: linear mixed-effects models using Eigen and S4. http://CRAN.R-project.org/package=lme4. R package version 1.1-7

  32. Koning RH (2000) Balance in competition in Dutch soccer. J R Stat Soc, Ser D, Statist 49:419-431

    Article  Google Scholar 

  33. Parent E, Rivot E (2012) Introduction to hierarchical Bayesian modeling for ecological data. CRC Press, Boca Raton

    Google Scholar 

  34. Henningsen A, Toomet O (2011) maxLik: a package for maximum likelihood estimation in R. Comput Stat 26(3):443-458

    Article  MathSciNet  MATH  Google Scholar 

  35. Albert J, Koning RH (2010) Statistical thinking in sports. CRC Press, Boca Raton

    Google Scholar 

  36. Taiji M, Ikegami T (1999) Dynamics of internal models in game players. Physica D 134(2):253-266

    Article  MathSciNet  MATH  Google Scholar 

  37. Macy MW, Flache A (2002) Learning dynamics in social dilemmas. Proc Natl Acad Sci USA 99(Suppl 3):7229-7236

    Article  Google Scholar 

  38. Masuda N, Ohtsuki H (2009) A theoretical analysis of temporal difference learning in the iterated prisoner’s dilemma game. Bull Math Biol 71(8):1818-1850

    Article  MathSciNet  MATH  Google Scholar 

  39. Masuda N, Nakamura M (2011) Numerical analysis of a reinforcement learning model with the dynamic aspiration level in the iterated prisoner’s dilemma. J Theor Biol 278(1):55-62

    Article  MathSciNet  MATH  Google Scholar 

  40. Neiman T, Loewenstein Y (2011) Reinforcement learning in professional basketball players. Nat Commun 2:569

    Article  Google Scholar 

  41. Bar-Eli M, Avugos S, Raab M (2006) Twenty years of “hot hand” research: review and critique. Psychol Sport Exerc 7(6):525-553

    Article  Google Scholar 

  42. Bialkowski A, Lucey P, Carr P, Yue Y, Matthews I (2014) “Win at home and draw away”: automatic formation analysis highlighting the differences in home and away team behaviors. In: MIT Sloan sports analytics conference (SSAC)

  43. Bialkowski A, Lucey P, Carr P, Yue Y, Sridharanand S et al. (2014) Large-scale analysis of soccer matches using spatiotemporal data. In: IEEE international conference on data mining (ICDM)

  44. Stan Development Team (2014) Rstan: the R interface to Stan, version 2.5.0. http://mc-stan.org/rstan.html

Download references

Acknowledgements

This work is supported by JST, CREST.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Naoki Masuda.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors designed the study. KT collected and analyzed the data. All authors wrote the manuscript.

Appendices

Appendix A: Analysis of the J-League data on the basis of yearly seasons

From 1993 to 2004, except for 1996, each season of J-League, spanning a year, was subdivided into two half seasons. In the main text, we regarded each half season as a season. To examine the robustness of our results, we carried out the same analysis when we regarded one entire season (i.e., one year), not one half season, as a season. The results were qualitatively the same as those shown in the main text (Tables 57) except that winning in the second last match significantly decreased the probability of formation change.

Table 5 Results of the GLMM analysis for the J-League data when a year was regarded as a season
Table 6 Results of the GLMM analysis with all independent variables for the J-League data when a year was regarded as a season
Table 7 Effects of variables on match results for the J-League data when a year was regarded as a season

Appendix B: Fussballdaten

We analyzed data on Bundesliga from another website, Fussballdaten [28]. In the Fussballdaten data, each field player was assigned to one of the three positions (i.e., DF, MF, or FW) registered for an entire season. We defined the formation by counting the number of each type of field player in the same manner as that for the J-League data set.

First, to examine possible existence of WSLS behavior by managers, we applied the GLMM analysis to the Fussballdaten data. The results shown in Tables 8 and 9 are largely consistent with those for the Kicker-online data (Tables 1 and 2). In particular, winning and losing in the previous match significantly decreased and increased the probability of formation change in the next match, respectively, consistent with WSLS behavior.

Table 8 Results of the GLMM analysis for the Fussballdaten data when the results of the previous match were used as the sole independent variables
Table 9 Results of the GLMM analysis for the Fussballdaten data when all the independent variables were considered

Second, we also investigated the effect of formation change and other factors on the match result using the ordered probit model. Table 10 indicates that formation changes have decreased the probability of winning. This result is not consistent with those for the two data sets shown in the main text. In addition, winning in the previous match decreased the probability of winning in the next match, indicating the presence of the negative persistence effect. This result is consistent with that for the Kicker-online data (Table 4).

Table 10 Effects of variables on match results for the Fussballdaten data

Appendix C: Burstiness and memory coefficient of interevent time series

To capture temporal properties of formation changes, in this section we calculated burstiness, B, and memory coefficient, M, [29] on the basis of the interevent time series \(\{ \tau_{i} \}\) defined as follows. We calculated B and M for each manager. As in the main text, we treated a manager as different data points when he directed different teams. In the main text, we used interevent time series for individual seasons without concatenating different seasons. In this section, however, we use \(\{\tau^{k}_{i}\}\) obtained by concatenating all seasons. For a given manager, we denote by \(t_{0}, t_{1}, \ldots, t_{N}\) (\(2\le t_{0} < t_{1} < \cdots< t_{N}\)) the times when the manager changes the formation. The number of formation change summed over all seasons is equal to \(N+1\). We counted the time in terms of the number of match rather than the day to exclude effects of variable intervals between consecutive matches in terms of the real time. If \(t_{2} = 5\), for example, the manager changed the formation to play the fifth match, and it was the third change for the manager since the first match in the data set.

Managers sometimes moved from one team to another or did not direct any team. Because formation changes occurring as a result of a manager’s move or after a long absence were not considered to be strategic, we discarded the corresponding intervals. It should be noted that we did not mix interevent time series for a manager leading different teams. Then, a time series \(\{t_{i}\}\) for each manager was partitioned into \(N_{k}\) segments by either a manager’s move or absence. We denoted by \(N_{k}+1\) and K the total number of formation changes in the kth segment and the total number of segments, respectively. It holds true that \(\sum_{k=1}^{K} (N_{k} + 1) = N+1\). We also denoted by \(t^{k}_{i}\) (\(0\le i\le N_{k}\)) the time of the ith formation change in the kth segment. The interevent time for formation changes was defined by \(\tau^{k}_{i} = t^{k}_{i}-t^{k}_{i-1}\) (\(1\le i\le N_{k}\)). It holds true that \(\sum_{k=1}^{K} \sum_{i=1}^{N_{k}}\tau^{k}_{i}=N+1-K\). In the following analysis, we used managers who directed at least 100 matches in a team and has \(N+1-K \ge10\).

The burstiness is defined by

$$ B=\frac{\sigma/m-1}{\sigma/m+1}=\frac{\sigma-m}{\sigma+m}, $$
(8)

where \(m=\sum_{k=1}^{K} \sum_{i=1}^{N_{k}} \tau^{k}_{i}/(N+1-K)\) and \(\sigma=\sqrt{ \sum_{k=1}^{K} \sum_{i=1}^{N_{k}} (\tau^{k}_{i} - m)^{2} /(N+1-K)}\) are the mean and standard deviation of interevent time, respectively. B ranges between −1 and 1. A large value of B indicates that a sequence of formation change events is bursty in the sense that the interevent time obeys a long-tailed distribution. The Poisson process yields the exponential distribution and hence yields \(B = 0\).

The memory coefficient quantifies the correlation between two consecutive interevent times and is defined by

$$ M=\frac{1}{N-K} \sum_{k=1}^{K} \sum _{i=1}^{N_{k}-1}\frac{(\tau _{i}^{k}-m_{1})(\tau_{i+1}^{k}-m_{2})}{\sigma_{1}\sigma_{2}}, $$
(9)

where

$$\begin{aligned}& m_{1}= \sum_{k=1}^{K} \sum _{i=1}^{N_{k}-1} \frac{\tau^{k}_{i}}{N-K}, \end{aligned}$$
(10)
$$\begin{aligned}& m_{2}= \sum_{k=1}^{K} \sum _{i=2}^{N_{k}}\frac{ \tau^{k}_{i}}{N-K}, \end{aligned}$$
(11)
$$\begin{aligned}& \sigma_{1}=\sqrt{\sum_{k=1}^{K} \sum_{i=1}^{N_{k}-1} \frac{(\tau^{k}_{i} - m_{1})^{2}}{N-K}}, \end{aligned}$$
(12)

and

$$ \sigma_{2}=\sqrt{\sum_{k=1}^{K} \sum_{i=2}^{N_{k}} \frac{(\tau^{k}_{i} - m_{2})^{2}}{N-K}}. $$
(13)

An uncorrelated sequence of interevent times yields \(M=0\).

To examine the statistical significance of B value for each manager, we generated 103 sequences of interevent times from the exponential distribution whose mean was equal to that of the original data. Each synthesized sequence had the same length (i.e., N) as that of the original data. We calculated B for each synthesized sequence. We regarded that the value of B for the original data was significant if it was not included in the 95% confidential interval (CI) on the basis of the distribution generated by the 103 sequences corresponding to the Poisson process. We calculated the CI for M in the same manner except that we generated synthesized sequences by randomizing the original sequence of interevent times, instead of sampling sequences from the exponential distribution.

Figure 6 shows histograms of burstiness, B, and memory coefficient, M, for the managers. The average values of B and M for interevent times of formation changes in the J-League data were equal to 0.145 and −0.137, respectively. Those for the Bundesliga data set were equal to 0.022 and −0.120, respectively. For both data sets, the average values of B were positive, and those of M were negative. The fraction of managers yielding significantly positive and negative B values were equal to 0.385 and 0, respectively, for the J-League data. Those for the Bundesliga data were equal to 0.372 and 0.244, respectively. These results suggest that in both data sets, a moderate fraction of managers changed formations in a bursty manner. In the Bundesliga data, however, some managers changed formations more regularly than expected from the Poisson process. The fractions of significantly positive and negative M values were equal to 0.077 and 0.077, respectively, for the J-League data. Those for the Bundesliga data were equal to 0.026 and 0.026, respectively. In both data sets, the fractions of managers with significant M values were small, indicating that two consecutive interevent times were uncorrelated for a majority of managers.

Figure 6
figure 6

Distributions of burstiness, B , and memory, M , across managers. (a) B for the J-League data. (b) M for the J-League data. (c) B for the Bundesliga data. (d) M for the Bundesliga data. The colored bars correspond to managers who have statistically significant values.

Appendix D: Hierarchical Bayesian model

In the main text, we used the fraction of matches that a team won in a season to define the strength of the team. In this section, we analyze a model in which the strength of a team is assumed to be a latent variable. We used the hierarchical Bayesian ordered probit model combined with the Markov Chain Monte Carlo (MCMC) method [33]. The model is the same as that used in the main text except for the derivation of the team strength. We assumed that the prior of the strength of team i in a season, denoted by \(r_{i}\), obeyed the normal distribution with mean 0 and variance \(\sigma^{2}\). The priors of \(\beta_{\mathrm{f}}\), \(\beta_{\mathrm{ h}}\), \(\beta_{\mathrm{ w}}\), and \(\beta_{\ell}\) obeyed the normal distribution with mean 0 and variance 102. The prior of \(\sigma^{2}\) obeyed the uniform distribution on \([0, 10^{4}]\). We conducted MCMC simulations for four independent chains starting from the same prior distributions. The total iterate per chain was set to 25,000, and the first 5,000 iterates were discarded as transient. The thinning interval was set to 20 iterates. A final coefficient was regarded to be significant if the 95% credible interval did not bracket zero. We excluded the matches that were the first game in a season at least for either team. We performed the analysis using R 3.1.2 [30] and RStan package [44].

Table 11 summarizes the results obtained from the Bayesian probit model. For both data sets, the credible interval of the coefficient representing the effect of the formation change brackets zero. Therefore, we conclude that formation changes have not affected the probability of winning.

Table 11 Effects of variables on match results obtained from the hierarchical Bayesian ordered probit model

Appendix E: Cross-correlation analysis

To further investigate possible relationships between formation changes and match results, we measured the cross-correlation between the two. In this analysis, we did not exclude the first match in each season. We used the teams that played at least 100 matches. We set \(f_{i,t}=1\) if team i changes the formation in the tth match (\(2\le t\le T_{i}\)), where \(T_{i}\) is the number of matches played by team i, and \(f_{i,t}=0\) otherwise; \(w_{i,t} = 1\) if team i wins in the tth match, and \(w_{i,t}=0\) otherwise; \(\ell_{i,t} = 1\) if team i loses in the tth match, and \(\ell _{i,t}=0\) otherwise. We defined the cross-correlation between two time series \(\{x_{i,t}\}\) and \(\{y_{i,t}\}\) by

$$ \rho(x,y,\tilde{\tau}) = \frac{\sum_{i=1}^{N_{\mathrm{ team}}}\sum _{t=2}^{T_{i}-\tilde{\tau}} (x_{i,t+\tilde{\tau}}-\bar{x})(y_{i,t}-\bar {y})}{\sqrt{\sum_{i=1}^{N_{\mathrm{team}}}\sum_{t=2}^{T_{i}-\tilde{\tau }} (x_{i,t+\tilde{\tau}}-\bar{x})^{2}}\sqrt{\sum_{i=1}^{N_{\mathrm{ team}}}\sum_{t=2}^{T_{i}-\tilde{\tau}} (y_{i,t}-\bar{y})^{2}}}, $$
(14)

where \(\bar{x}=(1/N_{\mathrm{ team}})\times\sum_{i=1}^{N_{\mathrm{ team}}}\sum _{t=2}^{T_{i}-\tilde{\tau}}x_{i,t+\tilde{\tau}}/(T_{i}-\tilde{\tau }-1)\), \(\bar{y}=(1/N_{\mathrm{ team}})\times\sum_{i=1}^{N_{\mathrm{ team}}}\sum _{t=2}^{T_{i}-\tilde{\tau}}y_{i,t}/(T_{i}-\tilde{\tau}-1)\), \(N_{\mathrm{ team}}\) is the number of teams, and \(\tilde{\tau}\) is the lag. We measured the cross-correlation between formation changes and wins by

$$\begin{aligned}& \rho(f,w,\tilde{\tau}) \quad \mbox{if } \tilde{\tau} \geq0, \end{aligned}$$
(15)
$$\begin{aligned}& \rho(w,f,-\tilde{\tau}) \quad \mbox{if } \tilde{\tau} < 0. \end{aligned}$$
(16)

Replacing w by in Eqs. (15) and (16) defines the cross-correlation between formation changes and losses.

To examine the statistical significance of the cross-correlation obtained from the original data, we generated 103 randomized sequences of formation changes as follows. For given team i and positive lag \(\tilde{\tau}\), we randomly shuffled the original sequence of formation changes, \(\{f_{i,2+\tilde {\tau}}, \ldots, f_{i,T_{i}} \}\), by assigning 1 (i.e., formation change) to each match with the equal probability such that the number of 1s in the synthesized sequence was the equal to that in the original sequence. We generated a randomized sequence for each team. Then, we measured the cross-correlation between the randomized sequences of formation changes and \(\{ w_{i,2}, \ldots, w_{i,T_{i}-\tilde {\tau}}\}\) or \(\{ \ell_{i,2}, \ldots, \ell_{i,T_{i}-\tilde{\tau}}\}\) using Eq. (15). We repeated this procedure 103 times to obtain 103 cross-correlation values. The cross-correlation for the original data was considered to be significant for a given \(\tilde{\tau}\) if it was not included within the 95% CI calculated on the basis of the 103 correlation coefficient values for the randomized samples. We also examined the statistical significance of the cross correlation for a negative lag on the basis of 103 cross-correlation values between randomized sequences of \(\{ f_{i,2}, \ldots, f_{i,T_{i}-\tilde {\tau}}\}\) and \(\{w_{i,2+\tilde{\tau}}, \ldots, w_{i,T_{i}} \}\) or \(\{ \ell_{i,2+\tilde{\tau}}, \ldots, \ell_{i,T_{i}} \}\) using Eq. (16).

The cross-correlation measured for various lags is shown in Figure 7. The cross-correlation value was the largest in the absolute value at \(\tilde{\tau} = 1\). The effect of past matches on formation changes was mildly significant between \(\tilde{\tau} = 2\) and \(\tilde{\tau} \approx5\). The sign of the effect (i.e., positive or negative) was the same for different lag values, which complies with the concept of WSLS. When \(\tilde{\tau} \leq0\), the cross-correlation was insignificant or weakly significant even if \(\tilde{\tau} \approx0\). This result is suggestive of a causal relationship, i.e., a match result tends to cause a formation change.

Figure 7
figure 7

Cross-correlation between temporal patterns of formation changes and match results. Ranges between the dashed lines represent 95% CIs on the basis of the randomized sequences of formation changes. Cross-correlation between (a) formation changes and wins for the J-League data, (b) formation changes and losses for the J-League data, (c) formation changes and wins for the Bundesliga data, and (d) formation changes and losses for the Bundesliga data.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tamura, K., Masuda, N. Win-stay lose-shift strategy in formation changes in football. EPJ Data Sci. 4, 9 (2015). https://doi.org/10.1140/epjds/s13688-015-0045-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1140/epjds/s13688-015-0045-1

Keywords