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