set.seed(123)
# simulasi Poisson
kedatangan <- rpois(100, lambda = 5)
# ringkasan
mean(kedatangan)
## [1] 5.01
table(kedatangan)
## kedatangan
## 0 1 2 3 4 5 6 7 8 9 10 11
## 1 1 9 15 19 16 15 10 8 4 1 1
# grafik
barplot(table(kedatangan),
main = "Distribusi Kedatangan Pelanggan",
xlab = "Jumlah pelanggan per jam",
ylab = "Frekuensi",
col = "lightgreen")
# Jika rata-rata simulasi mendekati 5, maka model Poisson sesuai untuk
menggambarkan pola kedatangan.
set.seed(123)
# simulasi normal
waktu <- rnorm(200, mean = 10, sd = 2)
# ringkasan
mean(waktu)
## [1] 9.982859
sd(waktu)
## [1] 1.88632
# histogram
hist(waktu,
main = "Distribusi Waktu Pelayanan",
xlab = "Menit",
col = "skyblue")
# proporsi > 12 menit
mean(waktu > 12)
## [1] 0.14
#STUDI KASUS SIMULASI VARIABEL RANDOM # Judul: Simulasi Waktu Pelayanan dan Jumlah Kedatangan Pelanggan pada Sistem Antrian Minimarket
# =================================
# SIMULASI SISTEM MINIMARKET
# =================================
set.seed(123)
n <- 1000
# -----------------------------
# 1. Simulasi jumlah pelanggan
# -----------------------------
lambda <- 8
pelanggan <- rpois(n, lambda)
# -----------------------------
# 2. Simulasi waktu pelayanan
# -----------------------------
mu <- 5
sigma <- 1
pelayanan <- rnorm(n, mean = mu, sd = sigma)
# pastikan tidak negatif
pelayanan[pelayanan < 0] <- 0.1
# -----------------------------
# Statistik deskriptif
# -----------------------------
mean(pelanggan)
## [1] 7.968
var(pelanggan)
## [1] 7.834811
mean(pelayanan)
## [1] 5.011933
sd(pelayanan)
## [1] 1.001516
# -----------------------------
# Probabilitas penting
# -----------------------------
# P(jumlah pelanggan > 12)
prob_pelanggan <- mean(pelanggan > 12)
# P(waktu pelayanan > 7 menit)
prob_pelayanan <- mean(pelayanan > 7)
prob_pelanggan
## [1] 0.063
prob_pelayanan
## [1] 0.026
# -----------------------------
# Visualisasi
# -----------------------------
# Poisson
barplot(table(pelanggan),
main = "Distribusi Jumlah Pelanggan",
xlab = "Jumlah pelanggan",
ylab = "Frekuensi",
col = "lightgreen")
# Normal
hist(pelayanan,
probability = TRUE,
main = "Distribusi Waktu Pelayanan",
xlab = "Menit",
col = "lightblue")
lines(density(pelayanan), col = "red", lwd = 2)
# Interpretasi Hasil # Hasil simulasi menunjukkan bahwa: