Common before-after accident study on a road site: a low-informative Bayesian method

<jats:title>Abstract</jats:title><jats:sec><jats:title>Purpose</jats:title><jats:p>This note aims at providing a Bayesian methodological basis for routine before-after accident studies, often applied to a single road site, and in conditions of limited resources in terms of time and expertise.</jats:p></jats:sec><jats:sec><jats:title>Methods</jats:title><jats:p>A low-informative Bayesian method is proposed for before-after accident studies using a comparison site or group of sites. As compared to conventional statistics, the Bayesian approach is less subject to misuse and misinterpretation by practitioners. The low-informative framework seems appropriate in situations of limited expertise. The proposed approach gives the possibility of correcting for regression to the mean. Examples illustrate the application of this method.</jats:p></jats:sec><jats:sec><jats:title>Results and conclusions</jats:title><jats:p>It is shown that a relatively simple method, based on the Jeffreys’s rule prior considered as a “reasonable standard”, can be implemented without major difficulties. Posterior distributions are proper. The numerical calculation of posterior probabilities can be done without using Monte-Carlo simulations nor specialised software tools.</jats:p></jats:sec>

section_content: It is common that road sites are modified in order to achieve improvements from various points of view (traffic conditions, better integration of various uses and users of the road and public space, reduction of noise and air pollution, traffic safety, etc. ). A few years after a site has been modified, local engineers generally have to study the effects of this road change, regarding various aspects including road safety. Thus, a retrospective before-after accident study is often needed. 

In such routine situations, resources are limited in terms of time and expertise, and the risk of misuse of conventional statistical methods is increased. Even among people who are more experienced in statistics, like researchers, erroneous uses of conventional methods are common: misuse of tests of significance, erroneous understanding of p-values, misinterpretation of confidence intervals (as pointed out by many authors [15, 18, 19, 24, 27, 32] ; see also [5, 11, 28] ). For example, the p-value is often erroneously regarded as the probability that the null hypothesis is true, and the 95% confidence interval obtained is wrongly assumed to contain the true parameter with a 95% chance. The Bayesian approach to statistics is more in accordance with the expectations and intuitions of non-specialists. In particular, the posterior distribution can be legitimately used to give the probabilities that the parameter of interest is contained in various regions of the parameter space (a 95% credible interval, for example), or exceeds a particular value, given the data observed and prior knowledge. Some authors consider that teaching Bayesian statistics is easier than teaching frequentist statistics [10, 31]. Nevertheless, aids to practitioners are necessary to implement Bayesian methods, since the calculations in these approaches are sometimes complex. 

In this paper, we will not deal with studies based on large samples of sites and using multivariate modelling, for which Bayesian approaches were proposed in the recent period [4, 30, 35, 37]. Bayesian methods adapted to meta-analyses or to overviews of several studies (see, for example, [12] ) will not be considered here. We will focus on methods applicable to a single site and transferable to engineers for common practice. 

In the case we deal with here (routine evaluation, single site), the methods currently used and recommended are conventional statistical methods (see, for example, [23] ), even though they sometimes make use of empirical Bayes estimates of the expected accident number on the treated site in order to cope with 'regression to the mean' bias. The principle of a 'full' Bayesian approach was described by Hauer [21, 22] for studying the index of effectiveness θ of a road measure: the prior probability density function of the parameter θ, reflecting the prior knowledge concerning this parameter, is combined with the likelihood function (probability of the data given the parameter) to obtain the posterior probability density function. The posterior probabilities reflect the revised knowledge about the parameter, given previous knowledge and the data analysed. The method proposed by Hauer, however, is an informative (subjective) Bayes method and presupposes expertise or previously formalised knowledge: the prior probabilities are based on the "elicitation of prevailing opinion about the effectiveness of a treatment" ( [22], p. 289), or possibly on the results of previous studies or meta-analyses. Road safety expertise is limited, however, in the routine situations we consider here, since the study is often carried out by a local road engineer, and not by a road safety specialist. Moreover, the site modification is often singular and not generic (it may combine several treatments, for example: redesigning of islands, resurfacing and marking at a junction site). Therefore, it may be difficult to make use of results from previous meta-analyses. A method coping with this problem was described by Al Masaied et al. [3] : prior probabilities were estimated using a part of the accident data, for both the before and the after periods. In the case of a single site, however, this may lead to very small accident numbers for each data subset. Another way is to use the 'objective' or 'low-informative' Bayesian framework [6, 7, 17, 25, 26] where the prior probabilities are chosen in order to be neutral in some way as regards the possible parameter values, reflecting the lack of previous knowledge. Besides, it can be argued that results based on low-informative approaches are generally easier to communicate to a diverse or uninitiated audience, since, as mentioned by Box and Tiao [13], they represent "what someone who a priori knew very little about an unknown parameter should believe in light of the data" (p. 22). 

