hiseliusestimatingrelationshipbetweenaccident2004.md 44 KB


category: literaturenote citekey: hiseliusestimatingrelationshipbetweenaccident2004 title: Estimating the relationship between accident frequency and homogeneous and inhomogeneous traffic flows authors: "Hiselius, Lena Winslott" year: 2004 date: 2004-00-00 2004 doi: 10.1016/j.aap.2003.11.002 publication: "Accident Analysis & Prevention" url: "https://doi.org/10.1016/j.aap.2003.11.002" zotero_key: WHT97JG4 zotero_storage: VA46EBG7 collections: transpordi_planeerimine folder: 001_artiklid/Road Safety Indicators firstAuthor: "Hiselius, Lena Winslott"

status: converted

Accident Analysis and Prevention 36 (2004) 985–992

Estimating the relationship between accident frequency and homogeneous and inhomogeneous traffic flows

Lena Winslott Hiselius∗

Department of Economics or Department of Technology and Society, Lund University, P.O. Box 118, Lund 22100, Sweden Received 5 January 2002; received in revised form 19 November 2003; accepted 24 November 2003

Abstract

This paper estimates the relationship between accident frequency and the traffic flow empirically treating the hourly traffic flow in two different ways, as consisting of homogenous vehicles and as consisting of cars and lorries. Rural roads in Sweden are studied using Poisson and Negative Binominal regression models. It is found that important information is lost if no consideration is taken to differences between vehicle types when estimating the marginal effect of the traffic flow. The accident rate decreases when the traffic flow is treated as if homogeneous. However, when cars are studied separately the result suggests that the accident rate is constant or increases. The result with respect to lorries is reversed, indicating a decreasing number of accidents as the number of lorries increases. © 2004 Elsevier Ltd. All rights reserved.

Keywords: Accident prediction models; Road accidents; Hourly traffic flow; Poisson; Negative Binomial

1. Introduction

Studies that analyse the relationship between accident frequency and traffic flows show a great variation in their results. An early report from Vickery (1968) suggests that the marginal accident rate is 1.5 times the average accident rate. On the other hand, Vitaliano and Held (1991) cannot detect any significant increase in the accident rate as the traffic flow increases and, according to Hauer and Bamfo (1997), and a majority of the results reviewed in Ardekani et al. (1997), the accident rate even decreases with an increasing number of vehicles. Furthermore, only a few empirical studies analyse the effect of the traffic flow separated for different road user groups. One example is the study by Jovanis and Chang (1986), where they found that the accident frequency increased with the number of car and lorry miles travelled. However, it is difficult to distinguish the marginal effect of another car or lorry as the number of miles travelled varies in their study.

There is, thus, an interest in studying the effect of the traffic flow more thoroughly, and the aim of this article is to estimate the relationship between the accident frequency and the traffic flow treated as both homogeneous, i.e. consisting of homogeneous vehicles, and inhomogeneous, i.e.

E-mail address: lena.hiselius@tft.lth.se (L.W. Hiselius).

consisting of cars and lorries. The data is restricted to accidents on sections in rural areas of Sweden.

2. Model

The model seeks to describe the relationship between the number of vehicles per hour and the accident frequency. It is developed under the assumption that a limited and homogeneous road system is studied and only the traffic flow affects the number of accidents. In the first case the flow of vehicles per hour, q, is assumed to be homogeneous, and with all vehicles influencing the occurrence of an accident the same way. In the second case the traffic flow is defined as the number of cars and lorries per hour, qc and *q*l. The accident frequency, Z, is calculated as the expected number of accidents per unit of time and kilometer in order to apply the model to time periods and road systems of different length.

3. Data

Data is collected from 83 road sections in rural areas of Sweden, where the number of passing vehicles is continuously counted by the Swedish National Road Administration. The assumption is made that the counted traffic flow at a stationary place is valid along the section. Information is collected on police reported accidents with personal

