tiffany's question re gambler's fallacy

die_rolls <- 1:6
storage <- c()

for(i in 1:1000)
  {
  rolls_for_the_week <- sample(die_rolls, 7, replace = TRUE)
  storage[i] <- sum(rolls_for_the_week == 1) == 0 # if there's no 1, then TRUE
}

cat("\nthe probability of not rolling a 1 all week is estimated =", 
    sum(storage)/1000, "\n")

cat("\nthe probability of rolling at least a single 1 all week is estimated =",
    1 - sum(storage)/1000, "\n")