In before-after accident studies, it is important to be able to control for regression to the mean bias, which can be done by incorporating some limited information into the prior distribution concerning one component of the vector of parameters (see Section 4). Besides, although such studies are retrospective and not experimental, one should seek to control for the confusing influence of factors other than the road change. To this end, it can be useful to take into consideration a comparison group of similar sites, for example. The method described by Hauer [22] uses a comparison sample, but the calculations are based on approximations which presuppose that the accident counts in the comparison sample are large. The method proposed by Al-Masaied et al. [3] is a simple before-after method without comparison sites. 

In this methodological note, we describe a low-informative Bayesian method adapted to the current practice of beforeafter accident studies concerning a single treated site (or a group of sites considered as a whole). A comparison site (or group of sites) is used in order to control for factors other than the road modification. Practical means of calculation, for a commonly available spreadsheet software package, will also be provided on the author's webpage (http://www. inrets. fr/ur/ ma/Brenac. html). 

2 Data structure and parameters for the before-after study with comparison sites When a comparison site or group of sites is used, the basic data take the form of a 2×2 contingency table (Table 1 ) containing the observed accident counts x i. These counts are considered as observations of independent Poisson variables X i with expected values μ i (unknown). 

Under the assumption of a strong similarity between the treated site and the comparison site, and if the evolution of traffic does not differ between them, the effect of the measure can be represented by the odds ratio

θ expresses the ratio of the 'accidentality' on the treatment site during period II (after modification) to what this 'accidentality' would have been during the same period II, had the site not been modified-here we use the term 'accidentality' in the somewhat unusual sense of the expected value of the accident count. From a practical point of view, an odds ratio of 0. 8, for example, would mean that the effect of the treatment is a 20% reduction in accidentality. The ratio h ¼ m 4 =m 3 reflects the effect of other factors on the evolution from period I to period II, assumed to be common to both the treated and comparison sites (η can be considered as a trend parameter). In other terms, μ 2 and μ 4 can be expressed as follows: μ 2 =μ 1 θ η and μ 4 =μ 3 η. Thus, we are in the presence of a problem with four observations x 1, x 2, x 3, x 4 from four independent random Poisson variables X 1, X 2, X 3 and X 4, and four unknown parameters θ, η, μ 1 and μ 3 with the following relationships:

3 The Bayesian framework

According to the Bayesian approach to statistics, the unknown parameters θ, η, μ 1 and μ 3 are considered as instantiations of variables Θ, Η, M 1 and M 3, treated as random variables, but which in fact reflect our uncertainty about the values of these parameters. Given the observed data x=(x 1, x 2, x 3, x 4 ), given the likelihood function L(x | θ, η, μ 1, μ 3 ) 1 and the joint prior probability density function of the parameters π(θ, η, μ 1, μ 3 ), the application of Bayes' theorem leads to the joint posterior distribution

The joint prior distribution π(θ, η, μ 1, μ 3 ) represents our previous assumptions or knowledge (or lack of knowledge) regarding the parameters (see Section 4). The joint posterior distribution represents our revised knowledge about the parameters, after the observations are taken into account. The likelihood function can be easily derived from the problem formulation given in Section 2:

The parameter of interest is θ. Its posterior probability density function can be obtained by integrating the joint posterior distribution with respect to the three other parameters:

From a practical point of view, however, the most useful result is the posterior cumulative distribution function of Θ,

This cumulative distribution function makes it possible to calculate credible intervals and the probability that the effect studied is lower or higher than a particular value, given the data and prior probabilities.
section_name: Introduction

