curve(dnorm(x, mean = 0, sd = 1), from = -4, to = 4,
main = "Kurva Distribusi Normal",
xlab = "X", ylab = "Kepadatan (Density)",
col = "blue", lwd = 2)

# Parameter untuk 3 distribusi Binomial yang berbeda
n1 = 10; p1 = 0.5 # mu = n*p, varians = n*p*(1-p)
n2 = 20; p2 = 0.3
n3 = 30; p3 = 0.2
x1 = 0:n1
x2 = 0:n2
x3 = 0:n3
y1 = dbinom(x1, n1, p1) # dbinom = fungsi pmf binomial
y2 = dbinom(x2, n2, p2)
y3 = dbinom(x3, n3, p3)
mu1 = n1*p1; var1 = n1*p1*(1-p1)
mu2 = n2*p2; var2 = n2*p2*(1-p2)
mu3 = n3*p3; var3 = n3*p3*(1-p3)
# Plot distribusi pertama
plot(x1, y1, type = "h", lwd = 2, col = "blue",
xlim = c(0, 30), ylim = c(0, 0.35),
main = "Distribusi Binomial dengan Mu dan Varians Berbeda",
xlab = "X", ylab = "P(X = x)")
points(x1, y1, pch = 16, col = "blue")
# Tambahkan distribusi kedua dan ketiga di plot yang sama
lines(x2, y2, type = "h", lwd = 2, col = "red")
points(x2, y2, pch = 16, col = "red")
lines(x3, y3, type = "h", lwd = 2, col = "darkgreen")
points(x3, y3, pch = 16, col = "darkgreen")
legend("topright",
legend = c(paste0("n=", n1, ", p=", p1, " (μ=", mu1, ", σ²=", round(var1,2), ")"),
paste0("n=", n2, ", p=", p2, " (μ=", mu2, ", σ²=", round(var2,2), ")"),
paste0("n=", n3, ", p=", p3, " (μ=", mu3, ", σ²=", round(var3,2), ")")),
col = c("blue", "red", "darkgreen"), lwd = 2)

# Parameter untuk 3 distribusi normal yang berbeda
mu1 = 0; sigma1 = 1 # varians = sigma1^2 = 1
mu2 = 0; sigma2 = 2 # varians = sigma2^2 = 4
mu3 = 3; sigma3 = 1.5 # varians = sigma3^2 = 2,25
x = seq(-10, 10, length = 500)
y1 = dnorm(x, mean = mu1, sd = sigma1)
y2 = dnorm(x, mean = mu2, sd = sigma2)
y3 = dnorm(x, mean = mu3, sd = sigma3)
# Plot kurva pertama
plot(x, y1, type = "l", lwd = 2, col = "blue",
ylim = c(0, 0.45),
main = "Distribusi Normal dengan Mu dan Varians Berbeda",
xlab = "X", ylab = "Kepadatan (Density)")
# Tambahkan kurva kedua dan ketiga di plot yang sama
lines(x, y2, lwd = 2, col = "red")
lines(x, y3, lwd = 2, col = "darkgreen")
legend("topright",
legend = c(paste0("μ=", mu1, ", σ²=", sigma1^2),
paste0("μ=", mu2, ", σ²=", sigma2^2),
paste0("μ=", mu3, ", σ²=", sigma3^2)),
col = c("blue", "red", "darkgreen"), lwd = 2)

# Data sesuai tabel
kurva = c(1, 2, 3, 4, 5)
rata_rata = c(0, 2, 5, 10, 15)
simpangan_baku = c(1, 1, 5, 3, 0.1)
# Membagi jendela grafik menjadi 2 baris x 3 kolom (5 grafik + 1 kosong)
par(mfrow = c(2, 3))
for (i in 1:5) {
mu = rata_rata[i]
sigma = simpangan_baku[i]
# Rentang x disesuaikan per kurva (mean +/- 4*sd) agar bentuk kurva terlihat jelas
x = seq(mu - 4*sigma, mu + 4*sigma, length = 200)
y = dnorm(x, mean = mu, sd = sigma)
plot(x, y, type = "l", lwd = 2, col = "blue",
main = paste0("Kurva ", kurva[i], " (μ=", mu, ", σ=", sigma, ")"),
xlab = "X", ylab = "Kepadatan (Density)")
}

mu <- c(0, 2, 5, 10, 1.5) # Nilai Rata-rata (mu)
sigma <- c(1, 1, 5, 3, 0.1) # Nilai Simpangan Baku (sigma)
x <- seq(-5, 20, length.out = 500)
plot(x, dnorm(x, mean = mu[1], sd = sigma[1]), type = "n",
main = "Simulasi Kurva Distribusi Normal Tertumpuk",
xlab = "X", ylab = "Densitas", ylim = c(0, 4), xlim = c(-4, 18))
warna <- c("blue", "red", "green", "purple", "darkorange")
for (i in 1:length(mu)) {
y <- dnorm(x, mean = mu[i], sd = sigma[i])
lines(x, y, col = warna[i], lwd = 2)
polygon(c(x[1], x, x[length(x)]), c(0, y, 0),
col = adjustcolor(warna[i], alpha.f = 0.2),
border = warna[i])
}
legend("topright",
legend = paste("Kurva", 1:5, "(mu =", mu, ", sigma =", sigma, ")"),
fill = adjustcolor(warna, alpha.f = 0.2),
border = warna,
cex = 0.7)
