The data I have generated is as follows:
set.seed(100) n = 100 c1_prob = 0.8 X = matrix(0, nrow = n, ncol = 2) y = matrix(0,nrow=n,ncol=1) for (i in 1:n){ if(ru