U.S. Department of Transportation
Federal Highway Administration
1200 New Jersey Avenue, SE
Washington, DC 20590
202-366-4000


Skip to content
Facebook iconYouTube iconTwitter iconFlickr iconLinkedInInstagram

Federal Highway Administration Research and Technology
Coordinating, Developing, and Delivering Highway Transportation Innovations

 
REPORT
This report is an archived publication and may contain dated technical, contact, and link information
Back to Publication List        
Publication Number:  FHWA-HRT-17-069    Date:  December 2017
Publication Number: FHWA-HRT-17-069
Date: December 2017

 

Safety Evaluation of Edge-Line Rumble Stripes on Rural Two-Lane Horizontal Curves

Chapter 6. Development of Safety Performance Functions

This chapter presents the SPFs developed for each crash type. The SPFs were used in the EB methodology to estimate the safety effectiveness of this strategy.(36) Generalized linear modeling was used to estimate model coefficients assuming a negative binomial error distribution, which is consistent with the state of research in developing these models. In specifying a negative binomial error structure, the dispersion parameter, k, was estimated iteratively from the model and the data. For a given dataset, smaller values of k indicate relatively better models.

Crash Spillover and Migration

Before developing SPFs, the project team analyzed the separate reference groups to identify potential crash migration and spillover effects. An SPF was developed using data from both reference groups in order to develop yearly multipliers for each group. The form of the SPF is provided in figure 8, with parameter estimates presented in table 17.

Crashes per mile per year equals e to the a power times AADT to the b power time L to the c power times e to the power of open parenthesis invradius times d plus right_shoulder times e close parenthesis.

Figure 8. Equation. SPF model form for crash migration or spillover effects.

Where:

AADT = annual average daily traffic volume.
L = segment length (mi).
invradius = inverse of the horizontal curve radius (ft).
right_shoulder = right shoulder width (ft).
a, b, c, d, e = parameters estimated in the SPF calibration process.
k = overdispersion parameter of the model.

Table 17. Reference group SPF for total crashes—parameter estimates.

State

a (SE)

b (SE)

c (SE)

d (SE)

e (SE)

k

Kentucky

−9.063
(0.843)

1.288
(0.111)

0.391
(0.219)

27.738
(37.716)

−0.298
(0.159)

2.183

Ohio

−5.226
(0.442)

0.634
(0.046)

0.390
(0.057)

206.883
(47.521)

−0.069
(0.019)

1.022

 

Table 18 presents the observed crashes versus predicted crashes for each of the two reference groups in Kentucky. Group 1 was the reference group more than 5 mi downstream, and group 2 was the reference group immediately downstream of the treatment sites. Group 1 also included reference sites that were on different roadways than the treatment sites. The yearly factors were the ratio of observed crashes to predicted crashes for the given group within the given year. Spillover and crash migration effects would be apparent if the yearly factors became drastically different between group 1 and group 2 after treatment application (2010 for most sites). These effects would also become apparent if the yearly factor for group 2 increased or decreased markedly after treatment application. However, neither of these scenarios appeared to be the case. Figure 9 provides a graphical representation of the yearly factors from table 18.

Table 18. Observed and predicted crashes for reference groups in Kentucky.

Year

Group 1
Observed

Group 2
Observed

Group 1
Predicted

Group 2
Predicted

Group 1
Factor

Group 2
Factor

2004

5

9

6.072

12.806

0.823

0.703

2005

4

13

6.074

12.817

0.659

1.014

2006

5

14

6.077

12.827

0.823

1.091

2007

8

11

6.103

12.875

1.311

0.854

2008

5

12

6.166

12.928

0.811

0.928

2009

5

12

6.265

13.001

0.798

0.923

2010

7

14

6.368

13.065

1.099

1.072

2011

3

23

5.875

12.515

0.511

1.838

2012

13

16

5.446

12.007

2.387

1.333

 