∗ Fax: +46-46-12-32-72.

Table 1 Number of accidents in the data set

| Road
type | | Number of accidents | Total number of | | | |--------------|-------|---------------------|-----------------|-------------------|--| | | Total | Single-vehicle | Multi-vehicle | daytime accidents | | | I | 59 | 29 | 30 | 49 | | | II | 83 | 52 | 31 | 51 | | | III | 179 | 107 | 72 | 111 | | | IV | 186 | 121 | 65 | 129 | |

casualties that occurred on the studied road sections from 1989 to the middle of 1995. Given the time and date of the accidents, the hourly traffic flow that prevailed at the time of each accident is obtained. Furthermore, in order to calculate traffic flow frequencies, information on the number of hours that each traffic flow was observed during 1990 is also collected on the assumption that this is a representative year.

Only accidents occurring on sections without intersections are included. Accidents involving animals are excluded together with accidents that may be considered specific for each road section. Since it is not possible to obtain information on driving speed, the analysis will be made bearing in mind that the estimated effect of the traffic flow may also be an effect of speed adjustment. Moreover, there are several other factors likely to influence the occurrence of accidents, e.g. weather, road conditions, type of vehicle and drivers' characteristics. In order to take some of these factors into account, daylight accidents are studied separately.

Four road types are distinguished. Road type I with speed limit 70 km/h and road width 6–9.7 m, road type II with speed limit 90 km/h and road width 6–7.9 m, road type III with speed limit 90 or 110 km/h and road width 8–13 m, without separated road lanes, and road type IV, motorways, with speed limit 90 or 110 km/h. The numbers of accidents that occurred on the studied sections are presented in Table 1.

For the road types I, II and III the traffic flow is counted as the number of vehicles per hour driving in both directions. Thus, there is no information on the traffic flow per road lane. This is a problem since depending on the distribution of the traffic flow, different types of accidents may occur. The same situation occurs when studying road type IV. Even though the traffic flow is counted for each direction for this road type, there is no information on the traffic in each lane.

4. Accident frequency

In order to analyse the accident frequency, Z, the distribution model of the number of accidents, Y, is first discussed. Several authors use a Poisson distribution to approximate the number of accidents occurring during a certain time interval, e.g. Hall (1986) and Nicholson and Wong (1993). The assumption is based on the fact that although the total number of accidents is reasonably stable over a period, each accident is unpredictable. A reasonable assumption is furthermore that accidents are independent events. However, it may be argued that individual collisions in a multi-collision accident are not mutually independent events. This problem is avoided here, since a traffic accident in this study is defined without consideration being given to the number of vehicles involved. Let Yi denote the number of accidents occurring on a specific road site during a given time period. If we assume that the number of accidents follows a Poisson distribution with expected value, and thus variance, equal to λi, the probability that the number of accidents during this period will be equal to *y*i, may be written as

$$P(Y_i = y_i) = \frac{\lambda_i^{y_i} e^{-\lambda_i}}{y_i!}$$ $y_i = 0, 1, 2, 3, ...$

There is, however, no reason to believe that λi is same at all sites, and there will usually be a variation of the λi around a mean, λ. The variation of the Yi (the number of accidents occurring at different sites) will thus comprise both the variation of yi around λi and the variation of λi around λ. The variance will therefore tend to be larger than the mean, i.e. we can observe so called over dispersion with respect to the Poisson distribution. The variation in λi is often assumed to be Gamma distributed. Assuming λi to have a Gamma distributed parameter κ, it can be shown that the number of accidents, Y, follows a Negative Binomial distribution, e.g. in Kulmala (1995). The probability density function for Y can now be written as

$$P(Y = y) = \frac{\Gamma(\kappa + y)}{\Gamma(\kappa)y!} \left(\frac{\kappa}{\kappa + \lambda}\right)^{\kappa} \left(\frac{\lambda}{\kappa + \lambda}\right)^{y}$$ $$y = 0, 1, 2, 3, \dots$$