section_content: In this paper we assume a lack of previous knowledge or sufficient expertise regarding the parameters. Thus, the prior distributions should be low-informative or neutral as regards these parameters. This choice also tends to "let the data speak for themselves", giving a higher importance to the likelihood function in the calculation of posterior probabilities. Two situations should be distinguished, however, according to whether regression to the mean bias is likely or not. Regression to the mean (see, for example, [1] ) occurs when the site was chosen for treatment in consideration of a high accident record. In this case, the count x 1 gives only biased information on the expected value μ 1, and a lowinformative prior distribution for μ 1 would lead to biased results, overestimating the treatment effect. In this situation, other data or information are needed and should be taken into account in the prior distribution of μ 1 (see point 4. 2).
section_name: Low-informative prior distributions

section_content: In many circumstances, regression to the mean bias is unlikely: for example, when the site modification was not decided for safety reasons, but for other purposes (really independent from accident counts). In this case, a lowinformative joint prior distribution can be chosen for the four parameters θ, η, μ 1 and μ 3. The way of selecting lowinformative priors (also called non-informative, objective, default or reference priors) is widely discussed in Bayesian statistics (see the review by Kass and Wassermann [29] ; see also [2, 8, 9, 17, 20, 25, 26, 39] ). We will not enter this debate here since, as mentioned by Ghosh et al. [17], "even though there is no unique objective prior, the posteriors will usually be very similar even with a modest amount of data" (p. 147). In the present paper, for the sake of simplicity, we will only consider the prior obtained by the Jeffreys's general rule 2 [26], which is widely accepted as a "reasonable standard" [29]. For a vector of parameters ξ, the Jeffreys's rule prior is proportional to the square root of the determinant of the Fisher information matrix:

where ∝ denotes proportionality. In this expression, I(ξ) is the Fisher information matrix defined by

where l is the log-likelihood. Applied to our problem, using Eq. 4, this rule leads to the joint prior

Like many non-informative priors, this prior is improper since it does not integrate to a finite value over the parameter space. In Bayesian statistics, however, this is not regarded as a problem, provided that the posterior distribution is proper (i. e., the integral in the denominator of Eq. 3 converges to a finite value).
section_name: Case where regression to the mean bias is unlikely

section_content: In this situation, conventional methods correct for regression to the mean by considering that the site is taken from a population of comparable sites and extracting complementary information from a sample of such sites 3. Each of the accident counts x 1j at these sites, during period I, is considered as an observation from a Poisson variable with mean μ 1j. The μ 1j are assumed to be distributed like a Gamma variable with shape parameter α and scale parameter l (some empirical justifications can be found in the literature [1, 34] ). This Poisson-Gamma structure leads to a negative binomial distribution of the counts x 1j among this sample of sites. Based on the mean m and variance s 2 of this distribution, estimated from the x 1j, it is possible to estimate 4 α and l: α = m 2 /(s 2 -m) and l = m/(s 2 -m). Conventional evaluation methods then replace x 1, the usual estimate of μ 1, by the empirical Bayes estimate μ 1 *=m 2 /s 2 +x 1 (s 2 -m)/s 2 =(α+x 1 )/ (1+l) for the calculation of the odds ratio [16, 23, 36]. This technique has proved to be effective for correcting for regression to the mean bias [38]. 

The equivalent in the 'full' Bayes approach consists in taking the Gamma(α,l) prior distribution for the parameter μ 1 :

In this situation, a joint prior distribution can be obtained by calculating π(θ, η, μ 3 ) with the Jeffreys's rule applied while holding μ 1 fixed (see [29] ), which gives π(θ, η, μ 3 ) ∝ (1/θ) ½ and leads to

(Constants are not taken into account since they would be cancelled anyway as common factors in the denominator and numerator of Eq. 3). This prior is also improper. The estimates of α and l are drawn from accident data at a sample of similar sites (independent from the group of comparison sites), or from an accident model, as described above for conventional methods using empirical Bayes estimates. Although this joint distribution (Eq. 10) uses some prior information concerning μ 1 (through α and l), it remains, however, low-informative in a relative sense, since no prior knowledge is used concerning the parameter of interest θ and the two other parameters η and μ 3. 

5 Posterior probabilities
section_name: Case where regression to the mean bias is likely

section_content: Applying the likelihood function (Eq. 4) and the joint prior distribution given in Eq. 8 to the calculation of the joint posterior distribution (Eq. 3) gives the following expression, after simplification (cancelling of factors present both at the numerator and the denominator):

