For you Stata users, I'm wondering what's the best SAS procedure to use to replicate the output from Stata xtprobit when using two-way panel data (repeated cross section/time series). I've tried PROC LOGISTIC with a probit link function and two class variables. Also PROC GLIMMIX with double random intercept effects. I've also tried PROC NLMIXED, but I haven't been able to figure out if it's possible to use more than one random effect. None of the results I've gotten from SAS are very compatible with the output from Stata. Also the SAS functions run for a long time compared with the procedure in Stata. Any help would be much appreciated!