Hi, I am wondering if anybody has already written a GHK Simulator function for an ordered probit model in R? We have a system of three equations where the variables can take ordered categorial values. Thanks a million, Pensi