where

. This latter integral converges to a finite value even when some (or all) of the x i equal zero. Therefore, a proper posterior distribution can always be obtained. The terms in μ 1 and μ 3 are proportional to Gamma density functions, which makes it possible to integrate the expression given in Eq. 11 with respect to μ 1 and μ 3, leading to the joint posterior of θ and η

Þwhere B denotes the Beta function 5. The posterior cumulative distribution function of Θ is then given by

The calculation of this integral is generally not possible by analytical means. We describe in the appendix of this paper a way of calculating it numerically.
section_name: Case where regression to the mean bias is unlikely

section_content: For the prior given in Eq. 10, the same kind of calculations as those described in point 5. 1 lead to the following expression for the posterior cumulative distribution of Θ :

where K

. The calculation leading to this result is not valid if x 1 +α ≤ ½ (which is unlikely: α is a positive parameter and we are in a situation where the treated site was chosen in consideration of a high accident count x 1 ). For the numerical calculation of this integral, see the appendix.
section_name: Case where regression to the mean bias is likely

section_content: From a practical point of view, various useful results can be obtained using the function F Θ (t | x). For example, the lower limit θ LL and upper limit θ UL of a 95% symmetrical credible interval are defined by F Θ (θ LL | x)=0. 025 and F Θ (θ UL | x)=0. 975; the probability, given the data, that θ is contained in this interval is 95%. The median θ med defined by F Θ (θ med | x)=0. 5 gives a point estimate of the odds ratio for which the posterior risks of overestimation and underestimation are equal. The value F Θ (1 | x) represents the posterior probability that θ is lower than 1, i. e. the probability that the treatment is beneficial to safety, given the data and initial assumptions (see Section 2). 

6 Particular cases 6. 1 Group of comparison sites instead of a single comparison site

In this situation, the group of q comparison sites is considered as a whole, with x 3 =Σ x 3k and x 4 =Σ x 4k (where x 3k and x 4k are the accident counts during periods I and II on each comparison site k, with k=1 to q). The aggregated counts x 3 and x 4 are observations from random variables X 3 and X 4 which are Poisson variables (since they are obtained by summing the independent Poisson variables X 3k and X 4k ) with means μ 3 =Σ μ 3k and μ 4 =Σ μ 4k. The calculations described in Sections 3 to 5 are then applied by simply using the aggregated counts x 3 and x 4 and the aggregated means μ 3 and μ 4. The low informative joint prior is given by Eq. 8 or 10. The posterior cumulative distribution function of Θ is then given by Eq. 13 or 14 (with x 3 =Σ x 3k and x 4 =Σ x 4k ).
section_name: Practical uses of the posterior cumulative distribution function of Θ

section_content: The general case of several treated sites, considered independently, with possibly different odds ratios θ i due to heterogeneity in the treatment effect is beyond the purpose of this paper and will be the subject of further publications. Nevertheless, in the simpler situation where a group of treated sites is considered as a whole (with a focus on the overall effect of treatment), the methods described above can be easily adapted. 

Let us consider n treated sites with accident counts x 1i and x 2i (i = 1 to n) during periods I and II, with corresponding means μ 1i and μ 2i, and q comparison sites with accidents counts x 3k and x 4k (k=1 to q) during periods I and II, with corresponding means μ 3k and μ 4k. 

When regression to the mean bias is unlikely, and if we consider the treated sites as a whole (and the comparison sites as a whole), the calculations and results described in Sections 3 to 5 can be applied by simply using the aggregated counts x 1 = Σ x 1i, x 2 = Σ x 2i, x 3 = Σ x 3k, x 4 = Σ x 4k and the corresponding aggregated means, with the prior given in Eq. 8. In this case, the parameter θ represents the overall effect of the programme of treatment. The posterior probabilities are given by Eq. 13. 

When regression to the mean bias appears likely, if the same prior Gamma(α,l) distribution can be assumed for the mean μ 1i of each treated site i, the prior distribution of the overall mean μ 1 = Σμ 1i is a Gamma(nα, l) distribution (using the classical property of the sum of independent Gamma variables with same scale parameter l). 

Considering the treated sites as a whole (and the comparison sites as a whole), and considering θ as the overall effect, the joint prior distribution becomes

