Skip to main content

A decade of child pedestrian safety in England: a Bayesian spatio-temporal analysis



Child pedestrian injury is a public health and health equality challenge worldwide, including in high-income countries. However, child pedestrian safety is less-understood, especially over long time spans. The intent of this study is to understand factors affecting child pedestrian safety in England over the period 2011–2020.


We conducted an area-level study using a Bayesian space-time interaction model to understand the association between the number of road crashes involving child pedestrians in English Local Authorities and a host of socio-economic, transport-related and built-environment variables. We investigated spatio-temporal trends in child pedestrian safety in England over the study period and identified high-crash local authorities.


We found that child pedestrian crash frequencies increase as child population, unemployment-related claimants, road density, and the number of schools increase. Nevertheless, as the number of licensed vehicles per capita and zonal-level walking/cycling increase, child pedestrian safety increases. Generally, child pedestrian safety has improved in England since 2011. However, the socio-economic inequality gap in child pedestrian safety has not narrowed down. In addition, we found that after adjusting for the effect of covariates, the rate of decline in crashes varies between local authorities. The presence of localised risk factors/mitigation measures contributes to variation in the spatio-temporal patterns of child pedestrian safety.


Overall, southern England has experienced more improvement in child pedestrian safety over the last decade than the northern regions. Our study revealed socio-economic inequality in child pedestrian safety in England. To better inform safety and public health policy, our findings support the importance of a targeted system approach, considering the identification of high-crash areas while keeping track of how child pedestrian safety evolves over time.

Peer Review reports


Road safety is a global public health concern and one of the leading causes of death for children over the age of 5 years [1, 2]. Child pedestrians are particularly vulnerable road users due to their limited physical, cognitive-perceptual, and social development [3]. According to the UK Department for Transport, only in 2019 in England, 4700 child pedestrians under the age of 15 sustained traffic-related injuries, out of which 1200 were killed or seriously injured. Ensuring children’s safety on roads is a major public health priority as it can prevent various adverse physical, mental, and social consequences and can promote walking among children, increasing childhood physical activity.

Adopting active modes of travel, including walking and cycling, from early ages has positive impacts on both personal and planetary health. However, the number of children walking has in general declined over the last decades. For example, based on the National Travel Survey statistics, in England in 1998/2000, an estimated 49% of children under 16 walked to school while in 2017 this rate fell to 43%. Research suggests that this reduction is partly due to traffic safety concerns [4, 5]. To help address this issue and reverse this trend, it is essential that policies are geared towards making roads safer and enjoyable for children.

Child Pedestrian injury is a multi-faceted problem with many contributory factors such as driver characteristics [6], vehicle features [7, 8], road configuration [9, 10], environmental [11], built-environment [12], and socio-economic and ethnic features [13,14,15,16,17,18]. In addition to these factors, there are geographical and temporal variations influenced by place- and time-specific factors and interventions such as reduced vehicles circulations around schools, implementation of 20 mph speed limit zones, school walking bus, School Street schemes, among the others.

Previous research has mainly focused on identifying high-crash areas and linking the risk to explanatory variables to quantify the effect of several risk factors through a combination of geographic information systems and statistical models [19,20,21,22,23,24]. Due to the nature of road crash data, there might be spatial and temporal dependencies between observations; and therefore, statistical models need to accommodate these dependencies [25]. However, considering spatial-temporal dependencies [26] as well as child safety [15] is relatively limited in the crash literature, especially in England. In this study we investigate the spatiotemporal patterns of child pedestrian safety at Lower Tier Local Authorities (LTLA) level in England from 2011 to 2020. In England, the LTLA is a subdivision of administrative boundaries that include local authority districts, unitary authorities, metropolitan districts, and London boroughs. The main aims are to (i) explain the association between child pedestrian crash frequencies and LTLA level characteristics, including a host of deprivation, transport, and built-environment variables, (ii) identify LTLAs with particularly high crash frequencies for child pedestrians across the study period, (iii) evaluate the persistence of spatial patterns of child crashes over time, and (iv) pinpoint local time trends for each LTLA.

We conducted our analysis within a Bayesian framework as the Bayesian approach is particularly appealing in handling spatio-temporal dependencies in the data, readily allowing for borrowing of strength across space and time that leads to more reliable estimates. The Bayesian approach, unlike the frequentist approach that obtains point estimates of parameters of interest, regards the unknown parameters as random variables and obtains posterior densities for all model parameters, addressing uncertainties more fully. For a discussion on the advantages of the Bayesian methods in the context of road safety analysis, see [27]. The posterior density is the result of data (i.e., the likelihood function) integrated with a prior probability distribution that represents our prior belief regarding a parameter of interest. When there is no prior knowledge of a parameter, non-informative priors can be used, allowing the data to decide the form of the posterior entirely. For estimating the model parameters, the Bayesian approach employs Markov chain Monte Carlo simulations [28].



