# ========================================================================
# Blue Lightnings
# ========================================================================
# size
n <- 500
# generations
gens <- 50
# colors
blues <- hsv(h = 0.65, s = runif(20, 0.8, 1), v = 1,
alpha = runif(20, 0.6, 0.9))
whites <- hsv(h = 0, s = 0, v = 1, alpha = runif(20, 0.1, 0.7))
png("blue_lightnings.png", width = 700, height = 400)
# set graphical parameters
op <- par(bty = "n", mar = c(0, 0.2, 0, 0.2), bg = 'black')
# open plot
plot(1:gens, rep(0, 50), type = 'l', xlim = c(1, gens), ylim = c(0.15, 0.85),
axes = FALSE)
#rect(0, 0, 50, 1, col = "black")
for (j in 1:25)
{
# initial probability
p <- 0.5
# empty vector of frequencies
freq <- c(p, rep(NA, gens-1))
# get generations
for (i in 2:gens)
{
aux <- rbinom(1, n, p)
p <- aux / n
freq[i] <- p
}
lwds <- sample(seq(0.5, 2.5, by = 0.15), 1)
lines(1:gens, freq, type = 'l', col = whites[j], lwd = lwds)
lines(1:gens, freq, type = 'l', col = blues[j], lwd = lwds)
}
# signature
legend("bottomleft", legend = "© Gaston Sanchez", bty = "n",
text.col = "gray70")
# reset par
par(op)
dev.off()