Car Park Simulation

Model simulating movements of cars in a Car Park

The simulation is about the cars entering a public car park area over a period of 12 hours.
We use the Poisson distribution to simulate the number of cars entering during each 30 minutes period.
Next, for each car that enters, we simulate the period it stays in the parking by using the Geometric distribution.

library(tidyverse)

1 - Create the DataFrame

# create the DataFrame structure
# we track the cars entering over each 30 min period, from 8:00 am to 8:00 pm
# date_obs <- as.Date("2023-12-10")
set.seed(25)
n_entering <- rpois(n = 24, lambda = 10)
cars_movements <- data.frame(ent_period = 0, id_car = 0, stay_period = 0)

# populate the DataFrame
for (ii in 1:24) {
    cc <- n_entering[ii]
    ent_cars <- expand_grid(ent_period = ii, id_car = 1:cc)
    ent_cars$id_car <- ent_cars$id_car + 100 * ii
    ent_cars$stay_period = rgeom(n = cc, prob = 0.25)
    cars_movements <- bind_rows(cars_movements, ent_cars)
}

# define the exit period as ent_period + stay_period
cars_movements <- cars_movements %>% 
    mutate(ext_period = ent_period + stay_period)
head(cars_movements, 12)

# define the starting time (the time when the car park opens)
start_time <- as.POSIXct("08:00", format = "%H:%M")
# convert the period number to time (N.B.: there are 1800 seconds in a period of 30 minutes)
cars_movements <- cars_movements %>% 
    mutate(ent_period_t = start_time + ent_period * 1800, ext_period_t = start_time + ext_period * 1800) %>% 
    filter(id_car != 0)

head(cars_movements, 12)

# range of entering and exiting periods
range(cars_movements$ent_period_t)

[1] "2024-04-07 08:30:00 CEST" "2024-04-07 20:00:00 CEST"

range(cars_movements$ext_period_t)

[1] "2024-04-07 09:00:00 CEST" "2024-04-08 04:30:00 CEST"

2 - Line Plot of cars entering and exiting over time

# set plotting area width in inches 
#options(repr.plot.width = 12)

# define the graphic object ent_plot
ent_plot <- cars_movements %>% group_by(ent_period_t) %>% 
  summarize(num_cars = n()) %>% 
  ggplot(aes(x = ent_period_t, y = num_cars)) + geom_line(linewidth = 2, color = "forestgreen")
# define the ext_data
ext_data <- cars_movements %>% group_by(ext_period_t) %>% 
  summarize(num_cars = n()) 
# plot the cars entering and exiting 
ent_plot + 
  geom_line(data = ext_data, aes(x = ext_period_t, y = num_cars), linewidth = 2, color = "red") + 
    scale_x_datetime(limits = c(min(cars_movements$ent_period_t), max(cars_movements$ent_period_t)), 
                     date_labels = "%H:%M", date_breaks = "1 hour")

# plot the cars entering and exiting (the latter with negative sign) 
ent_plot + 
  geom_line(data = ext_data, aes(x = ext_period_t, y = -num_cars), linewidth = 2, color = "red") + 
    scale_x_datetime(limits = c(min(cars_movements$ent_period_t), max(cars_movements$ent_period_t)), 
                     date_labels = "%H:%M", date_breaks = "1 hour") + geom_hline(aes(yintercept = 0))