The outcome of interest was the annual counts of crashes involving child pedestrians at local authority level in England. Crash data were obtained from the Department for Transport, which collects information on crashes that occurred on public roads, reported to the police, and recorded on STATS19 forms. As the focus of our study was to model child pedestrian crash frequencies, information relating to the location of crashes, road user types (e.g., pedestrian, cyclist, driver), and age of individuals involved in crashes were extracted from the STATS19 databases. We only included crashes between a motorised vehicle and a child pedestrian, who was less than 16 years old. This yielded 50,993 crashes over the period 2011–2020. Using the geographic coordinates of the crashes, we obtained yearly crash counts at LTLA level. We removed Isles of Scilly due to the sparsity of outcome data and the City of London because many explanatory variables for the latter were missing. Isle of Scilly and the City of London had one and ten cases of child pedestrian crashes over the study period, respectively. This resulted in 315 LTLAs with a mean crash of 161, standard deviation of 169.83, minimum of 6, and a maximum of 1808 over the study period. The geocoded location of crashes is mapped in Additional File 1.

To model child pedestrian crash frequencies, we considered several sociodemographic, transport, and built-environment features, based on literature, domain expertise, and data availability. As sociodemographic covariates, we considered percent of child (0–15 year) population [29], number of licensed vehicles per capita [30], percent of claimants, percent population who are White, and job density [31]. The number of licensed vehicles in 2020 per capita was calculated by dividing the total number of registered vehicles by local authority population. The percent of claimants for 2011–2020 was defined as the proportion of residents aged 16–64 claiming some form of unemployment-related benefit. From the annual population survey [31], we obtained yearly unemployment rate, and yearly percent of population in employment who are manager, directors, and senior officials for 2011–2020.

As transport-related variables indicating travel behaviour and exposure at LTLA level, we included the percent of adults who walk or cycle at least three times per week, the percent of adults who walk or cycle at least five times per week [32], and road density. To obtain road density, the latest road network data sourced from Ordnance Survey Meridian [33] was overlaid with LTLA boundaries. We then calculated road density by dividing the total length (in km) of A roads and B roads by the total land area (in km2) of each LTLA. The A roads are major roads linking town and cities and the B roads are distributor roads with lower traffic density than A roads and are intended to connect different areas [34].