where μ 3 = Σμ 3k. The posterior cumulative distribution function of Θ is then

where K 00 ¼

and where n is the number of treated sites, x 1 =Σ x 1i, x 2 =Σ x 2i, x 3 =Σ x 3k and x 4 =Σ x 4k. 

7 Examples of application
section_name: Multiple treated sites

section_content: We describe here the case of an urban section of road where the infrastructure was largely modified in order to enhance the quality of local urban life. Raised median islands, small roundabouts, speed humps and raised tables were implemented in 2000 on this section of a main urban road in a town of 40,000 inhabitants (length of the treated section: 700 m). All the unmodified sections of the main roads in this town were taken as a comparison group of sites. The comparability between the treated site and the comparison group of sites was verified by comparing the yearly injury accident counts for the 1989-1999 period. The 'before' period is the five-year period from 1995 to 1999. The 'after' period is the five-year period from 2001 to 2005. The presence of regression to the mean bias was considered to be unlikely for the following reasons: this project was not decided for safety reasons, and the proportion of accidents during the 1995-1999 period relative to 1989-1999 was not higher on the treated site as compared to all the unmodified sections of main roads in this town. For the 'before' period, 16 injury accidents occurred on the treated site and 61 injury accidents occurred on the comparison group of sites. For the 'after' period, 3 injury accidents occurred on the treated site, and 46 injury accidents occurred on the comparison group of sites. 

The calculations applied to these data (x 1 =16, x 2 =3, x 3 = 61, x 4 =46) with the low-informative prior based on the Jeffreys's rule (Eq. 8) give the following results based on the posterior cumulative distribution function of Θ : 95% symmetrical credible interval: 0. 062 to 0. 815 Median: 0. 259 Posterior probability that θ <1: 0. 990 Figure 1 shows an example of spreadsheet screen for the posterior probability calculation (see the appendix), applied to these data. 

These results suggest a beneficial effect on safety. They can be compared to the results that would be obtained by conventional statistical methods. Nevertheless, as mentioned in the introduction, Bayesian and non-Bayesian concepts (like credible interval and confidence interval) can not be interpreted in the same way 6. In this example, the usual unconditional maximum likelihood estimator of the odds ratio, with the related approximate 95% confidence interval (Woolf interval), would lead to the following results:

Woolf 95% confidence interval: 0. 068 to 0. 904

In this example, a practitioner would probably conclude in favour of a positive effect on safety, from both these Bayesian and non-Bayesian results.
section_name: Example 1: Safety effect of redesigning an urban road section

section_content: This example deals with a priority intersection on a main rural two-lane road. This crossroads was modified in 1986 (installation of median raised islands, marking) for safety reasons. Therefore, regression to the mean is likely to occur. At this junction, 14 injury accidents occurred during the three-year period before the treatment. During the threeyear period following the treatment, 4 injury accidents occurred. 

This evolution was compared to the evolution observed at a set of 11 similar intersections on main rural two-lane roads in the same region, used as a comparison group of sites. At these sites, considered as a whole, 33 injury accidents occurred during the before period and 22 injury accidents occurred during the after period. 

