Excercise 3

# 1. Running the simulation for 10000 times and increment by 1000.

datalist = list()
for (i in seq(min(1000), max(10000), by = 1000)) {
    
    dat <- data.frame(table(sample(c("Heads", "Tails"), i, rep = T)))
    dat$i <- i 
    datalist[[i]] <- dat 
}
big_data = do.call(rbind, datalist)
#big_data
newdata = subset(big_data, big_data$Var1=='Heads')
newdata$diff = newdata$Freq - (newdata$i/2)
knitr::kable(newdata)
Var1 Freq i diff
1 Heads 517 1000 17
3 Heads 994 2000 -6
5 Heads 1462 3000 -38
7 Heads 1985 4000 -15
9 Heads 2508 5000 8
11 Heads 3014 6000 14
13 Heads 3465 7000 -35
15 Heads 4003 8000 3
17 Heads 4436 9000 -64
19 Heads 5037 10000 37