Simulasi Variabel Random

Simulasi Distribusi Diskrit: Distribusi Binomial

n <- 30  # Jumlah percobaan
p <- 0.5  # Probabilitas sukses
binomial_dist <- rbinom(n, size = n, prob = p)

hist(binomial_dist, breaks = 30, main = "Histogram Distribusi Binomial", xlab = "Jumlah Sukses", col = "lightpink")

Simulasi Distribusi Kontinu: Distribusi Eksponensial

rate <- 1       # Parameter rate
n <- 1000
exp_dist <- rexp(n, rate)
hist(exp_dist, breaks = 30, main = "Histogram Distribusi Eksponensial", xlab = "Nilai", col = "lightcoral")

Latihan Studi Kasus

Studi Kasus: Simulasi Jumlah Kecelakaan Mobil

set.seed(48)
lambda <- 3     # Rata-rata kecelakaan per minggu
x <- 24          # Banyak pengamatan
crash_dist <- rpois(x, lambda)
crash_dist

##  [1] 3 0 0 2 2 1 7 0 4 3 4 2 5 3 3 5 4 4 2 3 2 2 8 3

mean_crash <- mean(crash_dist)
cat("Rata-rata jumlah kecelakaan per minggu:", mean_crash, "\n")

## Rata-rata jumlah kecelakaan per minggu: 3

prob_above_5 <- sum(crash_dist > 5) / x
cat("Probabilitas jumlah kecelakaan lebih dari 5:", prob_above_5, "\n")

## Probabilitas jumlah kecelakaan lebih dari 5: 0.08333333