Analysis for all phases

getwd()
## [1] "C:/Users/xuw12/OneDrive - Texas State University/NJDOT Paint Marking/0000_Final/Analysis"
setwd("C:/Users/xuw12/OneDrive - Texas State University/NJDOT Paint Marking/0000_Final/Analysis")

##Start
library(readxl)
library(skimr)
library(ggplot2)
library(ggstatsplot)

dat01= read_xlsx("MasterData_V04.xlsx")
dim(dat01)
## [1] 6780   25
names(dat01)
##  [1] "Id"           "SourceFile"   "RL"           "Max"          "Min"         
##  [6] "Std"          "Width"        "DC"           "RRPM"         "Date"        
## [11] "Driver"       "Length"       "Avg_Length"   "Time"         "Driven"      
## [16] "Route"        "Direction"    "Phase"        "Condition"    "Type"        
## [21] "RoutewithDir" "Color"        "LineType"     "MaterialType" "ftZERO"
skim(dat01)
Data summary
Name dat01
Number of rows 6780
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 40 0 6780 0
SourceFile 0 1 16 36 0 191 0
RL 0 1 1 4 0 628 0
Max 0 1 2 4 0 161 0
Width 0 1 1 18 0 366 0
DC 0 1 1 18 0 702 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 11 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 6780 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 240 0.96 53.07 30.23 2 40 43 52 282 ▇▂▁▁▁
Std 374 0.94 70.58 84.59 2 41 53 83 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 12:43:07 6113 
table(dat01$Site)
## < table of extent 0 >

Analysis by Route

Route GSP

