vcov for multiple linear model

[juniperlsimonis]

isnt calced automatically is a bit tricky, need to think about options a little more maybe this is why the one from multinom was sometimes asymmetric? maybe the asymmetry is important? see what's going on under the hood in the multinom vcov calc for multiple response