Calculate speed for modes from distance and duration

Author: usr110

code/scenarios/sao paulo/streamline_travel_survey.R

@@ -47,14 +47,12 @@ rd <- rename(rd, participant_id = ID_PESS ,
 # Recode modes as strings
 mode_df <- data.frame(
   mode_int = append(c(1:17), NA),
-  mode_string = c(rep('bus', 5), 'car_driver', 
-                  'car_passenger', 'taxi',
-                  rep('van', 3), 'subway',
+  mode_string = c(rep('bus', 5), rep('car', 2), 
+                  'taxi', rep('van', 3), 
+                  'subway',
                   'train', 'motorbike',
                   'bicycle', 'walk', 'others', NA)
 
-  
-  
 )
 
 # Convert numeric to string modes
@@ -86,7 +84,7 @@ ggplot(rd %>%
          filter(!is.na(trip_mode)) %>% 
          group_by(trip_mode) %>% 
          summarise(sum_trip_weights = sum(FE_VIA)) %>% 
-         mutate(perc = round(sum_trip_weights/sum_total_trip_weight * 100, 1)), 
+         mutate(perc = round(sum_trip_weights/sum(sum_trip_weights) * 100, 1)), 
        aes(x = trip_mode, y = perc)) + 
   geom_bar(position = 'dodge', stat='identity') +
   geom_text(aes(label = perc), position = position_dodge(width=0.9), vjust=-0.25, color = "blue") +
@@ -118,6 +116,16 @@ ggplot(rd %>%
   theme(axis.text.x = element_text(angle = 90, hjust = 1)) +
   labs(x = "", y = "percentage(%)", title = "Main Mode Distance distribution")
 
+# Calculate mode speed from the dataset, but using mean distance and duration
+# Remove all trips with multiple modes
+# Using only commute mode as a proxy
+mode_speed <- rd %>% filter(is.na(MODO2) & is.na(MODO3) & is.na(MODO4) &  
+              ((MOTIVO_O %in% c(1, 2, 3) & MOTIVO_D %in% 8) | 
+              (MOTIVO_D %in% c(1, 2, 3) & MOTIVO_O %in% 8))
+             ) %>% group_by(trip_mode) %>% 
+  summarise(mean (trip_distance), 
+            speed = (mean(trip_distance)) / (mean(trip_duration) / 60))
+
 
 # source_modes <- c('Bus', 'Walking')
 # target_modes <- c('Private Car')