I am running a heckprobit for selection and clustering the standard errors at the state level. I am wondering how the standard errors are calculated for the main effects since the square root of the diagonal of the variance covariance matrix does not give the same values that are reported in the output of the heckprobit. However, for the selection equation the standard error reported is the square root of the variance covariance matrix.