This function is an equivalent of survey:::svy.varcoef
. In the original approach estfun
is calculated from
glm's working residuals:
estfun <- model.matrix(glm.object) * resid(glm.object, "working") * glm.object$weights
In the hopit package, estfun is directly calculated as a gradient (vector of partial derivatives) of the log likelihood function.
Depending on detected design an appropriate survey
function is called.
svy.varcoef_hopit(vcovMat, estfun, design)
a variance-covariance matrix.
a gradient function of the log-likelihood function.
a survey.design
object.
svydesign
hopit