This line graph shows years from 2003 to 2013 on the x-axis and annual factor from 0 to 3 on the y-axis. Three lines are shown: Group 1, Group 2, and Combined. The Group 1 line shows approximate annual factors of 0.8 for 2004, 0.6 for 2005, 0.8 for 2006, 1.3 for 2007, 0.8 for 2008, 0.8 for 2009, 1.1 for 2010, 0.5 for 2011, and 2.4 for 2012. The Group 2 line shows approximate annual factors of 0.7 for 2004, 1.0 for 2005, 1.1 for 2006, 0.8 for 2007, 0.9 for 2008, 0.9 for 2009, 1.1 for 2010, 1.8 for 2011, and 1.3 for 2012. The Combined line shows approximate annual factors of 0.75 for 2004, 0.9 for 2005, 1.0 for 2006, 1.0 for 2007, 0.9 for 2008, 0.9 for 2009, 1.1 for 2010, 1.4 for 2011, and 1.6 for 2012.

Source: FHWA.

Figure 9. Graph. Yearly multiplier by year for reference groups in Kentucky.

 

The sample sizes for reference sites were too small to make a statistical observation of difference between the groups (i.e., fewer than 20 crashes per year were observed). However, both groups experienced an abnormally large number of total crashes in 2011 or 2012, and group 1 observed an abnormally low number of crashes in 2011. Therefore, data for 2011 and 2012 were removed for the reference sites for SPF development. Before-period data from treatment sites were combined with reference-site data to bolster sample size, and the project team estimated an interaction term for the pretreatment site indicator and AADT to determine whether there was a difference in the effect of traffic volume at treatment sites in the before period versus at the references sites. The results indicated that the reference data and pretreatment data could be combined for SPF estimation.

In addition, because of the lack of data at the reference curves in this dataset, the project team sought an alternative reference group for developing annual factors. The Kentucky reference data from Lyon et al. provided a robust dataset for estimating annual factors for the before and after periods from combined horizontal curves and tangents.(30) This reference set included sections that were eligible for resurfacing but had not yet received resurfacing and had texturing; this was consistent with the treatment group. The reference set also included those identified by Lyon et al. to be “resurfacing effort” sites but excluded sites that had SRSs in the prior condition.(30) In total, the reference group consisted of 401.21 roadway mi. The project team used the reference segments to develop predicted crashes, which were compared with observed crashes to develop annual factors. The after-period factor was 1.034 for 2009 installations, 0.993 for 2010 installations, and 0.982 for 2011 installations. Table 19 presents the observed crashes versus predicted crashes for each of the two reference groups in Ohio. As with Kentucky, group 1 represents the reference group more than 5 mi downstream, and group 2 represents the reference group immediately downstream of the treatment sites. Figure 10 provides a graphical representation of the yearly factors from table 19. The annual factors show that the reference sites, particularly those not adjacent to treatment sites, observed a substantial reduction, beginning in 2011, relative to the before period. This possibility was expected for group 2 but not for group 1.

The project team explored this finding further with statewide data for rural two-lane roadways and for horizontal curves on rural two-lane highways. Consistently, the trends showed an approximate 10- to 15-percent reduction in crashes beginning in 2011. The project team contacted ODOT to gain further insight into this finding to determine why this systematic reduction was observed statewide. ODOT noted that in 2010, the systemic program focused on upgrading signage for horizontal curves statewide. This included upgrading hundreds of curves. The upgraded signage included chevrons, curve ahead signs, and speed advisory signs, among others. Therefore, the project team considered using horizontal tangent sections as a potential reference group for developing annual factors to mitigate the impact of signage upgrades on the overall findings. Spillover effects were observed for short tangents; the annual multipliers for short tangents (i.e., tangents less than 0.5 mi in length) were found to match those of horizontal curves within 1.0 percent for each installation year. Tangents longer than 1.0 mi were used as a reference group for developing annual factors because no spillover effects were observed. The reference segments were used to develop predicted crashes, which were compared with observed crashes to develop annual factors. The after-period factor was 1.040 for 2010 installations, 1.043 for 2011 installations, and 1.038 for 2012 installations.

Table 19. Observed and predicted crashes for reference groups in Ohio.

Year

Group 1
Observed

Group 2
Observed

Group 1
Predicted

Group 2
Predicted

Group 1
Factor

Group 2
Factor

2005

81

49

84.82

46.24

0.954963

1.059689

2006

83

60

85.75

45.87

0.96793

1.308044

2007

72

38

84.47

45.37

0.852374

0.837558

2008

79

66

83.64

43.95

0.944524

1.501706

2009

87

49

81.15

44.11

1.072089

1.110859