GSP= subset(dat01, Route=="GSP")
dim(GSP)
## [1] 838  25
skim(GSP)
Data summary
Name GSP
Number of rows 838
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 24 0 838 0
SourceFile 0 1 16 19 0 4 0
RL 0 1 1 3 0 332 0
Max 0 1 3 3 0 4 0
Width 0 1 1 18 0 187 0
DC 0 1 1 18 0 287 0
RRPM 0 1 1 2 0 13 0
Driver 0 1 12 12 0 1 0
Length 0 1 3 3 0 1 0
Driven 0 1 1 18 0 212 0
Route 0 1 3 3 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 2 0
Condition 0 1 4 4 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 6 6 0 2 0
Color 0 1 6 6 0 1 0
LineType 0 1 4 4 0 1 0
MaterialType 0 1 13 13 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 838 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 42.73 3.70 40 40 40 42 49 ▇▃▁▁▃
Std 0 1 97.49 17.02 81 83 103 103 123 ▇▁▅▁▃
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2023-12-04 00:00:00 2
Time 0 1 1899-12-31 11:55:57 1899-12-31 12:44:52 1899-12-31 12:23:59 779 
GSP$`RL`= as.integer(GSP$`RL`)
ggbetweenstats(
  data = GSP,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route GSP Bridge

GSPBridge= subset(dat01, Route=="GSP Bridge")
dim(GSPBridge)
## [1] 539  25
skim(GSPBridge)
Data summary
Name GSPBridge
Number of rows 539
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 20 40 0 539 0
SourceFile 0 1 17 36 0 10 0
RL 0 1 2 4 0 272 0
Max 0 1 3 4 0 10 0
Width 0 1 1 18 0 172 0
DC 0 1 1 18 0 265 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 4 0
Driven 0 1 1 18 0 390 0
Route 0 1 10 10 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 2 0
Condition 0 1 4 4 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 13 13 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 22 22 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 539 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 48.74 14.18 42 42 43 43 103 ▇▁▁▁▁
Std 0 1 135.22 85.86 71 71 92 195 345 ▇▁▁▂▁
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2024-09-13 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 10:09:47 1899-12-31 13:50:27 1899-12-31 11:37:56 539 
GSPBridge$`RL`= as.integer(GSPBridge$`RL`)
ggbetweenstats(data=GSPBridge, x= Phase, y=`RL`,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route I 80

I80= subset(dat01, Route=="I 80")
dim(I80)
## [1] 1039   25
skim(I80)
Data summary
Name I80
Number of rows 1039
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 21 35 0 1039 0
SourceFile 0 1 19 31 0 18 0
RL 0 1 2 3 0 400 0
Max 0 1 2 4 0 16 0
Width 0 1 1 18 0 195 0
DC 0 1 1 18 0 258 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 63 0
Route 0 1 4 4 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 12 12 0 1 0
RoutewithDir 0 1 7 7 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 13 13 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 1039 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 107 0.90 59.66 29.98 5 44 49 73 135 ▁▇▃▁▂
Std 59 0.94 46.66 16.19 7 37 49 53 72 ▁▃▅▇▆
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 09:20:10 1899-12-31 15:46:55 1899-12-31 13:33:41 1039 
I80$`RL`= as.integer(I80$`RL`)
ggbetweenstats(data=I80, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 15

sh15= subset(dat01, Route=="SH 15")
dim(sh15)
## [1] 547  25
skim(sh15)
Data summary
Name sh15
Number of rows 547
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 34 0 547 0
SourceFile 0 1 20 30 0 18 0
RL 0 1 1 3 0 266 0
Max 0 1 3 3 0 18 0
Width 0 1 1 18 0 152 0
DC 0 1 1 18 0 232 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 4 0
Driven 0 1 1 18 0 71 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 547 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 20 0.96 45.03 7.48 40 40 41 45 66 ▇▁▁▁▁
Std 23 0.96 111.32 253.17 2 41 47 67 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 10:19:28 1899-12-31 12:24:49 1899-12-31 11:11:27 531 
sh15$`RL`= as.integer(sh15$`RL`)
ggbetweenstats(data=sh15, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 33

SH33= subset(dat01, Route=="SH 33")
dim(SH33)
## [1] 304  25
skim(SH33)
Data summary
Name SH33
Number of rows 304
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 304 0
SourceFile 0 1 20 32 0 12 0
RL 0 1 2 3 0 164 0
Max 0 1 2 3 0 12 0
Width 0 1 1 18 0 109 0
DC 0 1 1 18 0 238 0
RRPM 0 1 1 2 0 16 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 3 0
Driven 0 1 1 18 0 28 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 304 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 47.87 12.71 40 41 43 50 87 ▇▂▁▁▁
Std 78 0.74 47.01 7.55 31 42 50 52 58 ▃▁▅▇▂
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 14:43:59 1899-12-31 15:36:47 1899-12-31 15:11:53 296 
SH33$`RL`= as.integer(SH33$`RL`)
ggbetweenstats(data=SH33, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 37

SH37= subset(dat01, Route=="SH 37")
dim(SH37)
## [1] 411  25
skim(SH37)
Data summary
Name SH37
Number of rows 411
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 411 0
SourceFile 0 1 20 32 0 18 0
RL 0 1 2 3 0 224 0
Max 0 1 2 3 0 17 0
Width 0 1 1 18 0 130 0
DC 0 1 1 18 0 245 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 26 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 411 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 85.60 56.59 4 51 84 103 282 ▇▆▂▁▁
Std 0 1 52.47 24.91 7 42 52 65 99 ▂▂▇▂▃
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 13:25:07 1899-12-31 14:27:15 1899-12-31 14:04:58 401 
SH37$`RL`= as.integer(SH37$`RL`)
ggbetweenstats(data=SH37, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 50

SH50= subset(dat01, Route=="SH 50")
dim(SH50)
## [1] 274  25
skim(SH50)
Data summary
Name SH50
Number of rows 274
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 274 0
SourceFile 0 1 20 32 0 12 0
RL 0 1 2 3 0 194 0
Max 0 1 2 3 0 12 0
Width 0 1 1 18 0 113 0
DC 0 1 1 18 0 149 0
RRPM 0 1 1 2 0 15 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 24 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 274 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 23 0.92 49.76 24.9 4 40 42 83 87 ▁▂▇▁▅
Std 0 1.00 35.54 10.1 17 31 34 43 53 ▅▇▆▅▇
Avg_Length 0 1.00 528.00 0.0 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 09:17:41 1899-12-31 16:02:45 1899-12-31 10:15:17 269 
SH50$`RL`= as.integer(SH50$`RL`)
ggbetweenstats(data=SH50, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 53

SH53= subset(dat01, Route=="SH 53")
dim(SH53)
## [1] 375  25
skim(SH53)
Data summary
Name SH53
Number of rows 375
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 375 0
SourceFile 0 1 20 32 0 10 0
RL 0 1 2 3 0 241 0
Max 0 1 3 3 0 10 0
Width 0 1 1 18 0 131 0
DC 0 1 1 18 0 196 0
RRPM 0 1 1 2 0 20 0
Driver 0 1 12 12 0 1 0
Length 0 1 3 3 0 1 0
Driven 0 1 1 18 0 38 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 18 18 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 375 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 38 0.9 42.81 3.85 40 41 41 43.0 53 ▇▃▁▁▂
Std 75 0.8 76.90 25.41 31 61 91 99.5 107 ▂▃▅▃▇
Avg_Length 0 1.0 528.00 0.00 528 528 528 528.0 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 09:45:59 1899-12-31 12:56:46 1899-12-31 12:23:47 364 
SH53$`RL`= as.integer(SH53$`RL`)
ggbetweenstats(data=SH53, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 54

SH54= subset(dat01, Route=="SH 54")
dim(SH54)
## [1] 220  25
skim(SH54)
Data summary
Name SH54
Number of rows 220
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 220 0
SourceFile 0 1 20 32 0 12 0
RL 0 1 2 3 0 156 0
Max 0 1 3 3 0 12 0
Width 0 1 1 18 0 115 0
DC 0 1 1 18 0 152 0
RRPM 0 1 1 2 0 19 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 21 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 220 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.0 58.12 36.99 2 41 52 80 127 ▃▇▂▃▃
Std 21 0.9 58.16 15.23 40 50 53 65 97 ▇▇▅▁▂
Avg_Length 0 1.0 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 14:20:53 1899-12-31 17:03:38 1899-12-31 15:47:20 217 
SH54$`RL`= as.integer(SH54$`RL`)
ggbetweenstats(data=SH54, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 72

SH72= subset(dat01, Route=="SH 72")
dim(SH72)
## [1] 225  25
skim(SH72)
Data summary
Name SH72
Number of rows 225
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 34 0 225 0
SourceFile 0 1 20 30 0 12 0
RL 0 1 2 3 0 141 0
Max 0 1 2 3 0 12 0
Width 0 1 1 18 0 75 0
DC 0 1 1 18 0 147 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 20 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 225 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 18 0.92 98.44 76.78 2 40 98 168 252 ▇▂▃▃▂
Std 0 1.00 40.79 14.61 2 36 44 49 59 ▂▁▃▆▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 12:07:57 1899-12-31 13:27:13 1899-12-31 13:08:36 223 
SH72$`RL`= as.integer(SH72$`RL`)
ggbetweenstats(data=SH72, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 73

SH73= subset(dat01, Route=="SH 73")
dim(SH73)
## [1] 340  25
skim(SH73)
Data summary
Name SH73
Number of rows 340
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 340 0
SourceFile 0 1 20 32 0 12 0
RL 0 1 1 3 0 171 0
Max 0 1 2 3 0 12 0
Width 0 1 1 18 0 149 0
DC 0 1 1 18 0 193 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 30 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 340 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 40.17 0.38 40 40 40 40 41 ▇▁▁▁▂
Std 25 0.93 71.18 27.97 36 47 53 103 105 ▇▃▁▂▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 14:27:24 1899-12-31 16:57:52 1899-12-31 14:59:51 332 
SH73$`RL`= as.integer(SH73$`RL`)
ggbetweenstats(data=SH73, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 130

US130= subset(dat01, Route=="US 130")
dim(US130)
## [1] 449  25
skim(US130)
Data summary
Name US130
Number of rows 449
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 35 0 449 0
SourceFile 0 1 21 31 0 18 0
RL 0 1 1 3 0 267 0
Max 0 1 2 3 0 18 0
Width 0 1 1 18 0 166 0
DC 0 1 1 18 0 312 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 28 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 9 9 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 449 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 52.26 26.28 40 41 43 50 163 ▇▁▁▁▁
Std 25 0.94 73.42 39.76 28 40 64 95 164 ▇▅▃▁▂
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 10:56:04 1899-12-31 18:26:01 1899-12-31 16:09:13 444 
US130$`RL`= as.integer(US130$`RL`)
ggbetweenstats(data=US130, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 206

US206= subset(dat01, Route=="US 206")
dim(US206)
## [1] 548  25
skim(US206)
Data summary
Name US206
Number of rows 548
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 37 0 548 0
SourceFile 0 1 21 33 0 11 0
RL 0 1 1 3 0 235 0
Max 0 1 2 3 0 11 0
Width 0 1 1 18 0 126 0
DC 0 1 1 18 0 182 0
RRPM 0 1 1 2 0 21 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 3 0 3 0
Driven 0 1 1 18 0 91 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 9 9 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 548 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 43.67 7.56 40 40 40 41 62 ▇▁▁▁▂
Std 43 0.92 36.07 13.07 21 26 33 41 72 ▇▆▂▁▂
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 10:59:07 1899-12-31 14:45:28 1899-12-31 13:02:37 542 
US206$`RL`= as.integer(US206$`RL`)
ggbetweenstats(data=US206, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 322

US322= subset(dat01, Route=="US 322")
dim(US322)
## [1] 285  25
skim(US322)
Data summary
Name US322
Number of rows 285
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 30 0 285 0
SourceFile 0 1 21 26 0 12 0
RL 0 1 1 3 0 195 0
Max 0 1 2 3 0 12 0
Width 0 1 1 18 0 100 0
DC 0 1 1 18 0 167 0
RRPM 0 1 1 2 0 18 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 28 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 9 9 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 285 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 47.76 26.32 4 40 42 51 122 ▁▇▁▁▁
Std 25 0.91 51.37 18.72 7 43 57 63 76 ▂▁▇▇▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-12 00:00:00 3
Time 0 1 1899-12-31 11:45:51 1899-12-31 17:05:06 1899-12-31 13:33:25 283 
US322$`RL`= as.integer(US322$`RL`)
ggbetweenstats(data=US322, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 40

US40= subset(dat01, Route=="US 40")
dim(US40)
## [1] 386  25
skim(US40)
Data summary
Name US40
Number of rows 386
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 386 0
SourceFile 0 1 20 32 0 12 0
RL 0 1 3 3 0 226 0
Max 0 1 3 3 0 12 0
Width 0 1 1 18 0 137 0
DC 0 1 1 18 0 171 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 35 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 8 8 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 386 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 34 0.91 60.20 23.74 40 42 48 80 107 ▇▁▂▁▂
Std 0 1.00 58.56 27.00 7 42 45 89 94 ▂▃▆▃▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 09:27:59 1899-12-31 15:53:18 1899-12-31 10:18:06 383 
US40$`RL`= as.integer(US40$`RL`)
ggbetweenstats(data=US40, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by route with direction

Route GSP Bridge_NB

GSPBridge_NB= subset(dat01, RoutewithDir=="GSP Bridge_NB")
dim(GSPBridge_NB)
## [1] 242  25
skim(GSPBridge_NB)
Data summary
Name GSPBridge_NB
Number of rows 242
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 20 40 0 242 0
SourceFile 0 1 17 36 0 5 0
RL 0 1 2 4 0 165 0
Max 0 1 3 4 0 5 0
Width 0 1 1 18 0 132 0
DC 0 1 1 18 0 160 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 3 0 3 0
Driven 0 1 1 18 0 182 0
Route 0 1 10 10 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 2 0
Condition 0 1 4 4 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 13 13 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 22 22 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 242 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 52.90 18.62 43 43 43 49 103 ▇▁▁▁▁
Std 0 1 154.07 87.14 92 92 92 233 345 ▇▁▁▂▁
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2024-09-13 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 10:19:44 1899-12-31 13:50:27 1899-12-31 11:17:54 242 
GSPBridge_NB$`RL`= as.integer(GSPBridge_NB$`RL`)
ggbetweenstats(data=GSPBridge_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route GSP Bridge_SB

GSPBridge_SB= subset(dat01, RoutewithDir=="GSP Bridge_SB")
dim(GSPBridge_SB)
## [1] 297  25
skim(GSPBridge_SB)
Data summary
Name GSPBridge_SB
Number of rows 297
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 20 40 0 297 0
SourceFile 0 1 17 36 0 5 0
RL 0 1 2 4 0 184 0
Max 0 1 3 4 0 5 0
Width 0 1 1 18 0 113 0
DC 0 1 1 18 0 157 0
RRPM 0 1 1 1 0 9 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 3 0
Driven 0 1 1 18 0 234 0
Route 0 1 10 10 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 2 0
Condition 0 1 4 4 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 13 13 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 22 22 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 297 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 45.35 7.61 42 42 42 43 66 ▇▁▁▁▁
Std 0 1 119.87 81.79 71 71 71 164 332 ▇▁▁▁▁
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2024-09-13 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 10:09:47 1899-12-31 13:42:25 1899-12-31 11:44:10 297 
GSPBridge_SB$`RL`= as.integer(GSPBridge_SB$`RL`)
ggbetweenstats(data=GSPBridge_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route GSP_NB

GSP_NB= subset(dat01, RoutewithDir=="GSP_NB")
dim(GSP_NB)
## [1] 415  25
skim(GSP_NB)
Data summary
Name GSP_NB
Number of rows 415
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 24 0 415 0
SourceFile 0 1 16 19 0 2 0
RL 0 1 1 3 0 210 0
Max 0 1 3 3 0 2 0
Width 0 1 2 18 0 175 0
DC 0 1 1 18 0 195 0
RRPM 0 1 1 2 0 9 0
Driver 0 1 12 12 0 1 0
Length 0 1 3 3 0 1 0
Driven 0 1 1 18 0 208 0
Route 0 1 3 3 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 2 0
Condition 0 1 4 4 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 6 6 0 1 0
Color 0 1 6 6 0 1 0
LineType 0 1 4 4 0 1 0
MaterialType 0 1 13 13 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 415 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 45.51 3.50 42 42 49 49 49 ▇▁▁▁▇
Std 0 1 102.05 21.03 81 81 123 123 123 ▇▁▁▁▇
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2023-12-04 00:00:00 2
Time 0 1 1899-12-31 11:55:57 1899-12-31 12:17:41 1899-12-31 12:07:02 383 
GSP_NB$`RL`= as.integer(GSP_NB$`RL`)
ggbetweenstats(data=GSP_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route GSP_SB

GSP_SB= subset(dat01, RoutewithDir=="GSP_SB")
dim(GSP_SB)
## [1] 423  25
skim(GSP_SB)
Data summary
Name GSP_SB
Number of rows 423
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 24 0 423 0
SourceFile 0 1 16 19 0 2 0
RL 0 1 1 3 0 252 0
Max 0 1 3 3 0 2 0
Width 0 1 1 18 0 109 0
DC 0 1 1 18 0 170 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 3 3 0 1 0
Driven 0 1 1 18 0 212 0
Route 0 1 3 3 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 2 0
Condition 0 1 4 4 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 6 6 0 1 0
Color 0 1 6 6 0 1 0
LineType 0 1 4 4 0 1 0
MaterialType 0 1 13 13 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 423 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 40.00 0.00 40 40 40 40 40 ▁▁▇▁▁
Std 0 1 93.02 10.01 83 83 103 103 103 ▇▁▁▁▇
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2023-12-04 00:00:00 2
Time 0 1 1899-12-31 12:23:42 1899-12-31 12:44:52 1899-12-31 12:33:57 396 
GSP_SB$`RL`= as.integer(GSP_SB$`RL`)
ggbetweenstats(data=GSP_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route I 80_EB

I80_EB= subset(dat01, RoutewithDir=="I 80_EB")
dim(I80_EB)
## [1] 505  25
skim(I80_EB)
Data summary
Name I80_EB
Number of rows 505
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 21 35 0 505 0
SourceFile 0 1 19 31 0 9 0
RL 0 1 2 3 0 273 0
Max 0 1 2 4 0 8 0
Width 0 1 1 18 0 178 0
DC 0 1 1 18 0 182 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 60 0
Route 0 1 4 4 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 12 12 0 1 0
RoutewithDir 0 1 7 7 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 13 13 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 505 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 54 0.89 50.29 6.50 41 48 48 53.0 62 ▅▇▂▁▅
Std 59 0.88 50.33 12.83 31 49 49 63.5 70 ▃▁▇▁▃
Avg_Length 0 1.00 528.00 0.00 528 528 528 528.0 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 09:20:10 1899-12-31 15:46:55 1899-12-31 13:31:34 505 
I80_EB$`RL`= as.integer(I80_EB$`RL`)
ggbetweenstats(data=I80_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route I 80_WB

I80_WB= subset(dat01, RoutewithDir=="I 80_WB")
dim(I80_WB)
## [1] 534  25
skim(I80_WB)
Data summary
Name I80_WB
Number of rows 534
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 21 35 0 534 0
SourceFile 0 1 19 31 0 9 0
RL 0 1 2 3 0 298 0
Max 0 1 2 3 0 9 0
Width 0 1 1 4 0 134 0
DC 0 1 1 18 0 216 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 63 0
Route 0 1 4 4 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 12 12 0 1 0
RoutewithDir 0 1 7 7 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 13 13 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 534 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 53 0.9 68.46 39.29 5 40 55 114 135 ▂▇▅▁▅
Std 0 1.0 43.59 17.98 7 37 42 50 72 ▂▃▇▅▆
Avg_Length 0 1.0 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 09:27:12 1899-12-31 15:38:29 1899-12-31 13:40:01 534 
I80_WB$`RL`= as.integer(I80_WB$`RL`)
ggbetweenstats(data=I80_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 15_NB

SH15_NB= subset(dat01, RoutewithDir=="SH 15_NB")
dim(SH15_NB)
## [1] 284  25
skim(SH15_NB)
Data summary
Name SH15_NB
Number of rows 284
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 34 0 284 0
SourceFile 0 1 20 30 0 9 0
RL 0 1 1 3 0 181 0
Max 0 1 3 3 0 9 0
Width 0 1 1 18 0 120 0
DC 0 1 1 18 0 173 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 3 0
Driven 0 1 1 18 0 64 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 284 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 43.47 3.67 40 40 42 45 50 ▇▃▂▁▃
Std 0 1 164.02 334.50 37 46 61 61 1193 ▇▁▁▁▁
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 10:19:28 1899-12-31 12:13:43 1899-12-31 11:07:15 278 
SH15_NB$`RL`= as.integer(SH15_NB$`RL`)
ggbetweenstats(data=SH15_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 15_SB

SH15_SB= subset(dat01, RoutewithDir=="SH 15_SB")
dim(SH15_SB)
## [1] 263  25
skim(SH15_SB)
Data summary
Name SH15_SB
Number of rows 263
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 34 0 263 0
SourceFile 0 1 20 30 0 9 0
RL 0 1 2 3 0 160 0
Max 0 1 3 3 0 9 0
Width 0 1 2 18 0 116 0
DC 0 1 1 18 0 122 0
RRPM 0 1 1 2 0 11 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 3 0
Driven 0 1 1 18 0 71 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 263 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 20 0.92 46.86 9.98 40 40 40 60 66 ▇▁▁▂▁
Std 23 0.91 48.96 23.83 2 33 41 67 89 ▂▆▇▇▃
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 10:24:43 1899-12-31 12:24:49 1899-12-31 11:12:19 261 
SH15_SB$`RL`= as.integer(SH15_SB$`RL`)
ggbetweenstats(data=SH15_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 33_EB

SH33_EB= subset(dat01, RoutewithDir=="SH 33_EB")
dim(SH33_EB)
## [1] 158  25
skim(SH33_EB)
Data summary
Name SH33_EB
Number of rows 158
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 158 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 112 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 85 0
DC 0 1 1 18 0 133 0
RRPM 0 1 1 2 0 13 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 27 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 158 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 50.2 15.94 40 42.0 43 50 87 ▇▂▁▁▂
Std 54 0.66 46.5 8.41 36 40.5 42 58 58 ▇▇▁▇▇
Avg_Length 0 1.00 528.0 0.00 528 528.0 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 14:43:59 1899-12-31 15:29:47 1899-12-31 15:05:37 153 
SH33_EB$`RL`= as.integer(SH33_EB$`RL`)
ggbetweenstats(data=SH33_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 33_WB

SH33_WB= subset(dat01, RoutewithDir=="SH 33_WB")
dim(SH33_WB)
## [1] 146  25
skim(SH33_WB)
Data summary
Name SH33_WB
Number of rows 146
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 146 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 99 0
Max 0 1 2 3 0 6 0
Width 0 1 2 18 0 79 0
DC 0 1 1 18 0 112 0
RRPM 0 1 1 2 0 15 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 3 0
Driven 0 1 1 18 0 28 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 146 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 45.36 7.09 40 41 41 53 59 ▇▁▁▁▂
Std 24 0.84 47.44 6.75 31 45 50 52 52 ▂▁▁▃▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 14:52:41 1899-12-31 15:36:47 1899-12-31 15:19:29 143 
SH33_WB$`RL`= as.integer(SH33_WB$`RL`)
ggbetweenstats(data=SH33_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 37_EB

SH37_EB= subset(dat01, RoutewithDir=="SH 37_EB")
dim(SH37_EB)
## [1] 198  25
skim(SH37_EB)
Data summary
Name SH37_EB
Number of rows 198
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 198 0
SourceFile 0 1 20 32 0 9 0
RL 0 1 2 3 0 143 0
Max 0 1 2 3 0 9 0
Width 0 1 1 18 0 106 0
DC 0 1 1 18 0 148 0
RRPM 0 1 1 1 0 9 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 25 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 198 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 88.09 71.50 44 49 52 104 282 ▇▂▁▁▁
Std 0 1 60.33 21.47 31 52 52 65 99 ▃▇▂▁▃
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 13:25:07 1899-12-31 14:23:32 1899-12-31 14:01:09 195 
SH37_EB$`RL`= as.integer(SH37_EB$`RL`)
ggbetweenstats(data=SH37_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 37_WB

SH37_WB= subset(dat01, RoutewithDir=="SH 37_WB")
dim(SH37_WB)
## [1] 213  25
skim(SH37_WB)
Data summary
Name SH37_WB
Number of rows 213
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 213 0
SourceFile 0 1 20 32 0 9 0
RL 0 1 2 3 0 135 0
Max 0 1 2 3 0 9 0
Width 0 1 2 18 0 103 0
DC 0 1 1 18 0 145 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 26 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 213 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 83.29 37.83 4 58 98 103 141 ▂▃▂▇▂
Std 0 1 45.16 25.70 7 42 44 52 92 ▃▁▇▂▂
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 13:29:32 1899-12-31 14:27:15 1899-12-31 14:05:56 209 
SH37_WB$`RL`= as.integer(SH37_WB$`RL`)
ggbetweenstats(data=SH37_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 50_EB

SH50_EB= subset(dat01, RoutewithDir=="SH 50_EB")
dim(SH50_EB)
## [1] 136  25
skim(SH50_EB)
Data summary
Name SH50_EB
Number of rows 136
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 136 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 104 0
Max 0 1 2 3 0 6 0
Width 0 1 1 18 0 75 0
DC 0 1 1 18 0 99 0
RRPM 0 1 1 1 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 24 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 136 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 41.94 25.00 4 27 42 48 87 ▂▃▇▁▃
Std 0 1 33.55 7.69 21 31 32 39 46 ▃▇▃▃▃
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 09:44:36 1899-12-31 16:02:45 1899-12-31 10:42:16 134 
SH50_EB$`RL`= as.integer(SH50_EB$`RL`)
ggbetweenstats(data=SH50_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 50_WB

SH50_WB= subset(dat01, RoutewithDir=="SH 50_WB")
dim(SH50_WB)
## [1] 138  25
skim(SH50_WB)
Data summary
Name SH50_WB
Number of rows 138
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 138 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 110 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 80 0
DC 0 1 1 18 0 91 0
RRPM 0 1 1 2 0 14 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 24 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 138 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 23 0.83 59.00 21.47 40 41 44 83 87 ▇▁▁▁▅
Std 0 1.00 37.51 11.71 17 31 35 47 53 ▅▅▅▅▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 09:17:41 1899-12-31 15:38:42 1899-12-31 10:02:41 135 
SH50_WB$`RL`= as.integer(SH50_WB$`RL`)
ggbetweenstats(data=SH50_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 53_NB

SH53_NB= subset(dat01, RoutewithDir=="SH 53_NB")
dim(SH53_NB)
## [1] 188  25
skim(SH53_NB)
Data summary
Name SH53_NB
Number of rows 188
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 188 0
SourceFile 0 1 20 32 0 5 0
RL 0 1 2 3 0 149 0
Max 0 1 3 3 0 5 0
Width 0 1 2 18 0 93 0
DC 0 1 1 18 0 129 0
RRPM 0 1 1 2 0 18 0
Driver 0 1 12 12 0 1 0
Length 0 1 3 3 0 1 0
Driven 0 1 1 18 0 38 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 18 18 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 188 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 38 0.8 44.05 5.28 40 41 42 50.25 53 ▇▁▁▁▃
Std 37 0.8 70.76 26.48 31 62 91 94.50 98 ▃▁▃▁▇
Avg_Length 0 1.0 528.00 0.00 528 528 528 528.00 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 09:58:08 1899-12-31 12:56:46 1899-12-31 12:29:13 181 
SH53_NB$`RL`= as.integer(SH53_NB$`RL`)
ggbetweenstats(data=SH53_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 53_SB

SH53_SB= subset(dat01, RoutewithDir=="SH 53_SB")
dim(SH53_SB)
## [1] 187  25
skim(SH53_SB)
Data summary
Name SH53_SB
Number of rows 187
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 187 0
SourceFile 0 1 20 32 0 5 0
RL 0 1 2 3 0 148 0
Max 0 1 3 3 0 5 0
Width 0 1 1 18 0 102 0
DC 0 1 1 18 0 131 0
RRPM 0 1 1 2 0 20 0
Driver 0 1 12 12 0 1 0
Length 0 1 3 3 0 1 0
Driven 0 1 1 18 0 38 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 18 18 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 187 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.0 41.82 1.47 40 41 41 43 44 ▃▇▁▅▅
Std 38 0.8 83.13 22.71 58 58 104 107 107 ▇▁▁▁▇
Avg_Length 0 1.0 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 09:45:59 1899-12-31 12:49:41 1899-12-31 12:20:35 183 
SH53_SB$`RL`= as.integer(SH53_SB$`RL`)
ggbetweenstats(data=SH53_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 54_EB

SH54_EB= subset(dat01, RoutewithDir=="SH 54_EB")
dim(SH54_EB)
## [1] 107  25
skim(SH54_EB)
Data summary
Name SH54_EB
Number of rows 107
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 107 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 85 0
Max 0 1 3 3 0 6 0
Width 0 1 2 18 0 75 0
DC 0 1 1 18 0 85 0
RRPM 0 1 1 2 0 17 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 21 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 107 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.0 58.21 34.13 4 42 55 79 112 ▅▃▇▅▅
Std 21 0.8 52.92 10.26 40 44 52 56 69 ▇▁▇▁▅
Avg_Length 0 1.0 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 14:25:00 1899-12-31 17:03:38 1899-12-31 15:52:51 105 
SH54_EB$`RL`= as.integer(SH54_EB$`RL`)
ggbetweenstats(data=SH54_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 54_WB

SH54_WB= subset(dat01, RoutewithDir=="SH 54_WB")
dim(SH54_WB)
## [1] 113  25
skim(SH54_WB)
Data summary
Name SH54_WB
Number of rows 113
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 113 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 90 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 81 0
DC 0 1 1 18 0 88 0
RRPM 0 1 1 2 0 16 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 21 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 113 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 58.04 39.66 2 41 48 80 127 ▃▇▁▃▃
Std 0 1 62.14 17.12 44 50 61 65 97 ▇▆▁▁▃
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 14:20:53 1899-12-31 17:00:23 1899-12-31 15:47:09 112 
SH54_WB$`RL`= as.integer(SH54_WB$`RL`)
ggbetweenstats(data=SH54_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 72_EB

SH72_EB= subset(dat01, RoutewithDir=="SH 72_EB")
dim(SH72_EB)
## [1] 113  25
skim(SH72_EB)
Data summary
Name SH72_EB
Number of rows 113
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 34 0 113 0
SourceFile 0 1 20 30 0 6 0
RL 0 1 2 3 0 91 0
Max 0 1 2 3 0 6 0
Width 0 1 1 18 0 58 0
DC 0 1 1 18 0 92 0
RRPM 0 1 1 1 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 20 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 113 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 18 0.84 94.63 82.16 40 40 43 98 252 ▇▂▁▁▂
Std 0 1.00 46.31 6.65 40 40 45 50 59 ▇▇▃▁▃
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 12:07:57 1899-12-31 13:21:51 1899-12-31 13:08:36 112 
SH72_EB$`RL`= as.integer(SH72_EB$`RL`)
ggbetweenstats(data=SH72_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 72_WB

SH72_WB= subset(dat01, RoutewithDir=="SH 72_WB")
dim(SH72_WB)
## [1] 112  25
skim(SH72_WB)
Data summary
Name SH72_WB
Number of rows 112
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 34 0 112 0
SourceFile 0 1 20 30 0 6 0
RL 0 1 3 3 0 86 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 44 0
DC 0 1 3 18 0 88 0
RRPM 0 1 1 2 0 10 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 20 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 7 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 112 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 101.67 72.12 2 7 132 168 172 ▆▁▁▂▇
Std 0 1 35.22 18.00 2 26 36 49 53 ▃▁▂▃▇
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 12:16:29 1899-12-31 13:27:13 1899-12-31 13:10:53 111 
SH72_WB$`RL`= as.integer(SH72_WB$`RL`)
ggbetweenstats(data=SH72_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 73_NB

SH73_NB= subset(dat01, RoutewithDir=="SH 73_NB")
dim(SH73_NB)
## [1] 169  25
skim(SH73_NB)
Data summary
Name SH73_NB
Number of rows 169
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 169 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 1 3 0 104 0
Max 0 1 2 3 0 6 0
Width 0 1 1 18 0 110 0
DC 0 1 1 18 0 122 0
RRPM 0 1 1 2 0 9 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 30 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 169 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 40.35 0.48 40 40 40 41 41 ▇▁▁▁▅
Std 25 0.85 70.29 28.52 39 48 53 103 105 ▇▅▁▁▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 14:27:24 1899-12-31 16:52:16 1899-12-31 14:54:19 165 
SH73_NB$`RL`= as.integer(SH73_NB$`RL`)
ggbetweenstats(data=SH73_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route SH 73_SB

SH73_SB= subset(dat01, RoutewithDir=="SH 73_SB")
dim(SH73_SB)
## [1] 171  25
skim(SH73_SB)
Data summary
Name SH73_SB
Number of rows 171
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 171 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 2 3 0 118 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 95 0
DC 0 1 1 18 0 127 0
RRPM 0 1 1 1 0 9 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 30 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 171 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 40.00 0.00 40 40 40 40 40 ▁▁▇▁▁
Std 0 1 71.92 27.56 36 47 53 103 103 ▇▃▁▃▇
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 14:34:13 1899-12-31 16:57:52 1899-12-31 15:00:11 167 
SH73_SB$`RL`= as.integer(SH73_SB$`RL`)
ggbetweenstats(data=SH73_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 130_NB

US130_NB= subset(dat01, RoutewithDir=="US 130_NB")
dim(US130_NB)
## [1] 224  25
skim(US130_NB)
Data summary
Name US130_NB
Number of rows 224
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 35 0 224 0
SourceFile 0 1 21 31 0 9 0
RL 0 1 2 3 0 163 0
Max 0 1 2 3 0 9 0
Width 0 1 1 18 0 119 0
DC 0 1 2 18 0 184 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 28 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 9 9 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 224 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 61.39 34.73 42 43 49 62 163 ▇▁▁▁▁
Std 25 0.89 59.56 26.12 28 38 49 83 109 ▇▂▅▃▂
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 11:02:45 1899-12-31 18:26:01 1899-12-31 16:12:28 221 
US130_NB$`RL`= as.integer(US130_NB$`RL`)
ggbetweenstats(data=US130_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 130_SB

US130_SB= subset(dat01, RoutewithDir=="US 130_SB")
dim(US130_SB)
## [1] 225  25
skim(US130_SB)
Data summary
Name US130_SB
Number of rows 225
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 35 0 225 0
SourceFile 0 1 21 31 0 9 0
RL 0 1 1 3 0 162 0
Max 0 1 2 3 0 9 0
Width 0 1 1 18 0 122 0
DC 0 1 1 18 0 168 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 28 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 9 9 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 225 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 43.17 3.77 40 41 41 43 52 ▇▂▁▁▁
Std 0 1 85.68 45.39 30 50 76 123 164 ▇▅▃▂▅
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 10:56:04 1899-12-31 18:19:06 1899-12-31 16:04:21 223 
US130_SB$`RL`= as.integer(US130_SB$`RL`)
ggbetweenstats(data=US130_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 206_NB

US206_NB= subset(dat01, RoutewithDir=="US 206_NB")
dim(US206_NB)
## [1] 319  25
skim(US206_NB)
Data summary
Name US206_NB
Number of rows 319
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 37 0 319 0
SourceFile 0 1 21 33 0 6 0
RL 0 1 1 3 0 171 0
Max 0 1 2 3 0 6 0
Width 0 1 1 18 0 105 0
DC 0 1 1 18 0 152 0
RRPM 0 1 1 2 0 21 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 3 0 3 0
Driven 0 1 1 18 0 91 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 9 9 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 319 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 43.19 6.76 40 40 40 41 59 ▇▁▁▁▂
Std 0 1 37.52 14.98 26 26 31 41 72 ▇▅▁▁▂
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 11:13:08 1899-12-31 14:45:28 1899-12-31 13:24:05 315 
US206_NB$`RL`= as.integer(US206_NB$`RL`)
ggbetweenstats(data=US206_NB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 206_SB

US206_SB= subset(dat01, RoutewithDir=="US 206_SB")
dim(US206_SB)
## [1] 229  25
skim(US206_SB)
Data summary
Name US206_SB
Number of rows 229
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 37 0 229 0
SourceFile 0 1 21 33 0 5 0
RL 0 1 1 3 0 149 0
Max 0 1 3 3 0 5 0
Width 0 1 1 18 0 98 0
DC 0 1 1 18 0 122 0
RRPM 0 1 1 2 0 19 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 3 0 2 0
Driven 0 1 1 18 0 50 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 9 9 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 229 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 44.33 8.52 40 40 40 41 62 ▇▁▁▁▂
Std 43 0.81 33.58 8.33 21 21 33 39 43 ▅▁▅▁▇
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 3
Time 0 1 1899-12-31 10:59:07 1899-12-31 14:03:07 1899-12-31 12:37:07 227 
US206_SB$`RL`= as.integer(US206_SB$`RL`)
ggbetweenstats(data=US206_SB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 322_EB

US322_EB= subset(dat01, RoutewithDir=="US 322_EB")
dim(US322_EB)
## [1] 130  25
skim(US322_EB)
Data summary
Name US322_EB
Number of rows 130
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 30 0 130 0
SourceFile 0 1 21 26 0 6 0
RL 0 1 2 3 0 111 0
Max 0 1 2 3 0 6 0
Width 0 1 1 18 0 73 0
DC 0 1 2 18 0 101 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 27 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 9 9 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 130 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 33.83 15.02 4 40 40 42 43 ▂▁▁▁▇
Std 25 0.81 59.05 11.45 44 48 58 59 76 ▅▁▇▁▅
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-12 00:00:00 3
Time 0 1 1899-12-31 11:45:51 1899-12-31 16:56:55 1899-12-31 13:26:34 129 
US322_EB$`RL`= as.integer(US322_EB$`RL`)
ggbetweenstats(data=US322_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 322_WB

US322_WB= subset(dat01, RoutewithDir=="US 322_WB")
dim(US322_WB)
## [1] 155  25
skim(US322_WB)
Data summary
Name US322_WB
Number of rows 155
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 30 0 155 0
SourceFile 0 1 21 26 0 6 0
RL 0 1 1 3 0 123 0
Max 0 1 3 3 0 6 0
Width 0 1 2 18 0 70 0
DC 0 1 1 18 0 105 0
RRPM 0 1 1 2 0 17 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 28 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 9 9 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 155 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 59.44 28.06 41 41 49 56 122 ▇▁▁▁▂
Std 0 1 46.17 20.81 7 37 43 63 71 ▃▁▇▃▇
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-13 00:00:00 2024-09-12 00:00:00 3
Time 0 1 1899-12-31 11:54:28 1899-12-31 17:05:06 1899-12-31 13:35:02 154 
US322_WB$`RL`= as.integer(US322_WB$`RL`)
ggbetweenstats(data=US322_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 40_EB

US40_EB= subset(dat01, RoutewithDir=="US 40_EB")
dim(US40_EB)
## [1] 192  25
skim(US40_EB)
Data summary
Name US40_EB
Number of rows 192
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 192 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 3 3 0 138 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 98 0
DC 0 1 3 18 0 112 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 35 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 192 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 34 0.82 57.07 16.08 40 41 55.0 70 80 ▇▃▁▃▅
Std 0 1.00 54.25 30.49 7 32 53.5 86 94 ▅▃▅▅▇
Avg_Length 0 1.00 528.00 0.00 528 528 528.0 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 09:34:06 1899-12-31 15:53:18 1899-12-31 10:18:20 191 
US40_EB$`RL`= as.integer(US40_EB$`RL`)
ggbetweenstats(data=US40_EB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Route US 40_WB

US40_WB= subset(dat01, RoutewithDir=="US 40_WB")
dim(US40_WB)
## [1] 194  25
skim(US40_WB)
Data summary
Name US40_WB
Number of rows 194
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 194 0
SourceFile 0 1 20 32 0 6 0
RL 0 1 3 3 0 149 0
Max 0 1 3 3 0 6 0
Width 0 1 1 18 0 100 0
DC 0 1 1 18 0 130 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 1 0
Driven 0 1 1 18 0 35 0
Route 0 1 5 5 0 1 0
Direction 0 1 2 2 0 1 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 13 13 0 1 0
RoutewithDir 0 1 8 8 0 1 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 194 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1 62.74 28.28 42 42 42 101 107 ▇▁▁▁▃
Std 0 1 62.82 22.30 40 42 45 89 93 ▇▁▂▁▆
Avg_Length 0 1 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-04 00:00:00 2025-04-14 00:00:00 2024-09-15 00:00:00 3
Time 0 1 1899-12-31 09:27:59 1899-12-31 15:48:46 1899-12-31 10:12:55 192 
US40_WB$`RL`= as.integer(US40_WB$`RL`)
ggbetweenstats(data=US40_WB, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by marking condition

con= subset(dat01)
dim(con)
## [1] 6780   25
skim(con)
Data summary
Name con
Number of rows 6780
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 40 0 6780 0
SourceFile 0 1 16 36 0 191 0
RL 0 1 1 4 0 628 0
Max 0 1 2 4 0 161 0
Width 0 1 1 18 0 366 0
DC 0 1 1 18 0 702 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 11 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 6780 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 240 0.96 53.07 30.23 2 40 43 52 282 ▇▂▁▁▁
Std 374 0.94 70.58 84.59 2 41 53 83 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 12:43:07 6113 
con$`RL`= as.integer(con$`RL`)
ggbetweenstats(data=con, x= Condition, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by marking type

typ= subset(dat01)
dim(typ)
## [1] 6780   25
skim(typ)
Data summary
Name typ
Number of rows 6780
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 40 0 6780 0
SourceFile 0 1 16 36 0 191 0
RL 0 1 1 4 0 628 0
Max 0 1 2 4 0 161 0
Width 0 1 1 18 0 366 0
DC 0 1 1 18 0 702 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 11 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 6780 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 240 0.96 53.07 30.23 2 40 43 52 282 ▇▂▁▁▁
Std 374 0.94 70.58 84.59 2 41 53 83 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 12:43:07 6113 
typ$`RL`= as.integer(typ$`RL`)
ggbetweenstats(data=typ, x= Type, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by Material Type

MaterialType= subset(dat01)
dim(MaterialType)
## [1] 6780   25
skim(MaterialType)
Data summary
Name MaterialType
Number of rows 6780
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 40 0 6780 0
SourceFile 0 1 16 36 0 191 0
RL 0 1 1 4 0 628 0
Max 0 1 2 4 0 161 0
Width 0 1 1 18 0 366 0
DC 0 1 1 18 0 702 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 11 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 6780 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 240 0.96 53.07 30.23 2 40 43 52 282 ▇▂▁▁▁
Std 374 0.94 70.58 84.59 2 41 53 83 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 12:43:07 6113 
MaterialType$`RL`= as.integer(MaterialType$`RL`)
ggbetweenstats(data=MaterialType, x= MaterialType, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by Phase

phase= subset(dat01)
dim(phase)
## [1] 6780   25
skim(phase)
Data summary
Name phase
Number of rows 6780
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 40 0 6780 0
SourceFile 0 1 16 36 0 191 0
RL 0 1 1 4 0 628 0
Max 0 1 2 4 0 161 0
Width 0 1 1 18 0 366 0
DC 0 1 1 18 0 702 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 11 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 6780 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 240 0.96 53.07 30.23 2 40 43 52 282 ▇▂▁▁▁
Std 374 0.94 70.58 84.59 2 41 53 83 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 12:43:07 6113 
phase$`RL`= as.integer(phase$`RL`)
ggbetweenstats(data=phase, x= Phase, y=`RL`,
               type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by Condition

Good

Good= subset(dat01, Condition=="Good")
dim(Good)
## [1] 5444   25
skim(Good)
Data summary
Name Good
Number of rows 5444
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 40 0 5444 0
SourceFile 0 1 16 36 0 143 0
RL 0 1 1 4 0 609 0
Max 0 1 2 4 0 124 0
Width 0 1 1 18 0 348 0
DC 0 1 1 18 0 554 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 10 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 12 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 4 0 1 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 24 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 5444 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 220 0.96 54.79 32.54 2 40 43 55 282 ▇▂▁▁▁
Std 301 0.94 66.15 45.01 2 41 53 83 345 ▇▅▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 17:05:06 1899-12-31 12:44:26 4962 
Good$`RL`= as.integer(Good$`RL`)
ggbetweenstats(
  data = Good,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Mixed

Mixed= subset(dat01, Condition=="Mixed")
dim(Mixed)
## [1] 449  25
skim(Mixed)
Data summary
Name Mixed
Number of rows 449
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 35 0 449 0
SourceFile 0 1 21 31 0 18 0
RL 0 1 1 3 0 267 0
Max 0 1 2 3 0 18 0
Width 0 1 1 18 0 166 0
DC 0 1 1 18 0 312 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 2 0
Driven 0 1 1 18 0 28 0
Route 0 1 6 6 0 1 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 15 15 0 1 0
RoutewithDir 0 1 9 9 0 2 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 449 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 0 1.00 52.26 26.28 40 41 43 50 163 ▇▁▁▁▁
Std 25 0.94 73.42 39.76 28 40 64 95 164 ▇▅▃▁▂
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-05 00:00:00 2025-04-15 00:00:00 2024-09-13 00:00:00 3
Time 0 1 1899-12-31 10:56:04 1899-12-31 18:26:01 1899-12-31 16:09:13 444 
Mixed$`RL`= as.integer(Mixed$`RL`)
ggbetweenstats(
  data = Mixed,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Older

Older= subset(dat01, Condition=="Older")
dim(Older)
## [1] 887  25
skim(Older)
Data summary
Name Older
Number of rows 887
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 36 0 887 0
SourceFile 0 1 20 32 0 30 0
RL 0 1 1 3 0 306 0
Max 0 1 2 3 0 30 0
Width 0 1 1 18 0 193 0
DC 0 1 1 18 0 313 0
RRPM 0 1 1 2 0 12 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 2 0 4 0
Driven 0 1 1 18 0 71 0
Route 0 1 5 5 0 2 0
Direction 0 1 2 2 0 2 0
Phase 0 1 7 7 0 3 0
Condition 0 1 5 5 0 1 0
Type 0 1 7 15 0 2 0
RoutewithDir 0 1 8 8 0 4 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 6 0 3 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 887 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 20 0.98 43.13 6.30 40 40 40 43 66 ▇▁▁▁▁
Std 48 0.95 96.25 201.68 2 43 53 89 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 6
Time 0 1 1899-12-31 10:19:28 1899-12-31 16:57:52 1899-12-31 11:54:16 863 
Older$`RL`= as.integer(Older$`RL`)
ggbetweenstats(
  data = Older,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by Color

Yellow

Yellow= subset(dat01, Color=="Yellow")
dim(Yellow)
## [1] 3622   25
skim(Yellow)
Data summary
Name Yellow
Number of rows 3622
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 37 0 3622 0
SourceFile 0 1 16 33 0 88 0
RL 0 1 1 3 0 449 0
Max 0 1 2 4 0 80 0
Width 0 1 1 18 0 333 0
DC 0 1 1 18 0 449 0
RRPM 0 1 1 2 0 23 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 10 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 6 6 0 1 0
LineType 0 1 4 6 0 2 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 3622 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 107 0.97 48.1 19.82 2 40 42 49 172 ▁▇▁▁▁
Std 150 0.96 63.7 37.14 2 39 53 83 253 ▇▆▂▁▁
Avg_Length 0 1.00 528.0 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:20:10 1899-12-31 17:47:31 1899-12-31 12:34:39 3364 
Yellow$`RL`= as.integer(Yellow$`RL`)
ggbetweenstats(
  data = Yellow,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

White

White= subset(dat01, Color=="White")
dim(White)
## [1] 3158   25
skim(White)
Data summary
Name White
Number of rows 3158
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 22 40 0 3158 0
SourceFile 0 1 20 36 0 103 0
RL 0 1 1 4 0 585 0
Max 0 1 2 4 0 89 0
Width 0 1 1 18 0 330 0
DC 0 1 1 18 0 604 0
RRPM 0 1 1 2 0 15 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 3 0 7 0
Driven 0 1 1 18 0 71 0
Route 0 1 4 10 0 14 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 7 13 0 28 0
Color 0 1 5 5 0 1 0
LineType 0 1 4 4 0 2 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 3158 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 133 0.96 58.85 38.17 2 41 45 66 282 ▇▃▁▁▁
Std 224 0.93 78.72 117.77 2 42 59 89 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-14 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 13:34:08 3000 
White$`RL`= as.integer(White$`RL`)
ggbetweenstats(
  data = White,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Analysis by LineType

Edge

Edge= subset(dat01, LineType=="Edge")
dim(Edge)
## [1] 4097   25
skim(Edge)
Data summary
Name Edge
Number of rows 4097
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 18 35 0 4097 0
SourceFile 0 1 16 31 0 105 0
RL 0 1 1 3 0 568 0
Max 0 1 2 4 0 96 0
Width 0 1 1 18 0 343 0
DC 0 1 1 18 0 579 0
RRPM 0 1 1 2 0 21 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 5 0 10 0
Driven 0 1 1 18 0 576 0
Route 0 1 3 10 0 15 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 3 0
Type 0 1 7 18 0 5 0
RoutewithDir 0 1 6 13 0 30 0
Color 0 1 5 6 0 2 0
LineType 0 1 4 4 0 1 0
MaterialType 0 1 5 22 0 3 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 4097 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 134 0.97 49.09 22.59 4 40 42 49 172 ▁▇▁▁▁
Std 168 0.96 77.41 99.88 2 41 70 92 1193 ▇▁▁▁▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2024-09-13 00:00:00 11
Time 0 1 1899-12-31 09:17:41 1899-12-31 18:26:01 1899-12-31 12:17:17 3800 
Edge$`RL`= as.integer(Edge$`RL`)
ggbetweenstats(
  data = Edge,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)

Center

Center= subset(dat01, LineType=="Center")
dim(Center)
## [1] 1044   25
skim(Center)
Data summary
Name Center
Number of rows 1044
Number of columns 25
_______________________
Column type frequency:
character 19
logical 1
numeric 3
POSIXct 2
________________________
Group variables None

Variable type: character

skim_variable n_missing complete_rate min max empty n_unique whitespace
Id 0 1 23 37 0 1044 0
SourceFile 0 1 21 33 0 34 0
RL 0 1 1 3 0 293 0
Max 0 1 3 3 0 33 0
Width 0 1 1 18 0 191 0
DC 0 1 1 18 0 243 0
RRPM 0 1 1 2 0 22 0
Driver 0 1 12 12 0 1 0
Length 0 1 2 3 0 5 0
Driven 0 1 1 18 0 91 0
Route 0 1 5 6 0 8 0
Direction 0 1 2 2 0 4 0
Phase 0 1 7 7 0 3 0
Condition 0 1 4 5 0 2 0
Type 0 1 7 18 0 4 0
RoutewithDir 0 1 8 9 0 16 0
Color 0 1 6 6 0 1 0
LineType 0 1 6 6 0 1 0
MaterialType 0 1 5 5 0 1 0

Variable type: logical

skim_variable n_missing complete_rate mean count
ftZERO 1044 0 NaN :

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
Min 34 0.97 46.62 18.01 2 40 40 43 127 ▁▇▁▁▁
Std 75 0.93 45.77 15.21 21 31 45 56 89 ▇▇▆▃▁
Avg_Length 0 1.00 528.00 0.00 528 528 528 528 528 ▁▁▇▁▁

Variable type: POSIXct

skim_variable n_missing complete_rate min max median n_unique
Date 0 1 2023-12-02 00:00:00 2025-04-16 00:00:00 2025-04-13 00:00:00 8
Time 0 1 1899-12-31 10:12:56 1899-12-31 17:03:38 1899-12-31 14:09:52 1013 
Center$`RL`= as.integer(Center$`RL`)
ggbetweenstats(
  data = Center,
  x = Phase,
  y = RL,
  type = "parametric",
  pairwise.comparisons = FALSE,   # hides pairwise comparisons
  pairwise.display = "none",      # just to reinforce
  p.adjust.method = "none",       # disables p-value correction display
  results.subtitle = FALSE,       # hides subtitle with test results
  messages = FALSE                # suppress output messages in console
)+theme_bw(base_size=16)