and variance

$$\sigma^2 = \lambda + \frac{\lambda^2}{\kappa}$$

The Negative Binomial distribution is thus an extension of the Poisson distribution, where the gamma parameter, κ, determines the shape of the distribution. As κ → ∞, the rate of over dispersion declines and the distribution approaches the Poisson.

So far the distribution model of the number of accidents has been discussed. However, when analysing accident frequency, there is also need for a regression model that can describe this phenomenon. An exponential function is a common formulation, since this function ensures that the expected number of accidents is a positive number. Furthermore, measuring the explanatory variables on the logarithmic scale, Gj = ln(Xj), the accident frequency may be written as a multiplicative function for which the value of the exponents can be estimated directly

$$E(Y) = \exp\left[\sum \beta_j \ln(X_j)\right] = \Pi X_j^{\beta_j}$$

The result of a regression model is, however, often dependent on the choice of model function, Hauer and Bamfo (1997). Using an exponential regression model, as described

Table 2 Statistics on traffic flow data

| Road type | Number of traffic flow intervals | | Average traffic flow | | | Maximum traffic flow | | | | |-----------|----------------------------------|-------------------------------|-----------------------------|-------------------------------|---------|-----------------------------|-------------------------------|---------|--| | | Homogeneous
traffic flow | Inhomogeneous
traffic flow | Homogeneous
traffic flow | Inhomogeneous
traffic flow | | Homogeneous
traffic flow | Inhomogeneous
traffic flow | | | | | | | | Cars | Lorries | | Cars | Lorries | | | I | 98 | 408 | 153 | 139 | 12 | 2712 | 2328 | 312 | | | II | 38 | 159 | 94 | 81 | 8 | 987 | 938 | 132 | | | III | 131 | 595 | 250 | 232 | 24 | 1687 | 1645 | 267 | | | IV | 121 | 458 | 550 | 498 | 55 | 3088 | 2937 | 326 | |

Table 3 Statistics on total number of accidents (Tot), single-vehicle (Sing) and multi-vehicle accidents (Mult)

Road type Number of accidents per traffic flow interval, average Number of accidents per traffic flow interval, variance
Homogeneous traffic
flow
Inhomogeneous
traffic flow
Homogeneous traffic
flow
Inhomogeneous traffic
flow
Tot Sing Mult Tot Sing Mult Tot Sing Mult Tot Sing Mult
I 0.6 0.3 0.3 0.1 0.1 0.1 1.4 0.7 0.4 0.4 0.2 0.1
II 2.2 1.4 0.8 0.5 0.3 0.2 4.4 3.8 0.7 1.7 0.9 0.4
III 1.4 0.8 0.5 0.3 0.2 0.1 4.5 2.5 0.9 0.6 0.3 0.2
IV 1.5 1.0 0.5 0.4 0.3 0.1 3.8 2.2 0.7 1.0 0.5 0.2

previously, restricts the model to the description of monotonic functions, i.e. the functions may be either increasing or decreasing, not both. The use of this model function is supported by studying the data visually and by analysing the residuals.

4.1. Variables and regression models

In this study, the relationship between the accident frequency and the traffic flow is estimated empirically. In order to apply the result to time periods and road systems of different length, the accident frequency, Z, is defined as the expected number of accidents, E(Y), per hour and kilometre. A parameter, H, is also defined, see Appendix A. This parameter describes the exposure, i.e. the number of hours and kilometres that each road type has been studied.

The accident frequency cannot be estimated as a function of the traffic flow directly, since the variable describing the accident frequency is not Poisson or Negative Binomial distributed. The regression model is therefore rewritten as E(Yt) = Zt ×Ht so that the accident frequency may be estimated separately as a function of the traffic flow. The complete regression model for the case of homogeneous traffic flow is estimated as

