Meat of U sandwich estimator
sandwich_meat(y, X, b, link, w)
numeric A vector of pseudo observations.
matrix The pairwise difference of original design matrix.
numeric The PIM parameter.
character The link function: "logit", or "probit".
numeric The weights, default is NULL.
matrix A p x p
matrix.