Speed Distribution by Facility Types

NPMRDS Data

Data

library(data.table)
library(ggplot2)
library(scales)
library(lubridate)
library(DT)

setwd("C:/Users/subas/Syncplicity/MyProjects_IMP/TxDOT_07035/7035_TMC_All")
tmc_dat= fread("7035_TMC.csv")
dim(tmc_dat)
## [1] 17962560        7
tmc_dat= subset(tmc_dat, speed <121)
dim(tmc_dat)
## [1] 11054138        7
summary(tmc_dat)
##    tmc_code         measurement_tstamp     speed      average_speed  
##  Length:11054138    Length:11054138    Min.   : 3.0   Min.   :22.00  
##  Class :character   Class :character   1st Qu.:56.0   1st Qu.:53.00  
##  Mode  :character   Mode  :character   Median :62.0   Median :60.00  
##                                        Mean   :59.9   Mean   :58.02  
##                                        3rd Qu.:66.0   3rd Qu.:64.00  
##                                        Max.   :99.0   Max.   :76.00  
##                                                       NA's   :1099   
##  reference_speed travel_time_minutes data_density      
##  Min.   :58.00   Min.   :  0.060     Length:11054138   
##  1st Qu.:66.00   1st Qu.:  2.780     Class :character  
##  Median :70.00   Median :  4.710     Mode  :character  
##  Mean   :69.54   Mean   :  4.736                       
##  3rd Qu.:72.00   3rd Qu.:  6.320                       
##  Max.   :78.00   Max.   :215.980                       
## 
datatable(head(tmc_dat, 20000))

4M

tmc_dat1= subset(tmc_dat, tmc_code=="112-08977")

tmc_dat1$Time <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%H:%M")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%Y-%m-%d")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date1 <- as.Date(tmc_dat1$Date)
tmc_dat1$Day <- wday(tmc_dat1$Date1, label = TRUE)


ggplot(tmc_dat1, aes(x=Time1, y=speed, colour=Day)) + 
  geom_point(size=.5, alpha = 0.1) + 
  scale_x_datetime(breaks=date_breaks("2 hour"), labels=date_format("%H:%M")) +
  facet_grid(Day ~.) +geom_smooth(col="black")+
  theme(legend.position="none") +theme_bw()

4T

tmc_dat1= subset(tmc_dat, tmc_code=="111-08371")

tmc_dat1$Time <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%H:%M")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%Y-%m-%d")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date1 <- as.Date(tmc_dat1$Date)
tmc_dat1$Day <- wday(tmc_dat1$Date1, label = TRUE)



ggplot(tmc_dat1, aes(x=Time1, y=speed, colour=Day)) + 
  geom_point(size=.5, alpha = 0.1) + 
  scale_x_datetime(breaks=date_breaks("2 hour"), labels=date_format("%H:%M")) +
  facet_grid(Day ~.) +geom_smooth(col="black")+
  theme(legend.position="none") +theme_bw()

4D

tmc_dat1= subset(tmc_dat, tmc_code=="112+07700")

tmc_dat1$Time <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%H:%M")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%Y-%m-%d")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date1 <- as.Date(tmc_dat1$Date)
tmc_dat1$Day <- wday(tmc_dat1$Date1, label = TRUE)



ggplot(tmc_dat1, aes(x=Time1, y=speed, colour=Day)) + 
  geom_point(size=.5, alpha = 0.1) + 
  scale_x_datetime(breaks=date_breaks("2 hour"), labels=date_format("%H:%M")) +
  facet_grid(Day ~.) +geom_smooth(col="black")+
  theme(legend.position="none") +theme_bw()

2SL

tmc_dat1= subset(tmc_dat, tmc_code=="112-08496")

tmc_dat1$Time <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%H:%M")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date <- format(as.POSIXct(strptime(tmc_dat1$measurement_tstamp, "%Y-%m-%d  %H:%M:%S",tz="")) ,format = "%Y-%m-%d")
tmc_dat1$Time1 <- as.POSIXct(tmc_dat1$Time, format = "%H:%M")
tmc_dat1$Date1 <- as.Date(tmc_dat1$Date)
tmc_dat1$Day <- wday(tmc_dat1$Date1, label = TRUE)



ggplot(tmc_dat1, aes(x=Time1, y=speed, colour=Day)) + 
  geom_point(size=.5, alpha = 0.1) + 
  scale_x_datetime(breaks=date_breaks("2 hour"), labels=date_format("%H:%M")) +
  facet_grid(Day ~.) +geom_smooth(col="black")+
  theme(legend.position="none") +theme_bw()