$$E(Y_1) = \exp[\beta_1 + \beta_2 \ln(q_1) + \beta_3 \ln(H_1)] = [e^{\beta_1} \times (q_1)^{\beta_2}]H_1$$ with the restriction, $\beta_3 = 1$ . When treating the traffic flow as consisting of cars and lorries, the following regression model is estimated

$$\begin{split} E(Y_t) &= \exp[\beta_1 + \beta_2 \ln(q_c)_t + \beta_3 \ln(q_l)_t + \beta_4 \ln(H)_t] \ &= \left| e^{\beta_1} \times (q_c)_t^{\beta_2} \times (q_l)_t^{\beta_3} \right| H_t \end{split}$$

with the restriction β4 = 1.

The regression models are estimated applying both the Poisson and the Negative Binomial distribution. The statistical program LIMDEP version 7 is used. In order to evaluate the goodness-of-fit, four indicators are studied; the scaled deviance statistic, the over dispersion parameter, κ, estimated in the Negative Binomial regression model and two test functions for over dispersion.

5. Results

The regression analysis of this study is based on aggregates of hours with similar traffic flows, see Appendix A. Traffic flow statistics per road type are presented in Table 2 together with the number of traffic flow intervals used. In Table 3, statistics on the aggregated number of accidents per traffic flow interval is presented for various accident types1. The average number of accidents per traffic flow interval is generally low when inhomogeneous traffic flows or different accidents types are studied.

The results from the regression analysis generally suggest a good fit for both the Poisson and the Negative Binominal models. Distributional assumptions do not seem to affect the results since there are only small differences. The test statistics that are used for over dispersion occasionally contradict the size of the dispersion parameter. Contradictions tend to arise when the scaled deviance is almost the same for the Poisson and the Negative Binominal regression models, suggesting a good fit for both. Since the estimated coefficients are practically the same, only the results of the

1 Statistics for daytime traffic flows and accidents are left out since these statistics corresponds largely to the presented statistics.

Table 6Expected number of accidents occurring all day and in daylight, inhomogeneous traffic flowa

| | Ex
ted
mb
of
pec
nu
er | ide
nts
acc | | | Ex
ted
mb
of
ide
in
day
lig
ht
nts
pec
nu
er
acc | | | | | |----------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------------------|--| | | Ro
ad
e I
typ | Ro
ad
e I
I
typ | Ro
ad
e I
II
typ | Ro
ad
e I
V
typ | Ro
ad
e I
typ | Ro
ad
e I
I
typ | Ro
ad
e I
II
typ | Ro
ad
e I
V
typ | | | Co
nst
ant
ln(
)
qc
ln(
)
ql
ln(
H) | b
−3
1.5
4 (
−1
6.4
2)
b
3.6
2 (
7.1
0)
b
−2
.66
(−
5.4
0)
1.0
0 | b
−2
3.4
8 (
−3
3.0
3)
b
1.4
3 (
7.2
8)
b
−0
.94
(−
5.6
9)
1.0
0 | b
−2
1.1
6 (
−5
0.5
6)
b
1.0
4 (
9.9
6)
b
−0
.77
(−
9.0
7)
1.0
0 | b
−2
2.8
4 (
−3
5.3
8)
b
1.9
0 (
13.
04
)
b
−1
.84
(−
14.
44
)
1.0
0 | b
−3
1.1
4 (
−1
5.6
5)
b
3.6
0 (
6.9
5)
b
−2
.66
(−
5.4
1)
1.0
0 | b
−2
4.1
8 (
−2
3.8
4)
b
1.6
8 (
6.1
8)
b
−1
.26
(−
5.4
0)
1.0
0 | b
−2
1.4
2 (
−2
9.9
7)
b
1.1
6 (
7.3
6)
b
−0
.92
(−
7.8
5)
1.0
0 | b
−1
9.5
2 (
−2
1.6
0)
b
1.3
0 (
7.0
2)
b
−1
.61
(−
11.
39)
1.0
0 | | | d.f
Sc
ale
d d
evi
atio
n
LM
-st
at
t-S
isti
1
tat
cs
t-S
isti
2
tat
cs | 40
4
56
.31
61
c
.97
0.6
0
3.6
9c | 15
5
35
.44
17.
52c
−1
.21
0.0
3 | 59
1
16
0.2
3
12.
55c
1.4
6
0.9
4 | 45
4
15
3.9
7
14.
40
c
−0
.39
0.6
3 | 39
6
55
.49
54
c
.67
0.7
4
4.2
8c | 15
4
26
.54
7.4
6
−2
.11
c
−3
.26
c | 57
1
11
3.8
8
5.5
4
−2
.11
c
−2
.43
c | 45
4
17
8.1
9
19
.70
c
1.2
5
1.7
9 | |

