The GHK algorithm (Geweke, Hajivassiliou and Keane) is an importance sampling method for simulating choice probabilities in the multivariate probit model. These simulated probabilities can be used to recover parameter estimates from the maximized likelihood equation using any one of the usual well known maximization methods (Newton's method, BFGS, etc.). Train has well documented steps for implementing this algorithm for a multinomial probit model. What follows here will apply to the binary multivariate probit model. Consider the case where one is attempting to evaluate the choice probability of Pr ( y i | X i β , Σ ) {\displaystyle \Pr(\mathbf {y_{i}} |\mathbf {X_{i}\beta } ,\Sigma )} where y i = ( y 1 , . . . , y J ) , ( i = 1 , . . . , N ) {\displaystyle \mathbf {y_{i}} =(y_{1},...,y_{J}),\ (i=1,...,N)} and where we can take j {\displaystyle j} as choices and i {\displaystyle i} as individuals or observations, X i β {\displaystyle \mathbf {X_{i}\beta } } is the mean and Σ {\displaystyle \Sigma } is the covariance matrix of the model. The probability of observing choice y i {\displaystyle \mathbf {y_{i}} } is
… excerpt ends here. Continue reading the full article.
