Note
that for the robust WLS
variants, we use the diagonal of the weight matrix
for estimation, but we use the full weight matrix to correct the standard
errors and to compute the test statistic.
what can I do when data is endogenous ordinal data and is not multivariate normal?