Values in parenthesis are t-statistics.

Table 7Expected number of single and multi-vehicle accidents, inhomogeneous traffic flowa

| | Ex
ted
mb
of
pec
nu
er | sin
le-
veh
icl
cci
den
ts
g
e a | | | Ex
ted
mb
of
lti-
veh
icl
cci
den
ts
pec
nu
er
mu
e a | | | | | |--------------------------------------|------------------------------------------|----------------------------------------------------------|------------------------------------------|------------------------------------------|-------------------------------------------------------------------------------------------------|------------------------------------------|------------------------------------------|------------------------------------------|--| | | Ro
ad
e I
typ | Ro
ad
e I
I
typ | Ro
ad
e I
II
typ | Ro
ad
e I
V
typ | Ro
ad
e I
typ | Ro
ad
e I
I
typ | Ro
ad
e I
II
typ | Ro
ad
e I
V
typ | | | Co
nst
ant | b
−2
9.0
0 (
−1
2.5
7) | b
−2
2.8
5 (
−2
7.2
3) | b
−2
0.3
6 (
−4
2.9
9) | b
−2
2.0
8 (
−3
0.0
4) | b
−3
7.4
1 (
−1
0.8
0) | b
−2
6.5
8 (
−2
0.1
8) | b
−2
4.9
9 (
−2
9.0
2) | b
−2
6.8
9 (
−2
0.4
7) | | | ln(
)
qc | b
2.9
0 (
4.6
5) | b
1.1
4 (
4.8
2) | b
0.7
6 (
6.1
0) | b
1.7
0 (
9.9
2) | b
4.6
9 (
5.2
3) | b
1.9
6 (
5.5
7) | b
1.6
1 (
8.3
8) | b
2.4
0 (
8.6
4) | | | ln(
)
ql | b
−2
.29
(−
3.8
3) | b
−0
.85
(−
4.3
4) | b
−0
.73
(−
7.0
9) | b
−1
.82
(−
11.
91
) | b
−3
.15
(−
3.6
7) | b
−1
.11
(−
3.6
8) | b
−0
.83
(−
5.6
9) | b
−1
.90
(−
8.2
3) | | | ln(
H) | 1.0
0 | 1.0
0 | 1.0
0 | 1.0
0 | 1.0
0 | 1.0
0 | 1.0
0 | 1.0
0 | | | d.f | 40
4 | 15
5 | 59
1 | 45
4 | 40
4 | 15
5 | 59
1 | 45
4 | | | Sc
ale
d d
evi
atio
n | 59
.34 | 47
.30 | 22
8.5
5 | 17
4.2
7 | 78
.34 | 49
.39 | 19
2.3
0 | 19
0.3
1 | | | LM
-st
at | c
13
.54 | 95c
13. | 8.9
3 | c
11
.15 | 4c
9.9 | 7.1
9 | c
13
.67 | c
56
.60 | | | t-S
isti
1
tat
cs | 0.0
2 | −0
.03 | 0.0
3 | −0
.04 | −0
.01 | −0
.05 | 0.0
4 | 0.0
2 | | | t-S
isti
2
tat
cs | 0c
0.3 | 0.0
2 | 0c
0.2 | −0
.02 | 0.1
2 | c
−0
.22 | −0
.09 | 0.0
8 | |

