- Open Access
Exploring the emergence and evolution of population patterns of leisure-time physical activity through agent-based modelling
International Journal of Behavioral Nutrition and Physical Activityvolume 15, Article number: 112 (2018)
Most interventions aiming to promote leisure-time physical activity (LTPA) at population level showed small or null effects. Approaching the problem from a systems science perspective may shed light on the reasons for these results. We developed an agent-based model to explore how the interactions between psychological attributes and built and social environments may lead to the emergence and evolution of LTPA patterns among adults.
The modeling process consisted of four stages: (1) conceptual model development, (2) formulation of the agent-based model, (3) parametrization and calibration, and (4) consistency and sensitivity analyses. The model represents a stylized community containing two types of agents: persons and LTPA sites. Persons interact with each other (proximal network and perceived community) and with the built environment (LTPA sites) over time. Decision-making is based on the person’s intention to practice LTPA, conditioned to the perceived environment. Each iteration is equivalent to one week and we assessed a period of 10 years.
The model was able to reproduce population temporal trends of intention and LTPA reported in the literature. Sensitivity analyses indicated that population patterns and trends of intention and LTPA were highly influenced by the relationship between a person’s behavior in the preceding week and his current intention, the person’s access to built and social environment, and the density of LTPA sites.
The proposed agent-based model is suitable to explore the emergence and evolution of LTPA patterns among adults, considering the dynamic interaction between individuals’ psychological attributes and the built and social environments in which they live. The model is available at https://doi.org/10.17605/OSF.IO/J2KAS.
There have been significant efforts to promote physical activity in both the leisure and transportation domains [1,2,3,4]. However, systematic reviews and meta-analyses [5,6,7,8,9,10,11] have revealed that population-based initiatives have produced small, null, or even inconsistent effects. Part of the problem is that most of these initiatives focused on individual behavior solely and did not consider the dynamic relations among individuals and between individuals and their environments .
Despite the accumulated knowledge on what encourages or prevents people from being physically active , little progress has been observed on increasing population levels of physical activity. One of the reasons is that knowledge and action on the determinants of individual behavior do not always lead to population behavior changes . Leisure-time physical activity (LTPA) exemplifies this situation. Evidence about plausible determinants and interventions to change LTPA levels has grown exponentially . However, and despite the growing investment in population-based initiatives to promote LTPA in the last decade, temporal trends of LTPA largely remained stable or increased marginally [15,16,17,18,19,20,21,22].
Kohl and collaborators  advocate the use of systems science methods to address this mismatch and widen our understanding on population-level physical activity. Under the systems science perspective, population patterns and changes of physical activity cannot be inferred only by aggregating individual behaviors. These patterns and changes derive from a complex adaptive system, composed of dynamic interactions among heterogeneous elements, which are, at the same time, autonomous and interdependent. Incorporating the systems science framework would help researchers and practitioners to understand how population patterns and trends of physical activity emerge and design effective interventions.
One method for modeling complex adaptive systems is agent-based modeling. In this modeling approach, a system is represented as a composition of autonomous entities, called agents, and the environment in which they live. Each agent possesses decision-making capabilities, determined by a set of rules on how to act in face of the agents’ interaction with the surrounding environment and other agents [23, 24]. Agent-based modeling has been increasingly utilized to investigate how dynamic processes involving built and social environments affect population health , including health behaviors and chronic diseases [26, 27].
Reasons to develop an agent-based model of population patterns of LTPA are fourfold. First, LTPA is consistently associated with longevity and better health  and is one of the pillars for most population-based physical activity promotion initiatives . Second, agent-based models can help researchers understand the interrelations and impacts of factors at different levels — from individual-level factors to environments and policies — thus facilitating the identification of innovative ways to intervene. Third, all agent-based models related to population levels of physical activity published so far focused on transport-related physical activity . Finally, literature and data on LTPA are vast and rich, in quantity and range, providing a solid modelling foundation.
Our goal was to develop an agent-based model that can allow researchers to explore the emergence and evolution of population patterns of LTPA among adults, taking into consideration the interaction between individuals’ psychological attributes and the built and social environments in which they live.
The modeling process was composed of four stages: (1) conceptual model development, (2) formulation of the agent-based model, (3) parametrization and calibration, and (4) consistency and sensitivity analyses.
Conceptual model development
The conceptual model depicts the main psychological and environmental (built and social) elements and processes that may play a role in the emergence and evolution of population patterns of LTPA in adult populations.
A detailed account of the conceptual model development is available elsewhere . First, we drafted a version based on the expertise of all authors. Second, we iteratively updated the conceptual model using information obtained from a literature review on psychological attributes, and built and social environments related to LTPA. Third, an intermediate version of the model was assessed by 18 experts in at least one of the topics approached by the conceptual model. The final version of the model integrated expert assessments and additional information found in the literature.
Figure 1 summarizes the conceptual model. We did not intend to produce a comprehensive framework of factors influencing physical activity behavior but a rather focused one to inform and support the development of the agent-based model, which was delimited a priori to encompass the main psychological and environmental elements and processes shaping population patterns of LTPA. Therefore, some aspects, such as demographic attributes, were not included as we felt they were not highly relevant to our current research questions.
The agent-based model was developed based on the conceptual model. We used the Overview, Design Concepts, and Details + Decision (ODD+D) protocol  to formulate the structure and content of the agent-based model. Below we present a brief description of the model main aspects. For the full description, see Additional file 1.
The model has two types of agents: persons and LTPA sites. A person has four attributes: location, intention (ranging from 0.03 to 0.97), one favorite type of LTPA (among 10 available in the model), and behavior (did or did not practice LTPA in the previous week). An LTPA site also has four attributes: location, quality (ranging from 0 to 1), and number and types of activities available (among 10 available in the model).
A square grid (50 × 50 patches as default) with non-contiguous sides represents the physical space, not emulating any particular location or community. Patches are characterized as being or not a LTPA site. By default, the model contains two thousand people. Persons are placed in the grid in patches that are not LTPA sites. Time is discrete, and each iteration is equivalent to one week, a timeframe frequently used to investigate physical activity behavior and that, therefore, makes easier to calibrate and compare our model against empirical data. We assessed a period of 10 years (520 iterations).
Every week, each person’s intention is updated based on his behavior and the behavior of his proximal network and perceived community in the previous week. Then, each person decides whether he will practice LTPA during the current week, based on the new level of intention and conditional to the perceived built environment. The entire process is synchronic (i.e., every person updates the same attributes at the same time) and happens in a weekly basis.
At the initialization, a proportion of patches is defined as LTPA sites and randomly positioned in the grid. Each site gets a quality score drawn from a normal distribution. Next, the number of LTPAs available in each site is drawn from a uniform distribution. Then, the types of LTPAs available in each site are randomly drawn from a list of 10 possible activities using a uniform distribution.
Persons are placed randomly over patches that do not represent LTPA sites. Every person sets, randomly, none or one of the 10 available LTPA types as his favorite, drawing from a uniform distribution. Next, each person incorporates in his memory a list of LTPA sites within his perception radius — a delimited area within the grid that the person can observe — and their attributes:
Access (comprising factors such as traffic, safety, physical proximity, cost and ease of transportation to it), which is represented in the model as the Euclidian distance between the person and the site (the longer the distance, the worse the access);
Whether his favorite LTPA is available at the site.
Each person then calculates and stores in his memory the sites’ perceived utility, obtained by Eq. 1:
In which us,i is the utility of site s as perceived by person i, vi,s is the subjective assessment (scaling factor) given by person i to site s, qs is the quality of site s, zi,s is the access of person i to site s, and mi,s shows whether the site s offers the favorite LTPA of person i. A subjective assessment scaling factor is assigned to each person–LTPA-site dyad by drawing from a normal distribution.
Next, the level of intention is set for every person by drawing from two uniforms distributions, each one applied to one of two subsets of the population, split randomly. This procedure allows more flexibility for initializing the population distribution of intention. For this work, the first subset encompasses 25% of the population and has level of intention ranging from 0.03 to 0.30, whereas the level of intention of the other 75% ranges from 0.31 to 0.97.
Each person is also assigned a behavior (i.e., did or did not practice LTPA in the past week), with probability equal to his intention. Persons who do not have any LTPA site within their perception radius are assigned to not having practiced the behavior in the past week.
Finally, a proximal network and a perceived community are defined for each person, as follows:
Proximal network: represents those with whom the person has close relationships (such as friends and relatives). The network is formed in a two-stage process. First, the k closest persons in the grid are selected. Then, each link has a probability p, set by the modeler (0.15 for this work), to be exchanged for a link with any other person outside the initial proximal network. This procedure aims to ensure that the proximal network is mostly formed by people with similar social and environmental influences, but with some degree of variability.
Perceived community: formed by those people within the person’s perception radius. Each person infers the social norm of the entire community looking at the behavior of the perceived community.
Initial configurations of the same scenario may differ between replications as they mostly rely on stochastic processes. Figure 2 exemplifies a grid after the initialization.
Changes over time
In our model, persons tend to maintain their habitual intention while seeking to adapt it according to the behavior observed from their proximal network and community. Moreover, the translation of intention into behavior is conditional to the perceived built environment.
The person’s decision-making process is grounded in an extension of the Continuous Opinions and Discrete Actions model , which posits that a discrete action (in this case, practicing LTPA or not) is a function of a continuous internal opinion (in this case, intention), using Bayes’ theorem. Table 1 shows the equations to update the persons' intention every week. The level of intention in a given week (Eq. 7) depends on the behavior of those in the proximal network (Eq. 2) and perceived community (Eq. 3), the person’s behavior in the previous week (Eq. 4), current level of intention (Eq. 5), and the highest perceived utility amongst the LTPA sites in the person’s perception radius (Eq. 6).
In Eqs. 2, 3, 4, the log-odds of α can be interpreted as how confident someone is that the behavior of his neighbors (or his own behavior, in Eq. 4) represents the best option. Different values for α are possible, allowing for different weights of the person’s and neighbors’ behavior in affecting future levels of intention.
Lastly, level of intention in week t is transformed into the probability of practicing LTPA in that week by Eq. 8:
In which Prob(bi,t) is the probability of person i practicing LTPA in week t, and ei,t is the exponential of person i’s intention in week t. Prob(bi,t) is limited between 0.03 and 0.97 instead of 0 and 1 to account for the possibility that one could adopt the opposite behavior regardless of his current intention (for instance, being ill [for those with 0.97] or in vacation [for those with 0.03]). This probability is used to set each person’s behavior in week t. The entire cycle reinitiates in week t + 1 beginning with Eq. 2.
The model’s algorithm was implemented on NetLogo version 6.0.3 . A verification protocol was followed to ensure the accuracy of the algorithm, consistency with the conceptual model, and avoid errors, omissions and bugs. The protocol included continuous reviewing of spelling, syntax and commands; monitoring of unexpected elements, actions, calculation, and outputs; and stress tests (i.e., using inputs and parameter values out the plausible range) [34, 35].
Spatial distribution of LTPA practice, population distribution of intention, and proportion of people practicing LTPA over time can be visualized using the model’s user interface and/or recorded for future analyses.
Parametrization and calibration
We defined initial values and domains of each parameter in the model following the premises of pattern-oriented modeling , in which the patterns observed in real systems are used as a reference to obtain models that are more realistic in structure and, therefore, more useful and generalizable. Two population patterns were used as guidance for parametrization and calibration: temporal trends of LTPA (> 0 min/week) and population distribution of intention to practice LTPA.
Studies from around the world involving adult populations show temporal trends of LTPA are stable or increasing by very little, generally less than one percentage point a year, in periods varying from five to ten years (the period simulated in the model), with prevalence ranging from 35 to 50% [15,16,17,18,19,20].
Our search for temporal trends in the population distribution of intention to practice LTPA did not return any relevant information. We then searched for population-based cross-sectional studies and identified studies on population distributions of behavioral change stage — construct of the Transtheoretical Model  that is a proxy of intention. In general, distributions were U shaped, with higher proportions of people in the pre-contemplation (20–35%) and maintenance (25–45%) stages. Contemplation and preparation stages took turns as the third largest proportion (around 15% each), while the action stage presented the lowest values (around 5%) [37,38,39,40].
Therefore, the model was expected to reproduce scenarios in which LTPA was stable or presented slight increase over time, and a U-shaped distribution of intention. To monitor these patterns within the model, we extracted the proportion of people practicing LTPA and people with low (0.03 ≤ intention < 0.25), intermediate (0.25 ≤ intention ≤0.75), and high (0.75 < intention ≤0.97) intention.
After defining the expected population patterns for the model, we defined initial values and domains of parameters for which we found empirical data: size of proximal network [41,42,43] and perceived community [42, 44, 45], and amount/density of LTPA sites [46,47,48].
No empirical data was found on the size of persons’ perception radius. Therefore, this parameter was calibrated to match values reported in the literature for the size of perceived community [42, 44, 45] and amount/density of LTPA sites in persons’ neighborhoods [46,47,48]. Final value was set to 9.0 (domain = 0 to 50).
For the other parameters, we defined a plausible initial value and domain that would lead to the expected time trends of LTPA practice and the expected population levels of intention, since empirical data could not be found for parametrization. These parameters were calibrated simultaneously by trial and error, a procedure considered sufficient at the first stages to proceed with the analyses. In later stages, we identified the most relevant parameters and assessed sensitivity of outputs to input values.
Most of the empirical data used for parametrization and calibration came from urban settings within high-income Western countries. However, wide but plausible domains were set for all parameters to accommodate contexts not initially covered by the available data and allow thorough exploration of the parameters during sensitivity analyses.
Additional file 1: Table S1.2 displays all the parameters in the model and their domains and initial values.
Consistency and sensitivity analysis
Twenty sets were generated for each m number of replications (1, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100). From every replication we extracted the proportion of people practicing LTPA and with low, intermediate, and high intention. Each set contained a distribution with m values for each output. Using the Vargha-Delaney A Test  we identified whether the 20 distributions including m replications were consistent against each other.
The Vargha-Delaney A Test is not a hypothesis testing procedure but still generates an effect size measure that takes values between 0 and 1, with a value of 0.5 indicating that the medians of two distributions with m replications are similar, showing a high consistency of results (i.e., less uncertainty due to stochasticity). Scores close to 0.44 or 0.56 indicate small differences, values around 0.36 or 0.64 express medium effect size, and scores lower than 0.29 or higher than 0.71 indicate large differences [50, 51].
For each m replications, median and maximum A score were calculated from the scores obtained comparing the first and the remaining 19 distributions generated. Line graphs were plotted to identify the minimum required number of replications that result in low uncertainty due to stochasticity.
We conducted individual and global sensitivity analyses. For individual sensitivity analyses, we used the parameter robustness technique, in which every parameter of interest is disturbed individually, while the remaining ones are kept in their initial values [49, 50]. Fourteen parameters (two related to personal attributes, four to social environment, and eight to built environment) were investigated (Additional file 2: Table S2.1). Temporal trends of LTPA and distribution of intention obtained from the scenarios with the new conditions were compared against the data obtained from the baseline scenario through the Vargha-Delaney A Test . Scores close to 0.5 indicated that temporal trends were robust (i.e., not very sensitive) to alterations of the parameter’s values . The farther the score is from 0.5, the more sensitive the model is to that particular parameter. Scores higher than 0.71 or lower than 0.29 indicate significant sensitivity [50, 51]. Line graphs and boxplots were used to provide summaries and visualize results.
For the global sensitivity analysis, selected parameters of interest were simultaneously disturbed . We first selected the most influential parameters from the individual sensitivity analysis (Additional file 2: Table S2.2). Next, Latin Hypercube Sampling technique was employed to generated 100 scenarios by drawing values from each selected parameter to optimize the exploration and coverage of the whole domain of all parameters simultaneously . Partial rank correlation coefficients and scatter plots were generated to analyze the correlation between each parameter and outputs that can be observed despite the simultaneous changes in all other parameters . Results were summarized by minimum, maximum, mean, and standard-deviation. Effect sizes were interpreted as small (~ 0,2), average (~ 0,5) or high (~ 0,8) .
Outputs analyzed and statistical package
For all consistency and sensitivity analyses, we obtained and analyzed the proportion of people practicing LTPA and people with low, intermediate, and high intention at every 52 weeks in a total of 624 weeks, which is equivalent to yearly data for 12 years.
The statistical software R, version 3.2.2, and spartan package, version 2.3, were used to conduct the analyses.
Availability of data and materials
The model code and user interface, R codes, original outputs of all replications, and full results generated for this work are available at https://doi.org/10.17605/OSF.IO/J2KAS. Parts of this information are also described in Additional files 1, 2, 3 and 4.
Figure 3 shows temporal trends of persons practicing LTPA and with low, intermediate and, high intention, generated from the scenario using the parameter values in Additional file 1: Table S1.2. This scenario was able to reproduce the temporal trends reported in the literature, showing a stable prevalence of LTPA of approximately 48% throughout the years, while the population distribution of intention exhibited a U shape, with increasing proportions of people with high and low intention over time.
Results from the consistency analysis were summarized in 169 charts and 26 spreadsheets. Figure 4 displays the median and maximum scores from the Vargha-Delaney A Test for each number of replications tested, obtained at the 520th iteration (equivalent to the 10th year). Plots and spreadsheets of all other years (available at https://doi.org/10.17605/OSF.IO/J2KAS) have returned results similar to those portrayed in Fig. 4.
After 80 replications, the reduction of both the median and maximum A score became stable. Therefore, 80 replications were deemed as the minimum required to obtain time trends of LTPA practice and intention levels that are reliable and less influenced by stochasticity.
In all the other analyses conducted for this work, 80 replications were made for each scenario.
Additional file 3 provides a summary of the observed results from the individual sensitivity analysis. From the 14 parameters tested, two had high influence on the proportion of people practicing LTPA: the size of the person’s perception radius, and the proportion of patches in the grid that are LTPA sites. Three parameters (the size of the perception radius, the proportion of patches that are LPTA sites, and the influence of the person’s behavior in the previous week over his current intention — αb term in Eq. 4) had high influence on the proportion of people with low, intermediate, and high levels of intention.
Six parameters were selected for the global sensitivity analysis, two related to personal attributes (influence of the person’s behavior in the previous week over his current intention, size of the person’s perception radius), two to social environment (influence of the proximal network’s behavior over the person’s intention [αp term in Eq. 2], influence of the perceived community’s behavior over the person’s intention [αc term in Eq. 3]), and two to built environment (proportion of LTPA sites, mean quality score of LTPA sites) (Additional file 2: Table S2.2).
Table 2 displays a summary of the partial rank correlation coefficients, obtained from the global sensitivity analysis using the most influential parameters. Additional file 4, also generated from the global sensitivity analysis, exemplifies scatter plots obtained at the 520th iteration (equivalent to the 10th year) assessed in order to identify non-linear effects.
Global sensitivity analyses indicated that the model is highly sensitive to three parameters: the influence of the person’s behavior in the previous week over his current intention, the size of the person’s perception radius, and the proportion of patches in the grid that are LTPA sites (Table 2).
The influence of the person’s behavior in the previous week over his current intention had a high impact on the proportion of people in each level of intention. The stronger the influence of the previous behavior over current intention, the more persons presented high or low intention, while fewer persons remained in the intermediate range. However, the influence over the proportion of people practicing LTPA was small (Table 2).
The size of a person’s perception radius and the proportion of LTPA sites had a high impact on temporal trends of LTPA practice and a moderate to high effect on intention levels. The wider the perception radius and the higher the proportion of LTPA sites, the higher the proportion of people practicing LTPA and with intermediate and high intention, and the lower the proportion of people with low intention (Table 2). However, the relationship was non-linear. Small increments when the radius is narrower led to greater changes in the proportion of people practicing LTPA and in each level of intention (Additional file 4: Figure S4.2). As for the proportion of LTPA sites, changes in the proportion of people practicing LTPA and in each level of intention were observed only at the extremities of the values domain (Additional file 4: Figure S4.5).
The influence of the proximal network’s behavior over the person’s intention had a low impact on the proportion of people practicing LTPA and on the proportion with high intention, whereas an average impact on the proportion of those with low and intermediate intention. The influence of the perceived community’s behavior over a person’s intention had a low impact on the four proportions. For both parameters, the stronger the effect of the social environment, the higher the proportion of people with low intention, and the lower the proportion of those with intermediate intention, an unexpected result (Table 2).
The model displayed low sensitivity to the mean quality score of LTPA sites. As the quality increased, the proportion of people practicing LTPA and with high intention increased only slightly (Table 2).
Our goal was to develop an agent-based model that can enable researchers to explore the emergence and evolution of population patterns of LTPA among adults, emerging from the interaction between the individuals’ psychological attributes and the built and social environments in which they live. Analyses conducted with the model have shown that time trends of LTPA practice and population distribution of levels of intention are similar those reported in literature [15,16,17,18,19,20,21,22, 37,38,39,40], which suggests the model was capable of capturing the behavior of the emulated phenomenon. Our analyses also identified the elements and mechanisms that significantly influence the temporal trends in the model, particularly the influence of the person’s behavior in the previous week over his current intention, the size of the person’s perception radius, and the proportion of patches in the grid that are LTPA sites.
The model was based on a conceptual model designed from a literature review and expert assessment  in order to maintain coherence and consistency with the best evidence available in the field. It is strongly based on the premises of systems thinking, expressly suggesting the dynamic mechanisms and relationships between individual psychological attributes and social and built environments. In our model, psychological attributes are considered the strongest proximal determinants of LTPA, a relationship that is dynamically moderated by the built environment — moderation that depends on the person’s psychological attributes — and influenced by both the social environment and the behavior itself. The literature review, expert opinions , and the model’s initial results suggest that this is a plausible representation of reality.
Three elements and mechanisms exhibited stronger influence on time trends of people practicing LTPA and levels of intention: the influence of the person’s behavior in the previous week over his current intention, size of the person’s perception radius, and proportion of LTPA sites in the model. Other three elements and mechanisms had lower effect: proximal network’s and perceived community’s behaviors influence on the person’s intention, and mean quality score of LTPA sites. These six parameters — particularly the first three — should receive special attention when calibrating the model for future works.
Individual and global sensitivity analyses shown some unexpected results and non-linear dynamics. One example is the effect of the social environment. The influence of the proximal network’s and perceived community’s behavior over the person’s intention presented small effect sizes, but in an undesired way: the stronger the social influence, the higher the proportion of people with low intention to practice LTPA. This may have happened because, on average, persons had more contact with others who did not practice LTPA, as the prevalence of LTPA was usually less than 50% for the investigated scenarios. This situation is, however, common to several locations [15, 16, 18,19,20,21,22] and should be considered when planning campaigns and initiatives promoting LTPA using the social environment as one of the elements.
Another example is the positive effect of increasing the proportion of LTPA sites over the temporal trends of people practicing LTPA and levels of intention. This positive effect can be split into two distinct phases. When a certain proportion of sites is achieved, additional positive effects on LTPA practice may depend on a much stronger and sustained investment on infrastructure to translate it into significant changes on LTPA practice. Indeed, at least one meta-analysis , four systematic reviews of quantitative studies [54,55,56,57] and one of qualitative studies  shown positive associations between LTPA practice and distance to or density of LTPA sites. However, a possible non-linear relationship between density of LTPA sites and LTPA prevalence is still to be investigated.
The most important contribution of our agent-based model is the systems perspective it provides on the issue, being more focused on the system’s structure [59, 60]. As a general model, the usefulness of this model to LTPA promotion does not hinge on its use as a tool for evaluation of policy interventions in a particular location, but as an additional means to obtain deeper comprehension on the factors influencing population patterns of LTPA and how they interact shaping and sustaining the observed patterns. It is noteworthy, however, that the model can, in principle, be extended and validated to represent real locations and test policy options specific to these places.
There is space for improvements in future versions of the model, especially with the advent of more evidence that can inform the model structure and parameters. Future versions could:
Include perception and decision capacities to LTPA sites, allowing them to adapt to the characteristics of people in their surroundings;
Allow certain parameters to vary among persons, such as perception radius or the relative weights given to quality, access, and availability of favorite activity when raking LTPA sites;
Break down intention into its psychological precedents, such as attitude and self-efficacy;
Consider attributes of LTPA sites, for example, financial cost of use and distance, separately.
However, the current version of the model is already capable of reproducing population time trends of LTPA and distributions of intention observed in real settings and can be used to explore some of the dynamic mechanisms generating them.
The agent-based model developed in this work is appropriate to investigate the emergence and evolution of population patterns of LTPA among adults, emerging from the interaction between the individuals’ psychological attributes and the built and social environments in which they live. Initial results showed that three elements and mechanisms displayed stronger influence on time trends of people practicing LTPA and levels of intention within the model: the influence of the person’s behavior in the previous week over his current intention, size of the person’s perception radius, and proportion of patches in the model that are LTPA sites.
Asides from improving the model, future work includes testing scenarios and interventions — either hypothetical or equivalent to real situations — to gain a better understanding of the conditions that generate either observed or desired population patterns. The model is available to download (https://doi.org/10.17605/OSF.IO/J2KAS) and can be used, adapted, and improved by others for this purpose. We expect that work deriving from our model can inform researchers and policymakers on the design of more effective population-level physical activity interventions.
Population-based initiatives will probably continue to present results below the desired level until we seek to understand the dynamics and structures that support population patterns of LTPA. The results achieved by this work show that it is possible to approach population patterns of LTPA in a way that is more in line with a systems science perspective and with the idea that physical activity is a complex, multidimensional and multidetermined behavior.
Leisure-time physical activity
Brasil. Ministério da Saúde. Plano de ações estratégicas para o enfrentamento das doenças crônicas não transmissíveis (DCNT) no Brasil 2011-2022. Brasília: Ministério da Saúde; 2011.
Global Advocacy Council for Physical Activity. The Toronto charter for physical activity: a global call to action. Toronto: International Society For Physical Activity And Health; 2010.
World Health Organization. A guide for population-based approaches to increasing levels of physical activity: implementation of the WHO global strategy on diet, physical activity and health. Geneva: WHO Press; 2007.
World Health Organization Regional Office for Europe. Physical activity strategy for the WHO European region 2016–2025. Copenhagen: WHO Europe; 2015.
Abioye AI, Hajifathalian K, Danaei G. Do mass media campaigns improve physical activity? A systematic review and meta-analysis. Arch Public Health. 2013;71(1):20.
Baker PR, Francis DP, Soares J, Weightman AL, Foster C. Community wide interventions for increasing physical activity. Cochrane Database Syst Rev. 2015;1:CD008366.
Brand T, Pischke CR, Steenbock B, Schoenbach J, Poettgen S, Samkange-Zeeb F, et al. What works in community-based interventions promoting physical activity and healthy eating? A review of reviews. Int J Environ Res Public Health. 2014;11(6):5866–88.
Cavill N, Bauman A. Changing the way people think about health-enhancing physical activity: do mass media campaigns have a role? J Sports Sci. 2004;22(8):771–90.
Cleland CL, Tully MA, Kee F, Cupples ME. The effectiveness of physical activity interventions in socio-economically disadvantaged communities: a systematic review. Prev Med. 2012;54(6):371–80.
Conn VS, Hafdahl AR, Mehr DR. Interventions to increase physical activity among healthy adults: meta-analysis of outcomes. Am J Public Health. 2011;101(4):751–8.
Heath GW, Parra DC, Sarmiento OL, Andersen LB, Owen N, Goenka S, et al. Evidence-based intervention in physical activity: lessons from around the world. Lancet. 2012;380(9838):272–81.
Kohl HW, Craig CL, Lambert EV, Inoue S, Alkandari JR, Leetongin G, et al. The pandemic of physical inactivity: global action for public health. Lancet. 2012;380(9838):294–305.
Epstein JM. Agent-based computational models and generative social science. Complexity. 1999;4(5):41–60.
Rhodes RE, Nasuti G. Trends and changes in research on the psychology of physical activity across 20 years: a quantitative analysis of 10 journals. Prev Med. 2011;53(1–2):17–23.
Cruz MC. Tendência da prática de atividade física no lazer no Brasil (2006–2013) [Thesis]. Sao Paulo: Universidade de Sao Paulo; 2015.
Juneau CE, Potvin L. Trends in leisure-, transport-, and work-related physical activity in Canada 1994-2005. Prev Med. 2010;51(5):384–6.
Knuth AG, Hallal PC. Temporal trends in physical activity: a systematic review. J Phys Act Health. 2009;6(5):548–59.
Mielke GI, Hallal PC, Malta DC, Lee IM. Time trends of physical activity and television viewing time in Brazil: 2006-2012. Int J Behav Nutr Phys Act. 2014;11:101.
Nunes AP, Luiz OC, Barros MB, Cesar CL, Goldbaum M. Domínios de atividade física e escolaridade em São Paulo, Brasil: estudo transversal seriado, 2003 e 2008. Cad Saude Publica. 2015;31(8):1743–55.
Román-Viñas B, Serra-Majem L, Ribas-Barba L, Roure-Cuspinera E, Cabezas C, Vallbona C, et al. Trends in physical activity status in Catalonia, Spain (1992-2003). Public Health Nutr. 2007;10(11A):1389–95.
Silva IC, Knuth AG, Mielke GI, Azevedo MR, Goncalves H, Hallal PC. Trends in leisure-time physical activity in a southern Brazilian city: 2003-2010. J Phys Act Health. 2014;11(7):1313–7.
Stamatakis E, Chaudhury M. Temporal trends in adults’ sports participation patterns in England between 1997 and 2006: the health survey for England. Br J Sports Med. 2008;42(11):901–8.
Bonabeau E. Agent-based modeling: methods and techniques for simulating human systems. Proc Natl Acad Sci U S A. 2002;99(Suppl. 3):7280–7.
Gilbert N, Troitzsch K. Simulation for the social scientist. New York: Open University Press; 2005.
Auchincloss AH, Roux AV. A new tool for epidemiology: the usefulness of dynamic-agent models in understanding place effects on health. Am J Epidemiol. 2008;168(1):1–8.
Luke DA, Stamatakis KA. Systems science methods in public health: dynamics, networks, and agents. Annu Rev Public Health. 2012;33:357–76.
Nianogo RA, Arah OA. Agent-based modeling of noncommunicable diseases: a systematic review. Am J Public Health. 2015;105(3):e20–31.
Dishman R, Heath G, Lee IM. Physical activity epidemiology. 2nd ed. Champaign, IL: Human Kinetics; 2012.
World Health Organization. Global action plan on physical activity 2018–2030: more active people for a healthier world. Geneva: WHO Press; 2018.
Garcia LMT, Diez Roux AV, Martins ACR, Yang Y, Florindo AA. Development of a dynamic framework to explain population patterns of leisure-time physical activity through agent-based modeling. Int J Behav Nutr Phys Act. 2017;14:111.
Müller B, Bohn F, Dreßler G, Groeneveld J, Klassert C, Martin R, et al. Describing human decisions in agent-based models–ODD+ D, an extension of the ODD protocol. Environ Model Softw. 2013;48:37–48.
Martins ACR. Continuous opinions and discrete actions in opinion dynamics problems. Int J Mod Phys C. 2008;19(4):617–24.
Wilensky U. NetLogo. 1999 [12 June 2018]. Available from: http://ccl.northwestern.edu/netlogo.
Railsback SF, Grimm V. Agent-based and individual-based modeling: a practical introduction. Princeton: Princeton University Press; 2012.
Wilensky U, Rand W. An introduction to agent-based modeling: modeling natural, social, and engineered complex systems with NetLogo. Cambridge: MIT Press; 2015.
Prochaska JO, Redding CA, Evers KE. The transtheoretical model and stages of change. In: Glanz K, Rimer BK, Viswanath K, editors. Health behavior and health education: theory, research, and practice. San Francisco: Wiley; 2008. p. 97–121.
Dumith SC, Gigante DP, Domingues MR. Stages of change for physical activity in adults from southern Brazil: a population-based survey. Int J Behav Nutr Phys Act. 2007;4:25.
Laforge RG, Velicer WF, Richmond RL, Owen N. Stage distributions for five health behaviors in the United States and Australia. Prev Med. 1999;28(1):61–74.
Bull FC, Eyler AA, King AC, Brownson RC. Stage of readiness to exercise in ethnically diverse women: a U.S. survey. Med Sci Sports Exerc. 2001;33(7):1147–56.
Kearney JM, Graaf C, Damkjaer S, Engstrom LM. Stages of change towards physical activity in a nationally representative sample in the European Union. Public Health Nutr. 1999;2(1A):115–24.
Roberts SGB, Dunbar RIM, Pollet TV, Kuppens T. Exploring variation in active network size: constraints and ego characteristics. Soc Netw. 2009;31(2):138–46.
Stiller J, Dunbar RIM. Perspective-taking and memory capacity predict social network size. Soc Netw. 2007;29(1):93–104.
Carroll J. Americans satisfied with number of friends, closeness of friendships. 2004 [18 Sept 2015]. Available from: https://news.gallup.com/poll/10891/americans-satisfied-number-friends-closenessfriendships.aspx. Accessed 16 Oct 2018.
Hill RA, Dunbar RIM. Social network size in humans. Hum Nat. 2003;14(1):53–72.
Kanai R, Bahrami B, Roylance R, Rees G. Online social network size is reflected in human brain structure. Proc Biol Sci. 2012;279(1732):1327–34.
Hino AAF. Medidas objetivas e percebidas do ambiente do bairro e sua associação com a atividade física de lazer em adultos de Curitiba [Dissertation]. Curitiba: Universidade Federal do Paraná; 2014.
Halonen JI, Stenholm S, Kivimaki M, Pentti J, Subramanian SV, Kawachi I, et al. Is change in availability of sports facilities associated with change in physical activity? A prospective cohort study. Prev Med. 2015;73:10–4.
Estabrooks PA, Lee RE, Gyurcsik NC. Resources for physical activity participation: does availability and accessibility differ by neighborhood socioeconomic status? Ann Behav Med. 2003;25(2):100–4.
Alden K, Read M, Andrews PS, Timmis J, Coles M. Applying spartan to understand parameter uncertainty in simulations. R J. 2014;6(2):63–80.
Read M, Andrews PS, Timmis J, Kumar V. Techniques for grounding agent-based simulations in the real domain: a case study in experimental autoimmune encephalomyelitis. Math Comput Model Dyn Syst. 2012;18(1):67–86.
Vargha A, Delaney HD. A critique and improvement of the CL common language effect size statistics of McGraw and Wong. J Educ Behav Stat. 2000;25(2):101–32.
Sullivan GM, Feinn R. Using effect size - or why the p value is not enough. J Grad Med Educ. 2012;4(3):279–82.
Duncan MJ, Spence JC, Mummery WK. Perceived environment and physical activity: a meta-analysis of selected environmental characteristics. Int J Behav Nutr Phys Act. 2005;2:11.
Wendel-Vos W, Droomers M, Kremers S, Brug J, van Lenthe F. Potential environmental determinants of physical activity in adults: a systematic review. Obes Rev. 2007;8(5):425–40.
Arango CM, Paez DC, Reis RS, Brownson RC, Parra DC. Association between the perceived environment and physical activity among adults in Latin America: a systematic review. Int J Behav Nutr Phys Act. 2013;10:122.
Bauman AE, Reis RS, Sallis JF, Wells JC, Loos RJ, Martin BW, et al. Correlates of physical activity: why are some people physically active and others not? Lancet. 2012;380(9838):258–71.
Humpel N, Owen N, Leslie E. Environmental factors associated with adults' participation in physical activity: a review. Am J Prev Med. 2002;22(3):188–99.
McCormack GR, Rock M, Toohey AM, Hignell D. Characteristics of urban parks associated with park use and physical activity: a review of qualitative research. Health Place. 2010;16(4):712–26.
Meadows DH. Thinking in systems: a primer. Hartford: Chelsea Green Publishing; 2008.
El-Sayed AM, Scarborough P, Seemann L, Galea S. Social network analysis and agent-based modeling in social epidemiology. Epidemiol Perspect Innov. 2012;9(1):1.
Garcia LMT. Exploring population patterns and trends of leisure-time physical activity through agent-based modelling. Open Science Framework. 2018. Available from:. https://doi.org/10.17605/OSF.IO/J2KAS.
LG was supported by a scholarship from the Brazilian Coordination for the Improvement of Higher Education Personnel (grant 1406604). LG has worked under the auspices of the Centre for Diet and Activity Research (CEDAR), a UKCRC Public Health Research Centre of Excellence which is funded by the British Heart Foundation, Cancer Research UK, Economic and Social Research Council, Medical Research Council, the National Institute for Health Research, and the Wellcome Trust (grant MR/K023187/1). ADR is supported by the Wellcome Trust (grant 205177/Z/16/Z). AAF receives a fellowship from the Brazilian National Council for Scientific and Technological Development (CNPq) (grant 306635/2016–0).
Availability of data and materials
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Full description of the agent-based model. (PDF 321 kb)
Parameters and values investigated in sensitivity analyses. (PDF 85 kb)
Summary of results of individual sensitivity analysis. (PDF 113 kb)
Results of global sensitivity analysis (year 10). (PDF 712 kb)