One page 26 of the Mata manual [M-5] for the function optimize() (and in the corresponding Stata help file) it says
Instead it should be
optimize_result_V_robust(S) returns H*invsym(S'S)*H, which is the robust estimate of variance
optimize_result_V_robust(S) returns invsym(-H)*(S'S)*invsym(-H), which is the robust estimate of variance
Comment