Values in parenthesis are t-statistics.

Significantly different from zero at the 5% level.

The null-hypothesis of over dispersion is significantly rejected at the 5% level.

Significantly different from zero at the 5% level.

The null-hypothesis of over dispersion is significantly rejected at the 5% level.

Fig. 1–4. Number of accidents per hour and kilometre (×1E7) at averge number of cars per hour +/−20%. Remark: The number of dots does not correspond to the number of accidents since more than one accident have occurred at some traffic flows. The total number of accidents in the figure of road type I–IV is 8, 18, 15 and 37.

Fig. 5–8. Number of accidents per hour and kilometre (×1E7) at averge number of lorries per hour +/−20%. Remark: The total number of accidents in the figure of road type I–IV is 16, 10, 21 and 14.

Poisson regression model will be presented. A complete report of the results is to be found in Winslott (1998).

5.1. Homogeneous traffic flow

The results of the regression analysis are presented in Table 4 for accidents occurring all day and in daylight only. All estimated parameter values are significantly different from zero. Furthermore, for all road types, except for road type I, the estimated exponent is significantly different from 1. Thus, we may reject the hypothesis that the expected accident frequency increases in proportion with the traffic flow for these road types, i.e. that the accident rate is constant. Instead, an additional vehicle lowers the accident rate and increases the traffic safety.

For the road types I, II and III the estimated parameters are lower than when studying accidents that have occurred throughout the day. The difference is, however, not significant. The result for road type IV suggests that the model of the expected number of accidents occurring in daylight is not a correct choice. The estimated value of the exponent should be interpreted with caution.2

The results for single and multi-vehicle accidents are presented in Table 5. The results differ in some respects from the case when all accident types are studied together. The estimated values of the exponent of the traffic flow are significantly larger for multi-vehicle accidents than single-vehicle accidents. For multi-vehicle accidents, we cannot reject the hypothesis that the expected number of multi-vehicle accidents increases in proportion with the traffic flow. However,

2 The chosen regression model may only describe a monotonic function, and when the graph of the number of daytime accidents versus the number of vehicles per hour is studied for road type IV, there is no sign of a monotonic relationship, either increasing or decreasing.

for single-vehicle accidents, the estimated exponents of the traffic flow are significantly less than one for all road types.

5.2. Inhomogeneous traffic flow

Table 6 shows the results for all accidents and accidents occurring in daylight only when separating the traffic flow for cars and lorries. The estimated coefficients are significantly different from zero for all road types. For road type I and IV the exponents are significantly different from 1, indicating that the expected number of accidents increases more than proportionally with the number of cars per hour, i.e. that the accident rate increases. However, for the road types II and III the null-hypothesis cannot be rejected. The estimated exponent for the flow of lorries is generally negative and different from zero. This suggests that an increasing number of lorries per hour will lower the expected number of accidents independently of the flow of cars. When studying accidents in daylight, the results are not significantly different.

The estimates when studying single and multi-vehicle accidents are presented in Table 7. The value of the exponent for the number of cars per hour is generally different from 1 when both single and multi-vehicle accidents are studied. Accordingly, the expected number of accidents increases more than proportionally for both single and multi-vehicle accidents. However, the exponent for cars is significantly lower for single-vehicle accidents than for multi-vehicle accidents. The exponent for the number of lorries per hour is also negative here. The value of the exponent is generally, however not significantly, lower when multi-vehicle accidents are studied compared to when all accidents are studied together.

The difference between the studied road types is also tested. The results for the road types I and IV are significantly higher and lower than the results for the road type II and III, respectively. This corresponds to the results when all accidents are studied together.