The calculations applied to these data (x 1 =14, x 2 =4, x 3 = 33, x 4 =22), using the low-informative prior distribution given by Eq. 8 (Jeffreys's rule prior), would lead to the following results based on the posterior cumulative distribution function of Θ : 95% symmetrical credible interval: 0. 117 to 1. 389 Median: 0. 439 Posterior probability that θ <1: 0. 917

Due to the high risk of regression to the mean in this case, however, these results are probably biased. In order to correct for this regression to the mean bias, it is necessary to use a more 'informed' prior, concerning the parameter μ 1 (see Section 4). To this end, the parameters α and l of a prior Gamma distribution for μ 1 have to be estimated. By applying an accident model (which was established at a national level [14] ) to the characteristics of this junction (traffic volumes, number of arms, number of traffic lanes), as mentioned in Section 4. 2., it is possible to obtain the mean m=3. 55 and variance s 2 =15. 90 of the accident counts for a virtual population of similar sites during the same period I. On this basis, we can calculate the estimates α=1. 02 and l=0. 29. The joint prior given by Eq. 10 is then precisely defined and leads to the following results, in terms of posterior probabilities: 95% symmetrical credible interval: 0. 151 to 1. 789 Median: 0. 566

Posterior probability that θ<1: 0. 828

These results show that, in this case, the safety effect is in reality smaller than indicated by the biased results obtained with the low-informative prior given by Eq. 8. The median of the posterior distribution (0. 566) can be used as a point estimate of the odds ratio (where the posterior probabilities of overestimation and underestimation are equal). This value corresponds to an accident reduction of approximately 43%. The 95% credible interval, however, is large and the beneficial effect of the treatment remains uncertain. 

Using the same data, a more conventional approach would lead, for example, to the maximum likelihood estimate θ ML *=0. 429 (without controlling for regression to the mean), or to a corrected estimate of 0. 515 based on the empirical Bayes estimate of μ 1 [16, 36].
section_name: Example 2: Safety effect of a rural crossroads modification

section_content: This example is based on some of the data published in an article by Leden et al. [33], dealing with the effect of resurfacing on friction, speeds and safety on main roads in Finland. The treated sites are all sections on main roads (in the south of Finland) which were resurfaced in 1991. The comparison sites are all the untreated main roads in the same region. Due to the particular road conditions in winter in Finland, only the effects on the non-winter period (from April 1 to September 30) are studied. Regression to the mean bias is considered to be unlikely, since "sections were selected for treatment on a routine maintenance base" [33] (p. 82) and not for safety reasons. We consider the treated sites as a whole, and the comparison sites as a whole. The parameter θ thus represents the overall effect of the Fig. 1 Example of spreadsheet screen for the posterior probability calculation treatment on the set of sites. The following data concern the 'before' period from April to September 1990 and the 'after' period from April to September 1992. Before and after injury accident counts are x 1 =80 and x 2 =74 on the treated sites, and x 3 =931 and x 4 =779 on the comparison sites. Based on the Jeffreys's rule prior, the results are as follows, in terms of posterior probabilities: 95% symmetrical credible interval: 0. 794 to 1. 537 Median:

1. 106

Posterior probability that θ<1: 0. 275

One can note the proximity of these results with the following results which would be obtained with a conventional frequentist approach:

Woolf 95% confidence interval: 0. 794 to 1. 537 This proximity is not surprising: posterior credible intervals based on the Jeffreys's rule prior are frequently close to frequentist confidence intervals in large sample conditions [17, 40] although they do not have the same meaning. 

Based on these results, the posterior median estimate of θ would suggest a slight detrimental effect on safety (increase of accidentality of approximately 11%), but no certain conclusion can be drawn since the 95% credible interval is large. Based on the posterior probability that θ<1 (approximately 28%), however, one could say that the probability that the treatment increases the accidentality, given the data and assumptions, is 72%. No equivalent result from a conventional statistical analysis could lead to this kind of interpretation, except if one wrongly interprets a p-value as a posterior probability. A possible increase of accidentality could be explained by the fact that resurfacing tends to increase the average speeds, at least when the road is dry, as shown by Leden et al. [33].
section_name: Example 3: Safety effect of resurfacing on main roads

section_content: In this note, we described a low-informative Bayesian method for before-after accident studies, using a comparison site or group of sites, and giving the possibility of correcting for regression to the mean bias. The aim was to provide a methodological basis for routine evaluation studies, often applied to a single treated site, and in conditions of limited resources in terms of time and expertise. As compared to conventional statistics, the Bayesian approach is less subject to misuse and misinterpretation by practitioners with limited statistical experience. The low-informative or objective Bayesian methods seem appropriate in routine evaluation studies, where expertise or previous knowledge are often limited or hard to formalise. As shown in Sections 2 to 6, a relatively simple method, based on the Jeffreys's rule prior considered as a "reasonable standard", can be implemented without major difficulties. Posterior distributions are proper. The numerical calculation of posterior probabilities can be done without using Monte-Carlo methods nor specialised software tools. The examples given in Section 7 show that the results can be analysed in a direct way, without the high risk of misinterpretation involved in the analysis of frequentist results. 

Further developments, however, are still needed. Although this method seems to be transferable to engineers for common practice, further work is necessary in order to provide a simple, didactic description of the Bayesian line of reasoning, with minimal use of mathematical formalisms, appropriate for communicating this approach to practitioners. Concerning the practical means for calculating the posterior probabilities, the spreadsheet mentioned in the appendix (for a common spreadsheet software package) will be made available on our website. 

The proposed method has limitations, of course. Retrospective before-after studies are not randomised experiments and the validity of their results is based on the assumption that the treated and comparison sites are similar. Before-after studies based on multivariate generalised linear models make it possible to better control for the influence of differences between treated and comparison sites. But such methodologies would generally involve a thorough data collection and analysis on a large sample of sites, which seems hard to implement by practitioners in the routine situations we considered in this paper. The comparability of treated and comparison sites, however, can be checked by examining their accident history, when accident data are available for a long period before the treatment (see [23] ). A Bayesian approach to this subject could be studied. Besides, other developments could contribute to extending the field of application of the proposed method: in this paper, we only dealt with the case of a single treated site (or a group of sites treated as a whole, with a focus on the overall effect of the programme of treatment), with a comparison site or group of sites. The case of multiple treated sites considered independently and with possibly different odds ratios remains to be dealt with. However, this would involve an increased complexity and more difficulties for practitioners. 

We hope this methodological note will contribute to an increased use of the Bayesian approach, which is more in accordance with the expectations and intuitions of nonstatisticians, in the current practice of before-after accident studies.
section_name: Discussion and conclusion

section_content: In the case where regression to the mean bias is unlikely, this function is given by Eq. 13. It can be written in a more generic form:

where the parameters a, b, c and d take the values a=x Beta cumulative distribution function and its inverse are commonly available in spreadsheet software, and this integral can be calculated without using specialised tools (see below). In the case where regression to the mean bias is suspected, the posterior cumulative distribution function of Θ is preferably calculated from Eq. 14. This equation can be written in the following form:

where the parameters a, b, c and d take the values a=x 2 +½, b=x 1 +α-½, c=x 4 +½, d=x 3 +½ for the prior given by Eq. 10, and the values a = Σx 2i +½, b = Σx 1i +nα-½, c = Σx 4k +½, d = Σx 3k +½ in case of multiple sites considered as a whole, with the prior given in Eq. 15. After the three successive changes of variable z(η)= η /(1+η), ω(θ) = θ z / ((θ-1-l) z+1 + l) and lastly u(z)=Betacdf c,d (z), this integral becomes

We will propose (on our website) a spreadsheet calculating F Θ (t | x) for any given value of t, given the accident counts and the prior choice (among the priors mentioned in this paper), for a common spreadsheet software package. The calculation of the integrals of Eqs. 18 and 20 is based on a classical trapezoidal quadrature method, with an increase of the number of calculation points in the vicinity of 0 and 1. The reliability of the results was carefully checked by comparing them to the results obtained with a more powerful software tool using other quadrature methods (adaptive Simpson quadrature and Lobatto quadrature), for a large range of possible values for accident counts x 1, x 2, x 3, and x 4. NB: Simplified calculations can be equivalently used in the special case where the counts x 3 and x 4 are very large, since the trend parameter η can then be considered as non-random and approximately equal to x 4 / x 3. The problem then reduces to a two-parameter problem (θ and μ 1 ), with two random observations x 1 and x 2. In this situation, the Jeffreys's rule prior is again π(θ, μ 1 ) ∝ (1/θ) ½. In the case where regression to the mean should be taken into account, the method described in Section 4 leads to the following prior: π(θ, μ 1 ) ∝ (1/θ) ½ μ 1 α-1 exp(-lμ 1 ). With the Jeffreys's rule prior, the same kind of calculations as those described in Section 3 to 6 lead to the posterior cumulative distribution function

This is a Beta-prime distribution function with three parameters x 2 +½, x 1 +½ and 1/η. After a change of variable z(θ) = θ /(θ +1/η), this integral can be written in the form of a Beta cumulative distribution function:

where a=x 2 +½ and b=x 1 +½. In the case where the prior is π(θ, μ 1 ) ∝ (1/θ) ½ μ 1 α-1 exp(-lμ 1 ), the calculation leads to the following result:

where a=x 2 +½ and b=x 1 +α-½.
section_name: Appendix: Numerical calculation of the posterior cumulative distribution function of Θ

section_content: Acknowledgements The author would like to thank Sylvie Després ( Université Paris-Nord ) and three anonymous referees for their helpful comments.