Discussion 8

resultr = 1
resulty = 1
for (i in seq(100)){
  r = .02/(100/i)
  w = (1/6)-2*r
  nw = (1/6)+r
  roll = c(1,2,3,4,5,6)
  pb = c(w,nw,nw,nw,nw,w)
  rsum=0
  for (q in seq(1000000)){
    e = sum(sample(roll,2,replace =TRUE, prob = pb))
    if (e == 7|e == 11){
      rsum = rsum + 1
    }else if(e == 2|e == 3|e == 12){
      rsum = rsum - 1
    }
  }
  resultr[i] = r
  resulty[i] = rsum
  sample = data.frame(cbind(resultr,resulty))
}

ggplot(sample, aes(x = resultr, y= resulty))+geom_point()

hist(resulty)