This function matches the estfun function of the sandwich package and
returns the estimating functions for the fitted model.
It can for instance be used for an OPG estimator of the sigma matrix.
For esreg, the dimension of the estimating functions is n x (kq + ke).