IIRC what we did in other cases is to convert to ndarray first, and then reshape if the shape is not 2d and corresponds to k_vars, self.exog.shape[0]

Skipper recently did a change like this, but I think we didn't check yet that this is really consistently done across models.