def _HCCM(results, scale):
 H = np.dot(results.model.pinv_wexog, (scale[:, None] * results.model.pinv_wexog.T))
 return H

def _HCCM ( results , scale ) : H = np . dot ( results . model . pinv_wexog , ( scale [ : , None ] * results . model . pinv_wexog . T ) ) return H
sandwich with pinv(x) * diag * pinv(x) .