table(replicate(1000000, sum(sample(c(1,0),100,replace=T,prob=c(0.2,0.8)))==25)) FALSE TRUE 955920 44080 table(replicate(1000000, sum(sample(c(1,0),100,replace=T,prob=c(0.2,0.8)))==25)) FALSE TRUE 955886 44114