are computed using the following equations.
, where X is the
matrix regression matrix, and V is an
diagonal matrix with ith diagonal term πi(1 – πi). That is, the matrix X is:


is then the jth diagonal term of the matrix
, and the off-diagonal terms are the covariance estimates
for
and
.