  This paper provides a closed form expression for the pairwise score vector
for the multivariate ordered probit model. This result has several implications
in likelihood-based inference. It is indeed used both to speed-up gradient
based optimization routines for point estimation, and to provide a building
block to compute standard errors and confidence intervals by means of the
Godambe matrix.