The result using an inhomogeneous traffic flow differs greatly to that of homogeneous traffic flow. Since the estimates raise certain questions the data is also studied visually. Figs. 1–4 present the number of accidents per hour and kilometre in relation to the flow of lorries per hour. Figs. 5–8 present the number of accidents per hour and kilometre in relation to the flow of cars per hour. The figures are drawn for the average number of cars (Figs. 1–4) and lorries (Figs. 5–8) per hour for each road type±20%.3 Generally, there are few accidents in the data set and when restricting the analysis to certain traffic flows, the number of accidents is reduced even further. In spite of a limited number of observations, the plots can be argued to indicate the type of relation between the accident frequency and the flow of cars and lorries.

Data presented in Figs. 1–4 indicates that the number of accidents is constant or decreasing with an increasing number of lorries on the road, whereas in Figs. 5–8 the number of accidents is increasing with an increasing number of cars. Consequently, the visual analysis is not in conflict with the estimates received in the regression analysis.

6. Discussion

The estimated relationship between the expected number of accidents per hour and kilometre, and the traffic flow differs considerably depending on whether different types of traffic modes are considered or not, i.e. if a homogeneous or a inhomogeneous traffic flow analysis is carried out. Differences appear both in relation to accidents occurring throughout the day, accidents occurring in daylight, singleand multi-vehicle accidents. In the homogeneous traffic flow analysis regarding accidents that have occurred throughout the day, we reject the null-hypothesis that the expected number of accidents increases proportionally with the traffic flow for the road types II, III and IV. The value of the exponent is less than 1. This result is in line with a number of studies reviewed in Satterthwaite (1981), Ardekani et al. (1997) and Hauer and Bamfo (1997).

As cars constitute the main part of the traffic flow, one may expect the outcome, with respect to the number of cars per hour, to be similar to that of the homogeneous traffic flow analysis. However, here the suggestion is that the expected number of accidents increases in proportion or more, to the number of cars per hour. Furthermore, studying the effect of an increasing number of lorries one might expect that the presence of more lorries on the road would increase the number of accidents since the incidence of possibly dangerous overtaking manoeuvres increases. The result of this study indicates, however, that at a given number of cars per hour, the expected number of accidents will decrease with increasing number of lorries per hour. This safety effect of lorries may be regarded as counter intuitive at first glance. It is, however, possible that the number of accidents decreases as an effect of a speed reduction as the speed limit for lorries is lower than for cars. An increasing number of lorries per hour will accordingly slow down the average speed. The result with respect to the flow of lorries may also be a result of people's unease when sharing the road-space with a lorry causing the attention to increase. There is, furthermore, a possibility that the number of lorries per hour is correlated with an unknown factor that lowers the number of accidents. Hours with many lorries may coincide with hours with good road conditions and hours when experienced drivers are driving.

It is obvious that there are shortcomings in analysis indicating that some of the results in this paper can be questioned. For instance, the model for inhomogeneous traffic flow analysis may be considered as a better model in relation to differences between traffic modes. As this model can be assumed to have more explanatory power, no over dispersion is expected in the Poisson regression model. However, in the inhomogeneous traffic flow model, over dispersion becomes

3 The same pattern is given irrespective of the flow of cars and lorries used.

evident, while in the homogeneous traffic flow model, when the traffic flow is treated without consideration to different types of traffic modes, there is generally no over dispersion.

In addition, the accident sample is rather small even though data on accidents for 6.5 years is used. According to Table 3, the average numbers of accidents per traffic flow interval is generally less than 1, indicating low statistical power in especially the analysis of inhomogeneous traffic flows. Random fluctuations may then be a source of error. The maximum hourly flow of lorries is furthermore a fraction of that of cars which may lead to unstable results.4

