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-075    Date:  March 2018
Publication Number: FHWA-HRT-17-075
Date: March 2018

 

Safety Evaluation of Profiled Thermoplastic Pavement Markings

Chapter 3. Methodology

The empirical Bayes (EB) methodology for observational before–after studies was used for the evaluation conducted in this study. This methodology is considered rigorous in that it accounts for regression-to-the-mean using a reference group of similar but untreated sites. In the process, safety performance functions (SPFs) were also used to do the following:

The methodology derived and documented in detail by Hauer is only summarized here. It also provides a foundation for developing guidelines for estimating the likely safety consequences of a contemplated strategy.(3) The SPFs for roadways without profiled thermoplastic pavement markings can be used with observed crash histories to estimate the number of crashes without treatment, and the crash modification factors (CMFs) developed can be applied to this number to estimate the number of crashes with treatment.

In the EB approach, the estimated change in safety for a given crash type at a site is given by the equation in figure 3.

Delta of safety equals lambda minus pi.

Figure 3. Equation. Estimated change in safety.

Where:

λ = Expected number of crashes that would have occurred in the after period without the strategy.
π = Number of reported crashes in the after period.

In estimating λ, the effects of regression-to-the-mean and changes in traffic volume were explicitly accounted for using SPFs, which relate crashes of different types to traffic flow and other relevant factors for each jurisdiction based on untreated sites (i.e., reference sites). Annual SPF multipliers were calibrated to account for temporal effects on safety (e.g., variation in weather, demography, and crash reporting).

In the EB procedure, the SPF is first used to estimate the number of crashes that would be predicted to occur in each year of the before period at reference sites having traffic volumes and other characteristics similar to the one being analyzed. The sum of these annual SPF estimates (P) is then combined with the count of crashes (x) in the before period at a strategy site to obtain an estimate of the predicted number of crashes (m) before strategy. This estimate of m is calculated using the equation in figure 4.

m equals the sum of w times open parenthesis P close parenthesis plus open parenthesis 1 minus w close parenthesis times open parenthesis x close parenthesis.

Figure 4. Equation. EB estimate of expected crashes.

Where w is estimated from the mean and variance of the SPF estimate using the equation in figure 5.

w equals 1 divided by the sum of 1 plus k times P.

Figure 5. Equation. EB weight.

Where k is an overdispersion parameter estimated from the SPF calibration.

In specifying the SPF, a negative binomial distributed error structure is assumed, with k being the overdispersion parameter of this distribution and that is estimated along with the other parameters of the SPF.

A factor is then applied to m to account for the length of the after period and differences in traffic volumes between the before and after periods. This factor is the sum of the annual SPF predictions for the after period divided by P, the sum of these predictions for the before-period. The result, after applying this factor, is an estimate of λ. The procedure also produces an estimate of the variance of λ.

The estimate of λ is then summed over all sites in a strategy group of interest (to obtain λsum) and compared with the count of crashes observed during the after period in that group (πsum). The variance of λ is also summed over all sites in the strategy group.

The index of effectiveness (θ) is estimated using the equation in figure 6.

Theta equals pi subscript sum divided by lambda subscript sum, all divided by 1 plus open parenthesis variance of open parenthesis lambda subscript sum close parenthesis divided by the square of lambda subscript sum close parenthesis.

Figure 6. Equation. Index of effectiveness.(3)

The standard deviation of θ is given by the equation in figure 7.

The standard deviation of open parenthesis theta close parenthesis equals the square root of theta squared times open parenthesis the sum of the quotient of the variance of open parenthesis pi subscript sum close parenthesis divided by the square of pi subscript sum, plus the quotient of the variance of open parenthesis lambda subscript sum close parenthesis divided by the square of lambda subscript sum, close parenthesis, all divided by the square of open parenthesis 1 plus the quotient of the variance of open parenthesis lambda subscript sum close parenthesis divided by the square of lambda subscript sum, close parenthesis.

Figure 7. Equation. Standard deviation of index of effectiveness.(3)

The percent change in crashes is calculated as 100(1 - θ); thus, a value of θ = 0.70 with a standard error (SE) of 0.12 indicates a 30-percent reduction in crashes with an SE of 12 percent.

 

 

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