Built-environment variables included the number of schools, the number of bus stops, and the number of business establishments (as a proxy of activity levels) in each LTLA. School information was provided by the Department for Education “get information about schools” register (downloaded in December 2021 from, bus stop locations were retrieved from Point of Interest Ordnance Survey data, and the number of business establishments were obtained from the Office for National Statistics. Note that some variables were available at a yearly basis. For other variables we considered the latest and/or the most relevant available data. This is because the data availability and LTLA boundaries have changed over the study period. The descriptive statistics of the explanatory variables and data sources and spatial distribution of a number of covariates are reported in Additional File 1.

Statistical analysis

We used a Bayesian space-time Poisson lognormal model to evaluate associations between various relevant contributory factors and annual child pedestrian crashes, while accommodating dependencies between adjacent LTLAs and years. We initially attempted to account for spatio-temporal heterogeneity in the data through available covariates. However, due to unobserved/unmeasured factors that affect the outcome, which can be themselves spatially and temporally correlated, there remains residual autocorrelation. We dealt with the residual autocorrelation through specifying random effects that can account for the spatial and temporal dependencies and act as surrogate for unobserved/unmeasured covariates, affecting the outcome. This is achieved by including spatial, temporal, and spatio-temporal weight matrices that specify the neighbourhood structure, which allows us to exploit information sharing between neighbouring LTLAs and years.

After accounting for the covariates, we incorporated residual spatial variability through spatially structured random effects represented by a conditional autoregressive (CAR) prior [35]. The residual spatial term accounts for the dependencies between neighbouring LTLAs, which are conceptualised as LTLAs that share a common border. To account for the temporal dependency, we included temporally structured random effects using a random walk of order 1 (RW1), which captures the national temporal trend. Additionally, we included a space-time interaction term modelled as independent random walk for each LTLA [36]. The interaction term adds additional flexibility to the model and allows capturing local temporal deviations from the national (overall) time trend. Therefore, the temporal patterns in each LTLA are assumed to be temporally smooth but independent across space.

Model specification

In our model specification, the number of observed crashes involving child pedestrians yit in LTLA i (i = 1, …,315) and year t (t = 1,..,10) follows a Poisson distribution with the mean μit. We then decompose the log (μit) as

$$\log \left({\mu}_{it}\right)=\alpha +{X}_{it}\beta +{\phi}_i+{\xi}_t+{\delta}_{it}$$

where exp(α) represents the overall expected child pedestrian crash frequency (i.e., an intercept term), ϕi are spatially structured random effects that capture the main spatial patterns and ξt are temporally structured random effects to describe the global temporal pattern of crashes (i.e., national time trend) in England. The term δit represents a space-time interaction term and allows each LTLA to have a temporal trend deviating from the national one. The term Xit represents the covariate matrix of the ith LTLA in time t and β are their respective regression coefficients. The spatially structured random effects ϕi are assigned a CAR prior as in

$${\phi}_i\mid {\phi}_{-i}\sim Normal\left(\frac{1}{n_i}\sum_{j\ne i}{\phi}_j,\frac{\sigma_{\phi}^2}{n_i}\right)$$

where ϕi is the set of ϕi except for the ith LTLA, ni is the number of adjacent neighbours sharing a boundary to area i, and \({\sigma}_{\phi}^2\) is the variance of spatially structured random effects. The temporally structured random effects ξt are modelled using a random walk of order 1 (RW1) and is described as

$${\xi}_t\mid {\xi}_{t-1}\sim Normal\left({\xi}^{t-1},{\sigma}_{\xi}^2\right)$$

where \({\sigma}_{\xi}^2\) is the variance of temporal effects. In addition to the main temporal effects, we were interested in identifying localised departures. This was achieved by including a space-time interaction structure δit defined as

$${\delta}_{it}\sim Normal\left({\delta}^{t-1},{\sigma}_{\delta}^2\right)$$

As for priors, we assigned a Gamma (0.5, 0.005) to the inverse of variances \({\sigma}_{\phi}^2\), \({\sigma}_{\xi}^2\), and \({\sigma}_{\delta}^2\), a flat prior uniform(−∞,+∞) to α, and a Normal (0,1000) to the regression coefficients β. All explanatory variables were centred to improve the convergence of the model. For better interpretability of the explanatory variables, we obtained the magnitude of the effects of various covariates on child pedestrian safety in terms marginal effects [37] given by

$$\frac{\partial E\left(y|x\right)}{\partial {x}_p}=E\left(y|x\right){\beta}_p=\frac{1}{N}\sum_{i=1}^N\exp \left(\alpha +{X}_{it}\beta +{\phi}_i+{\xi}_t+{\delta}_{it}\right){\beta}_{i,p}$$

where subscript p refers to explanatory variables, here p = 6, and N is the number of observations.

Inferences were performed through Markov Chain Monte Carlo (MCMC) simulations in the NIMBLE Package in R [38]. We checked the convergence of the parameters using the Gelman-Rubin statistic [39] and visually using trace plots. In total, 20,000 post burn-in samples were obtained from the posterior distribution of the model parameters. In addition to the above-described model, we fitted other competing models with different specifications for space, time, and space-time effects, and compared the model fit using the Watanabe–Akaike information criterion (WAIC). However, these models did not improve the fit. Details of the competing models and their WAIC are reported in additional file 2.

We report the posterior summary of the magnitude of the effects of various covariates on child pedestrian safety in terms of marginal effects [40]. Marginal effects provide a more straightforward interpretation of the effect of covariates on safety, revealing the change in expected (mean) child pedestrian crashes following one unit change in each covariate. To check the goodness of fit of the model, we conducted posterior predictive checks and estimated Bayesian p-values [41], which is based on quantifying the discrepancies between predicted data, using the proposed model, and the observed data.

Identifying high-crash areas (hotspots), spatial distribution of residuals and area-specific time trends

Using the expected crash frequency, we identified high-crash areas, where safety improvement programmes are most warranted. We classified LTLAs based on the posterior probability of the exponential of the spatial residuals in each LTLA being above one. If a probability was larger than 0.9 for an LTLA, it was classified as an LTLA with excess child pedestrian crash frequency (after accounting for explanatory variables that are in the model). Such areas are where other unknown/unmeasured risk factors (other than those in the model) have a negative impact on child pedestrian safety. Therefore, further investigation is needed to identify the reasons behind this, which in turn helps improve safety in those areas.

The inclusion of the space-time interaction term allowed each LTLA to have its own specific temporal trend, which is composed of the sum of the national temporal trend and the space-time interaction term. Similar to the approach adopted by Boulieri et al. [42], we report the probability that the estimated incidence of child pedestrian crashes in an LTLA represents an increase compared to the national one.


LTLA factors affecting child pedestrian safety

Table 1 reports the estimated model parameters, including regression coefficients and parameters relating to our space-time interaction specification. We retained only statistically important covariates, considering those that were not highly correlated with each other, and that improved model fit. Child population, unemployment-related claimants, road density, and the number of schools were found to be positively associated with child pedestrian crash frequencies. However, licensed vehicles per capita and levels of walking and cycling were negatively associated with child pedestrian crash frequencies. With respect to model performance, we estimated Bayesian p-values for checking model adequacy. These being satisfactory are reported in the Additional File 3.

Table 1 Posterior summary of regression parameters

To interpret the regression coefficients, Table 2 reports the magnitude of the impact of explanatory variables on child pedestrian crash frequencies in terms of marginal effects, indicating the magnitude of change in child crash frequency due to one unit change in an explanatory variable.

Table 2 Posterior summary of marginal effects

In relation to deprivation, Fig. 1 displays that number of child pedestrian crashes (on log scale) in 2011 and 2019 in relation to the percent of population claiming some sort of unemployment-related benefit. Since 2020 included lock down periods and major shift in travel behaviour due to the COVID-19 pandemic [43, 44], we restricted our comparison to 2011 and 2019 as the two extremes of the study period. As shown in Fig. 1, the expected number of crashes was positively associated with the deprivation level. The slope of the 2019 line, however, is slightly larger than that of 2011, suggesting that as the percentage of unemployment-related claimants increases, crash frequencies increase at a higher rate in 2019 than in 2011. This implies socio-economic inequalities in child pedestrian safety increased in 2019 compared to 2011, highlighting the importance of the need for addressing inequity issues in this context.

Fig. 1
figure 1

Expected crash frequency in relation to deprivation

High-crash locations

Figure 2 displays the spatial distribution of yearly expected child pedestrian crashes over the study period (the darker the colour, the higher the expected value). There is a relatively considerable spatial variation in expected child pedestrian crashes across England which remains visible across the years while we can also see the temporal evolution of child pedestrian safety in England over the last decade. Although we expect that urban areas have higher expected crash values, Fig. 2 can be used by local authorities to prioritise safety interventions and to inform resource allocation across England. For example, in 2019, Birmingham had the highest expected child pedestrian crash, followed by Leeds, Bradford, Liverpool, and Croydon, which all are characterized by low socioeconomic status.

Fig. 2
figure 2

Spatial distribution of expected child pedestrian crash frequencies from 2011 to 2020

Overall spatial and temporal effects

Figure 3a shows the map of exceedance probabilities of spatial residuals being greater than 1. This allows us to identify LTLAs with excess child pedestrian crash (shown in darker colour where the probability of exceedance is > 90%), after adjusting for the effect of the covariates. We found that 36% of the LTLAs (114 LTLAs) experienced excess crash from 2011 to 2020. These are mainly located in urban areas, especially in Northern England: Yorkshire and the Humber regions.

Fig. 3
figure 3

a Mapping posterior probability of spatial residuals being larger than 1. The map of Greater London is enlarged for better visualization. b Posterior median and 95% credible intervals of the overall (national) temporal trend

Figure 3b shows the posterior median, including the 95% uncertainty band of the temporal trend, over the study period, representing the average national time trend. We observed a decreasing time trend, with year 2020 showing a much steeper decline compared to the other years, perhaps reflecting the effect of the Covid19 pandemic and its associated lockdown and work from home policies in England. Note that such policies resulted in reduced exposure (traffic volume, and walking and cycling) in general [45]. Between the years 2013 and 2017, we observed a relatively moderate but consistent decreasing trend. Figure 4 confirms a non-linear behaviour of temporal patterns over time.

Fig. 4
figure 4

Probability of LTLA- specific time trend exceeding the national trend over the study period

LTLA-specific time trends through space-time interaction

By specifying a space-time interaction term as an independent random walk for each LTLA, we were able to capture local time trends. Such local trends are due to the fact that the effects of some unknown highly localised variables vary smoothly over time while operating independently with respect to their locations. In fact, previous research in the field of road safety indicates that the effect of contributory factors may vary over time [46]. In the presence of missing localised variables (e.g., climate), the interaction term can act as a surrogate measure for these unmeasured/unknown variables. This allows us to capture their effects to some extent, thereby addressing unobserved heterogeneity more fully. The time trend for all LTLAs exhibited downward trend (similar to the national trend shown in Fig. 3b) with different degrees of deviation from the national trend (see Additional File 4). The decline in crash incidence was slower in some local authorities when compared to the national trend.

Figure 4 displays the map of the probability that, after accounting for the covariates, the incidence of child pedestrian crashes is higher than the national one in each LTLA in each year. The incidence of child pedestrian crash in North Lincolnshire was below the average national incidence between 2011 and 2013, but it exceeded the national value (with probability > 80%) after 2018. This trend suggests that there are specific risk factors in North Lincolnshire that contribute to the deterioration of road safety for children, which requires further in-depth investigations. In contrast, the crash incidence in the local authority of Bournemouth, Christchurch and Poole in 2011–2013 resulted higher than the average national, but child pedestrian safety improved over the recent years such that after 2017 the incidence of crash became less than the average national child pedestrian crash incidence. For a detailed visualisation, the time trend for the highlighted local authorities (shown blue and red colours) in comparison to the national trend is illustrated in the Additional File 4.


Association between LTLAs characteristics and child pedestrian safety

Identifying contributory factors affecting zonal-level child pedestrian safety can provide useful insights toward designing and implementing effective large-scale countermeasures. For example, based on our findings, increased levels of walking and cycling helps increase road safety for children. This is an interesting finding that, in accordance with previous research (see; e.g., Stoker et al., [47] and Jacobsen et al., [48]) indicates that the higher the prevalence of walking and cycling in LTLAs, the safer the road network for child pedestrians. Also, this could be attributed to motorists adjusting their driving behaviours (e.g., by lowering driving speeds) in the presence of increased numbers of pedestrians and cyclists [49,50,51].

Specifically, based on marginal effects reported in the section of results, we found that one unit increase in road density can lead to an additional 11.01 child pedestrian crashes annually. Road density can act as a proxy exposure measure for motorised traffic, which is known to have a deteriorating effect on pedestrian safety [52]. Our finding regarding road density is consistent with the results of previous studies [53]. One unit increase in child population, on average, resulted in 1.8 additional child pedestrian crashes per year. For every 10 additional schools, expected child pedestrian crash frequencies increased by 1.27 per year. One potential explanation is that an increase in child population and the number of schools leads to an increase in exposure [26, 49, 54]. In terms of deprivation, one unit (here, 1 %) increase in unemployment-related claimants resulted in 0.48 additional child pedestrian crashes per year. In contrast, one unit increase in the number of licensed vehicles per capita decreased child pedestrian crash frequencies by 8.98 crashes per year. One explanation for this finding could be that children in areas with higher number of licensed vehicles are more likely to travel by car rather than walking, reducing their exposure. Also, the latter two variables often relate to deprivation and previous studies have also found similar results [15, 26]. Finally, 1 % increase of adults who walk or cycle at least three times per week decreased expected child pedestrian crash frequencies by 0.18 crashes per year.

One important factor is deprivation which has a negative impact on road safety. The socio-economic disparities in child pedestrian crashes might be partially driven by exposure disparities as children in deprived areas are more likely to walk to school [55, 56] Therefore, road safety policies need to target more deprived areas through safety improvement programmes such as reducing traffic volume and speed, designing walking-friendly infrastructures, education, and training programmes. In addition, whilst child pedestrian crash frequencies have declined over the last decade, its association with deprivation over time has not changed substantially. If road safety interventions successfully target deprived areas, the association between deprivation and child crash frequency would weaken in the future. With increasing interest in policies to encourage children to walk, efforts to improve child pedestrian safety is successful only when a system approach is adopted, emphasising on data-informed engineering interventions in conjunction with interventions that address deprivation.

Spatio-temporal variations in child pedestrian safety

Since 2011, child pedestrian crashes have decreased by more than 50% in England. However, some local authorities still struggle to improve road safety conditions for children. Overall, the southern part of England has experienced higher levels of improvement in child pedestrian safety over the last decade compared to the northern regions. Local authorities of Birmingham, Leeds, and Bradford had consistently the highest expected child pedestrian crash frequency throughout England over the study period. In 2011, there was a difference of 218.3 [95% CrI 197.1–240.9] crashes between the LTLA with the lowest (Rutland) and the highest child pedestrian crash frequencies (Birmingham). Although road safety has in general improved over the study period, in 2020, there was a gap of 112.22 [98.9–125.4] expected annual crashes between LTLAs (Rutland and Birmingham) with the highest and the lowest expected crash frequencies.

Estimating the LTLA-specific time trends suggests the presence of localised risk factors or may reflect the impact of local interventions and policies, which requires further in-depth investigations. This can be particularly important from a public health perspective and for implementing cost-effective safety interventions. We noticed a major reduction in expected child pedestrian crash frequencies in most local authorities in 2020, which is expected due to the recent pandemic (see, for example, the work by Katrakazas et al. [57] for a discussion on the effect of the Covid-19 pandemic on road safety). However, the expected child pedestrian crash frequency was less affected in certain regions such as Cornwall and Northumberland.

From a policy insight perspective and with the aim of improving child pedestrian safety across England, our results can be used to prioritise safety interventions and to inform resource allocation across England.. This can be achieved based on the identification of high-crash local authorities (so called hotspots), understanding factors affecting child pedestrian safety, and tracking how safety conditions, in terms of child pedestrian safety, at different local authorities have evolved over the study period.

Strengths and limitations of the study

We accommodated spatial dependencies that could potentially account for similarities in unobservable factors and travel patterns (e.g., children mobility and traffic volume) in neighbouring local authorities. Our specification of space-time interaction allowed each local authority to have its own temporal pattern. Consequently, this specification not only led to more reliable statistical inferences but also allowed us to provide further insights with the same set of data.

A limitation of our study is that we could not use some other potentially useful variables such as traffic volume and land use characteristics in our analysis as their LTLA boundaries did not match the boundaries associated with our outcome of interest. However, we showed that our model performs very well in replicating the observed data so this would not cause any major issue in this study. This is partly because we indirectly accounted for these omitted variables in the model. For example, we included road density as a proxy measure for traffic volume. Due to lack of data availability, we removed Isles of Scilly and City of London. However, since the number of observed crashes in these two local authorities was very small, we do not expect that their removal substantially affects the results.


To our knowledge, this study is the first to explore spatial-temporal patterns of child pedestrian crashes at local authority level in England over a long-time span,10 years (2011–2020). We used a Bayesian spatio-temporal model where we included relevant covariates, accounted for spatial and temporal dependencies in the data as well as allowing for a spatio-temporal interaction. This enabled us to (i) identify statistically important area-level variables that can explain child pedestrian safety, (ii) reveal spatial patterns and national trend in child pedestrian crashes across England over the last decade, and (iii) understand how road safety conditions evolved in each local authority over the study period. The results indicate that child pedestrian crashes have been gradually declining in England over the last decade. Some local authorities (mainly in urban areas of northern England) exhibited higher child crash frequencies than national average over the study period. More deprived local authorities have been experiencing a higher number of child pedestrian crashes and there is no evidence suggesting that socioeconomic-related inequality gap has narrowed from 2011. Efforts to improve child pedestrian safety would be more successful if emphasise is given to areas where safety improvements are most warranted and to evidence-based policy making in conjunction with interventions that can address social inequalities.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files (Additional file 5).



lower tier local authority


conditional autoregressive


random walk of order 1


Markov Chain Monte Carlo


Watanabe–Akaike information criterion


  1. World Health O. Ten strategies for keeping children safe on the road. Geneva: World Health Organization; 2015. Contract No.: WHO/NMH/NVI/15.3

    Google Scholar 

  2. Hyder AA, Peden M. Inequality and road-traffic injuries: call for action. Lancet. 2003;362(9401):2034–5.

    Article  Google Scholar 

  3. Schieber R, Vegega M. Education versus environmental countermeasures. Inj Prev. 2002;8(1):10–1.

    Article  CAS  Google Scholar 

  4. Francis J, Martin K, Wood L, Foster S. ‘I’ll be driving you to school for the rest of your life’: a qualitative study of parents’ fear of stranger danger. J Environ Psychol. 2017;53:112–20.

    Article  Google Scholar 

  5. Nikitas A, Wang JY, Knamiller C. Exploring parental perceptions about school travel and walking school buses: a thematic analysis approach. Transp Res A Policy Pract. 2019;124:468–87.

    Article  Google Scholar 

  6. Fridman L, Pitt T, Rothman L, Howard A, Hagel B. Driver and road characteristics associated with child pedestrian injuries. Accid Anal Prev. 2019;131:248–53.

    Article  Google Scholar 

  7. Paulozzi LJ. United States pedestrian fatality rates by vehicle type. Inj Prev. 2005;11(4):232.

    Article  CAS  Google Scholar 

  8. Henary BY, Crandall J, Bhalla K, Mock CN, Roudsari BS. Child and adult pedestrian impact: the influence of vehicle type on injury severity. In: Annual Proceedings/Association for the Advancement of Automotive Medicine: Association for the Advancement of Automotive Medicine; 2003.

    Google Scholar 

  9. Lee G, Park Y, Kim J, Cho G-H. Association between intersection characteristics and perceived crash risk among school-aged children. Accid Anal Prev. 2016;97:111–21.

    Article  Google Scholar 

  10. Yu C-Y. How differences in roadways affect school travel safety. J Am Plan Assoc. 2015;81(3):203–20.

    Article  CAS  Google Scholar 

  11. Amoh-Gyimah R, Aidoo EN, Akaateba MA, Appiah SK. The effect of natural and built environmental characteristics on pedestrian-vehicle crash severity in Ghana. Int J Inj Control Saf Promot. 2017;24(4):459–68.

    Article  Google Scholar 

  12. Merlin LA, Guerra E, Dumbaugh E. Crash risk, crash exposure, and the built environment: a conceptual review. Accid Anal Prev. 2020;134:105244.

    Article  Google Scholar 

  13. Christie N, Ward H, Kimberlee R, Towner E, Sleney J. Understanding high traffic injury risks for children in low socioeconomic areas: a qualitative study of parents’ views. Inj Prev. 2007;13(6):394–7.

    Article  CAS  Google Scholar 

  14. Green J, Muir H, Maher M. Child pedestrian casualties and deprivation. Accid Anal Prev. 2011;43(3):714–23.

    Article  Google Scholar 

  15. Graham DJ, Stephens DA. Decomposing the impact of deprivation on child pedestrian casualties in England. Accid Anal Prev. 2008;40(4):1351–64.

    Article  Google Scholar 

  16. Hwang J, Joh K, Woo A. Social inequalities in child pedestrian traffic injuries: differences in neighborhood built environments near schools in Austin, TX, USA. J Transp Health. 2017;6:40–9.

    Article  Google Scholar 

  17. Chakravarthy B, Anderson CL, Ludlow J, Lotfipour S, Vaca FE. A geographic analysis of collisions involving child pedestrians in a large Southern California County. Traffic Inj Prev. 2012;13(2):193–8.

    Article  Google Scholar 

  18. Hamann C, Peek-Asa C, Butcher B. Racial disparities in pedestrian-related injury hospitalizations in the United States. BMC Public Health. 2020;20(1):1459.

    Article  Google Scholar 

  19. Boulieri A, Liverani S, de Hoogh K, Blangiardo M. A space-;time multivariate Bayesian model to analyse road traffic accidents by severity. J R Stat Soc Series A (Statistics in Society). 2017;180(1):119–39.

    Article  Google Scholar 

  20. Jones AP, Haynes R, Kennedy V, Harvey IM, Jewell T, Lea D. Geographical variations in mortality and morbidity from road traffic accidents in England and Wales. Health Place. 2008;14(3):519–35.

    Article  CAS  Google Scholar 

  21. Montella A, Marzano V, Mauriello F, Vitillo R, Fasanelli R, Pernetti M, et al. Development of macro-level safety performance functions in the City of Naples. Sustainability. 2019;11(7):1871.

    Article  Google Scholar 

  22. Abdel-Aty M, Lee J, Siddiqui C, Choi K. Geographical unit based analysis in the context of transportation safety planning. Transp Res A Policy Pract. 2013;49:62–75.

    Article  Google Scholar 

  23. Liu C, Sharma A. Using the multivariate spatio-temporal Bayesian model to analyze traffic crashes by severity. Analytic Methods Accident Res. 2018;17:14–31.

    Article  Google Scholar 

  24. Quistberg DA, Hessel P, Rodriguez DA, Sarmiento OL, Bilal U, Caiaffa WT, et al. Urban landscape and street-design factors associated with road-traffic mortality in Latin America between 2010 and 2016 (SALURBAL): an ecological study. Lancet Planet Health. 2022;6(2):e122–e31.

    Article  Google Scholar 

  25. Quddus MA. Modelling area-wide count outcomes with spatial correlation and heterogeneity: an analysis of London crash data. Accid Anal Prev. 2008;40(4):1486–97.

    Article  Google Scholar 

  26. Graham DJ, McCoy EJ, Stephens DA. Quantifying the effect of area deprivation on child pedestrian casualties by using longitudinal mixed models to adjust for confounding, interference and spatial dependence. J Royal Stat Society Series A-Stat Soc. 2013;176:931–50.

    Article  Google Scholar 

  27. Heydari S, Fu L, Thakali L, Joseph L. Benchmarking regions using a heteroskedastic grouped random parameters model with heterogeneity in mean and variance: applications to grade crossing safety analysis. Analytic Methods Accident Res. 2018;19:33–48.

    Article  Google Scholar 

  28. Gelman A, Carlin JB, Stern HS, Rubin DB, Brix A. Bayesian data analysis. Blackwell Publishers. 2005. p. 251–2.

  29. Office for National Statistics. Estimates of the population for the UK, England and Wales, Scotland and Northern Ireland [Available from:

  30. Department for Transport Driver and Vehicle Licensing Agency. Vehicle licensing statistics data tables [Available from:

  31. Office for National Statistics. official census and labour market statistics [Available from:

  32. Department for Transport. Walking and cycling statistics [Available from:

  33. Ordnance Survey. In: Survey O, editor. Meridian 2; 2016.

    Google Scholar 

  34. Department for Transport. Road lengths in Great Britain statistics: Notes and definitions. 2022.

    Google Scholar 

  35. Besag J. Spatial interaction and the statistical analysis of lattice systems. J R Stat Soc Ser B Methodol. 1974;36(2):192–236.

    Google Scholar 

  36. Knorr HL. Bayesian modelling of inseparable space-time variation in disease risk. Stat Med. 2000;19(17–18):2555–67.

    Article  Google Scholar 

  37. Winkelmann R. Econometric analysis of count data. Springer Science & Business Media; 2008.

  38. de Valpine P, Turek D, Paciorek CJ, Anderson-Bergman C, Lang DT, Bodik R. Programming with models: writing statistical algorithms for general model structures with NIMBLE. J Comput Graph Stat. 2017;26(2):403–13.

    Article  Google Scholar 

  39. Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7(4):457–72, 16.

    Article  Google Scholar 

  40. Washington S, Karlaftis M, Mannering F, Anastasopoulos P. Statistical and econometric methods for transportation data analysis. Chapman and Hall/CRC; 2020.

  41. Meng X-L. Posterior Predictive p-Values. Ann Stat. 1994;22(3):1142–60.

    Article  Google Scholar 

  42. Boulieri A, Blangiardo M. Spatiotemporal model to estimate life expectancy and to detect unusual trends at the local authority level in England. BMJ Open. 2020;10(11):e036855.

    Article  Google Scholar 

  43. Noland RB. Mobility and the effective reproduction rate of COVID-19. J Transp Health. 2021;20:101016.

    Article  Google Scholar 

  44. Heydari S, Konstantinoudis G, Behsoodi AW. Effect of the COVID-19 pandemic on bike-sharing demand and hire time: evidence from Santander cycles in London. PLoS One. 2021;16(12):e0260969.

    Article  CAS  Google Scholar 

  45. Hadjidemetriou GM, Sasidharan M, Kouyialis G, Parlikad AK. The impact of government measures and human mobility trend on COVID-19 related deaths in the UK. Transp Res Interdiscip Perspect. 2020;6:100167.

    Google Scholar 

  46. Mannering F. Temporal instability and the analysis of highway accident data. Anal Methods Accid Res. 2018;17:1–13.

    Google Scholar 

  47. Stoker P, Garfinkel-Castro A, Khayesi M, Odero W, Mwangi MN, Peden M, et al. Pedestrian Safety and the Built Environment:A Review of the Risk Factors. J Plan Lit. 2015;30(4):377–92.

    Article  Google Scholar 

  48. Jacobsen PL. Safety in numbers: more walkers and bicyclists, safer walking and bicycling. Inj Prev. 2003;9(3):205–9.

    Article  CAS  Google Scholar 

  49. Heydari S, Miranda-Moreno L, Hickford AJ. On the causal effect of proximity to school on pedestrian safety at signalized intersections: a heterogeneous endogenous econometric model. Anal Methods Accid Res. 2020;26:100115.

    Google Scholar 

  50. Jacobsen PL, Ragland DR, Komanoff C. Safety in numbers for walkers and bicyclists: exploring the mechanisms. Inj Prev. 2015;21(4):217–20.

    Article  Google Scholar 

  51. Klieger J, Savage I. Motor-vehicle drivers’ behavioral response to increased bicycle traffic. J Saf Res. 2020;74:97–102.

    Article  Google Scholar 

  52. Heydari S, Fu L, Miranda-Moreno LF, Jopseph L. Using a flexible multivariate latent class approach to model correlated outcomes: a joint analysis of pedestrian and cyclist injuries. Anal Methods Accid Res. 2017;13:16–27.

    Google Scholar 

  53. Xie K, Ozbay K, Yang H. A multivariate spatial approach to model crash counts by injury severity. Accid Anal Prev. 2019;122:189–98.

    Article  Google Scholar 

  54. Bhat CR, Astroza S, Lavieri PS. A new spatial and flexible multivariate random-coefficients model for the analysis of pedestrian injury counts by severity level. Anal Methods Accid Res. 2017;16:1–22.

    Google Scholar 

  55. Brophy S, Cooksey R, Lyons RA, Thomas NE, Rodgers SE, Gravenor MB. Parental factors associated with walking to school and participation in organised activities at age 5: analysis of the millennium cohort study. BMC Public Health. 2011;11(1):14.

    Article  Google Scholar 

  56. Sonkin B, Edwards P, Roberts I, Green J. Walking, cycling and transport safety: an analysis of child road deaths. J R Soc Med. 2006;99(8):402–5.

    Article  Google Scholar 

  57. Katrakazas C, Michelaraki E, Sekadakis M, Yannis G. A descriptive analysis of the effect of the COVID-19 pandemic on driving behavior and road safety. Transp Res Interdiscip Perspect. 2020;7:100186.

    Google Scholar 

Download references


Infrastructure support for the Department of Epidemiology and Biostatistics was provided by the NIHR Imperial Biomedical Research Centre (BRC).


This work was supported by an MRC Early Career Research Fellowship awarded to NS via the MRC Centre for Environment and Health (Award Ref. MR/T502613/1).

Author information

Authors and Affiliations



N.S: Conceptualization, Investigation, Methodology, Formal analysis, Funding acquisition, Writing – original draft; S.H: Conceptualization, Investigation, Methodology, Writing – editing; M.B.: Conceptualization, Methodology, Writing – editing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Niloofar Shoari.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: 

information on data. Figure 1. Geocoded location of child pedestrian crashes across England from 2011 to 2020. Table 1. Summary statistics of the explanatory variables. Figure 2. Spatial distribution of road density in England (1/km). Figure 3. Spatial distribution of Number of schools in England. Figure 4. Spatial distribution of the number of vehicles per capita in England. Figure 5. A Spatial distribution of the percent of adults who walk/cycle at least three times per week in England. Figure 6. Spatial distribution of the percent of child population over years in England. Figure 7. Spatial distribution of the percent of population who claim unemployment-related benefit over years in England.

Additional file 2: 

Competing models and model comparison. Table 2. Specification of competing models. Table 3. The model fit and the estimated parameters with 95% credible intervals for competing models.

Additional file 3: 

Checking goodness-of-fit. Table 4. Average Bayesian p-values.

Additional file 4: 

LTLA-specific time trends. Figure 8. LTLA-specific time trends. Figure 9. Time trend of local authorities where the crash incidence exceeds the average national value over time. Figure 10. Time trend of local authorities where the crash incidence goes below the average national value over time.

Additional file 5.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Shoari, N., Heydari, S. & Blangiardo, M. A decade of child pedestrian safety in England: a Bayesian spatio-temporal analysis. BMC Public Health 23, 215 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Bayesian spatio-temporal
  • Child pedestrian safety
  • Road crashes
  • Hierarchical models