I am trying to run a probit model using
> X1 <- cbind(dat1$daklak, dat1$lamha, dat1$baolam, dat1$num_agesq, dat1$num_edu, dat1$ethnicmaj, dat1$num_exp,