2010

90

68

90.33

44.56

0.996347

1.526032

2011

58

52

90.33

44.61

0.64209

1.165658

2012

67

47

88.88

44.13

0.753825

1.065035

2013

86

38

89.71

43.88

0.958645

0.865998

 

This line graph shows years from 2004 to 2014 on the x-axis and annual factor from 0 to 1.8 on the y-axis. Three lines are shown: Group 1, Group 2, and Combined. The Group 1 line shows approximate annual factors of 0.9 for 2005, 0.9 for 2006, 0.8 for 2007, 0.9 for 2008, 1.1 for 2009, 1.0 for 2010, 0.6 for 2011, 0.7 for 2012, and 0.95 for 2013. The Group 2 line shows approximate annual factors of 1.1 for 2005, 1.3 for 2006, 0.8 for 2007, 1.5 for 2008, 1.1 for 2009, 1.5 for 2010, 1.2 for 2011, 1.1 for 2012, and 0.85 for 2013. The Combined line shows approximate annual factors of 1.0 for 2005, 1.1 for 2006, 0.8 for 2007, 1.1 for 2008, 1.1 for 2009, 1.2 for 2010, 0.8 for 2011, 0.85 for 2012, and 0.9 for 2013.

Source: FHWA.

Figure 10. Graph. Yearly multiplier by year for reference groups in Ohio.

 

Safety Performance Functions

Figure 11 shows the form of the SPFs for combined reference groups, which are presented in table 20.

Crashes per mile per year is equal to e to the a power times AADT to the b power time L to the c power times e to the power of open parenthesis invradius times d plus right_shoulder times e plus posted times f close parenthesis.

Figure 11. Equation. SPF model form for Kentucky and Ohio.

Where:

posted = posted speed.
f = parameter estimated in the SPF calibration process.

The Kentucky SPF was estimated using treatment site before period data and 2004 to 2010 reference site data. The Ohio SPFs were estimated using reference site data from 2005 to 2010.

Table 20. SPFs by crash type.

Parameters

Kentucky Total (SE)

Ohio Total (SE)

Ohio ROR (SE)

a

−10.338
(0.573)

−6.962
(0.813)

−6.443
(0.927)

b

1.249
(0.072)

0.697
(0.060)

0.553
(0.066)

c

0.128
(0.100)

0.365
(0.069)

0.433
(0.078)

d

N/A

236.022
(58.139)

282.222
(63.957)

e

N/A

−0.063
(0.024)

−0.099
(0.028)

f

N/A

0.021
(0.011)

0.0319
(0.013)

k

0.738

1.083

1.289

N/A = not applicable.

In addition, the project team considered crash sample size for reference sites in the development of SPFs. In Kentucky, total crashes per year ranged from 65 to 77. Other crash types had smaller sample sizes. The Highway Safety Manual recommends a minimum of 100 crashes per year to produce reliable SPFs.(40) Therefore, an SPF was developed for total crashes, and proportion factors relating other crash types to total crashes were used in place of separate SPFs. In Ohio, total crashes ranged from 110 to 158. There were sufficient crashes to develop separate SPFs for total crashes and ROR crashes. To relate other crash types to total crashes, proportion factors were used in place of separate SPFs. The prediction from the SPF was multiplied by the proportion factor to determine the number of predicted crashes of each specific crash type. The following is a list of crash type proportions for Kentucky:

The following is a list of crash type proportions for Ohio:

In addition, observed crashes and predicted crashes were used to develop annual factors for time-based trends at reference sites. Factors were used as multipliers for predicted crashes at treatment sites in the after period. Factors greater than 1.00 indicate an increase in crashes at reference sites, and factors less than 1.00 indicate a decrease in crashes at reference sites. Table 21 provides the annual factors based on total crashes in Kentucky and Ohio for each installation year.

Table 21. Time trend factors for predicted crashes.

Year

Kentucky
Total

Ohio
Total

2009

1.034

N/A

2010

0.993

1.040

2011

0.982

1.043

2012

N/A

1.038

N/A = not applicable.

 

 

Federal Highway Administration | 1200 New Jersey Avenue, SE | Washington, DC 20590 | 202-366-4000
Turner-Fairbank Highway Research Center | 6300 Georgetown Pike | McLean, VA | 22101