phed {tpm}

Calculate PHED Metric

Description

Calculate the CMAQ Traffic Congestion Measure in accordance with FHWA General Guidance and Step-by-Step Metric Calculation Procedures for National Performance Measures for Congestion, Reliability, and Freight, and CMAQ Traffic Congestion Requires the speed limits for all TMC segments. Other parameters are optional and (with defaults supplied from FHWA's guidance). Uses the same travel time readings file as score(..., metric = "LOTTR"). Outputs annual hours of delay.

Usage

phed(
  travel_time_readings,
  tmc_identification,
  speed_limits,
  urban_code = NA,
  pm_peak = 3,
  avo_cars = 1.7,
  avo_trucks = 1,
  avo_buses = 10.7,
  moy_factor = moy_factor_default,
  dow_factor = dow_factor_default,
  hod_profile = hod_profile_default,
  population = NA
)

Arguments

travel_time_readings

path to readings CSV with 15-minute travel time observations for all vehicles exported from RITIS.

tmc_identification

Path to TMC_Identification.csv file provided by RITIS with travel time download.

speed_limits

a data.frame-like object with speed limits for all TMCs with format tmc,speed_limit

urban_code

optional vector of one (or more) urban_code values. if specified, filters TMCs to the relevant urban_code

pm_peak

3 or 4. Peak Period is defined as weekdays from 6 am to 10 am and either 3 pm to 7 pm (3) or 4 pm to 8 pm (4)

avo_cars

Average vehicle occupancy for passenger vehicles

avo_trucks

Average vehicle occupancy for freight trucks

avo_buses

Average vehicle occupancy for buses

moy_factor

Month of year traffic adjustment factors for freeways and non-freeway facilities in format

month,freeway,non_freeway
1, 0.99, 0.98
... 
12, 1.01, 1.00

dow_factor

Day of week adjustment factors. Monday (2) through Friday (6). Format:

day,freeway,non_freeway
2, 1.05, 1.05
...
6, 1.1, 1.1

hod_profile

Hourly traffic percentages for peak hours. Non-directional. Format:

hour,freeway,non_freeway
6,0.045,0.050
...
18,0.052,0.048

population

Optional population value. If provided, function will print PHED value

Value

Annual hours of peak hour excessive delay per capita

Examples

## Not run: 
phed(travel_time_readings = "npmrds/all_vehicles_2021/Readings.csv",
     tmc_identification = "npmrds/all_vehicles_2021/TMC_Identification.csv", 
     speed_limits = fread("birmingham_tmc_speed_limits.csv"), 
     urban_code = 7786,
     pm_peak = 3, 
     output_file_name = "phed_2021.csv",
     avo_cars = 1.62, avo_trucks = 1.0, avo_buses = 5.1,
     moy_factor = fread("birmingham_moy_factors.csv"), 
     dow_factor = fread("birmingham_dow_factors.csv"), 
     hod_profile = fread("birmingham_hod_profile.csv"),
     population = 752898)

## End(Not run)