Consequently, the marginal effect of the flow of different vehicle types may have been miscalculated and exaggerated in the regression analysis due to shortcomings in the data set and in the regression models. However, visual analysis of the inhomogeneous traffic flow data, presented in Figs. 1–8, indicates that the flow of lorries affects the accident frequency differently than the flow of cars. The result implies that important information is lost if no consideration is taken to vehicle type. Literature of today shows limited interest in studying traffic safety models distinguishing between different vehicle types, though. This is a shortcoming since the characteristics of our vehicles and their driving pattern differ in a number of respects, undoubtedly affecting the accident frequency in various ways. Other model structures ought to be applied in future research including, for instance, speed as an explanatory variable.

Acknowledgements

Finical support from the Swedish National Road Administration is gratefully acknowledged. I am grateful for insightful suggestions from Risto Kulmala, and Karin Brundell-Freij. Any errors are my sole responsibility.

Appendix A

Traffic flow data from the stationary counters are applied to calculate traffic flow frequencies for each road section, i.e. set of hours of different traffic volumes. An interval level of 25 vehicles per hour is used when treating the traffic flow as homogeneous. When treating the traffic flow as inhomogeneous, an interval level of 25 cars per hour is used together with an interval level of five lorries per hour. The traffic flow frequency for interval t is denoted by τt. For each road type the number of kilometres, η, is calculated for each section, i.e. ηr kilometres for road type r = I, II, III and IV. The number of kilometres is then multiplied by the traffic flow frequencies for the studied road section. The analysis is based on four road types and each road type is represented on several road sections. If s is the number of sections of road type r, the exposure, H, is calculated as the total number of hours and kilometres studied.

$$H{r,t} = \sum{s=1}^{n} \eta{r,s} \times \tau{t,s}$$

for $t = 1, \dots, p$ and $r = I, \dots, IV$

The traffic flow variables are calculated as the average number of vehicles, q, and the average number of cars and lorries per hour, qc and *q*l, for each traffic flow interval.

The number of accidents throughout the day and in daylight is calculated by aggregation using the same traffic flow interval as when calculating the traffic flow frequencies. These variables are used as dependent variables. When calculating the exposure for daytime traffic, frequencies of traffic flows observed between 6 a.m. and 6 p.m. are used.

References

Ardekani, S., Hauer, E., Jamei, B., 1997. In: Gartner, N. et al. (Eds.), Traffic Impact Models in Traffic Flow Theory—A State-of-the-Art Report, Oak Ridge National Laboratory.

Hall, R.D.,1986. Accidents at Four-arm Single Carriageway Urban Traffic Signals. Transport and Road Research Laboratory Contractor Report, 65.

Hauer, E., Bamfo, J., 1997. In: Proceedings of the ICTCT conference on Two Tools for Finding what Function Links the Dependent Variable to the Explanatory Variables, Lund.

Jovanis, P.P., Chang, H.L., 1986. Modelling the relationship of accidents to miles travelled. Transportation Res. Rec. 1068, 42–51.

Kulmala, R., 1995. Safety at Rural Three- and Four-arm Junctions, VTT-Publications, 233.

Nicholson, A., Wong, Y.D., 1993. Are accidents Poisson distributed? a statistical test. Accid. Anal. Prev. 1, 91–97.

Satterthwaite, P.S., 1981. A Survey of Research into Relationship Between Traffic Accidents and Traffic Volumes, Transport and Road Research Laboratory Contractor Report, 692.

Vickery, W.S., 1968. Automobile Accidents, Tort Law, Externalities and Insurance. Journal of Law and Contemporary problems 33, 464–484. Vitaliano, D.F., Held, J., 1991. Road Accident External Effects: an

Empirical Assessment. Applied Economics 23, 373–378. Winslott, L., 1998. The external costs of traffic accidents: an empirical analysis of the traffic flow. Department of Economics Working Paper Series, 1.

4 In this type of data one could also expect a correlation between the flow of cars and lorries, leading to spurious effects in the regression analysis. No correlation can be detected, however. Even though high flows of lorries are more frequently found at high flows of cars there are also hours in the data set with high flows of lorries and low flows of cars.