#====================
# Binary event simulator function
# generates n 0's or 1's, with p
# being the probability of a 1.
#====================
outcomes=function(n,p) {
  u=runif(n)
  x=numeric(n)
  for (i in 1:n) {
    if (u[i]<=p) x[i]=1 else x[i]=0
  }
  return(x)
}