Rtist

 
# ======================================================================== 
# Green Bolts 
# ======================================================================== 
# population size 
size <- 250 
 
# number of bolts 
n <- 50 
 
# number of generations 
gens <- 70 
 
# colors 
greens <- hsv(h = 0.35, s = 1, v = 1, alpha = runif(n, 0.2, 0.8)) 
 
 
png("green_bolts.png", width = 700, height = 400) 
# set graphical parameters 
par(mar = rep(0, 4), bg = 'black') 
# plot 
plot(1:gens, rep(0, gens), type = 'l', ylim = c(0, 1), xaxt = 'n',  
     yaxt = 'n', xlab = "", ylab = "") 
for (j in 1:n) 
{ 
  # initial probability 
  p <- 0.5 
  freq <- c(rep(p, 10), rep(NA, gens - 10)) 
  for (i in 11:gens) 
  { 
    aux <- rbinom(1, size, p) 
    p <- aux / size 
    freq[i] <- p 
  } 
  lwds <- sample(seq(0.25, 3, by = 0.15), 1) 
  # tossing coin 
  coin <- rbinom(1, 1, 0.5) 
  if (coin == 0) 
  { 
    temp <- lowess(freq ~ 1:gens) 
    lines(temp$x, temp$y, col = greens[j], lwd = lwds) 
  } else { 
    lines(1:gens, freq, type = 'l', col = greens[j], lwd = lwds)           
  } 
} 
# signature 
legend("bottomleft", legend = "© Gaston Sanchez", bty = "n",  
       text.col = "gray70") 
# reset par 
par(op) 
dev.off()