Obtain the heuristic residuals.
hresid(y, svar, design.mat, o, initial.res.par = c(0, 0), ...)
vector of numeric. The outcomes.
vector of numeric. The strata variable.
matrix. The design matrix including intercept term.
matrix. Contains the estimated coefficients of the RO-logit and their SEs.
The initial values of the intercept and log(scale), to
be passed to the optim
function. The default values are set to
c(0, 0)
, yet users are recommended to try a few initial values to
make sure global optimum is reached.
Other parameters to be passed to the optim
function. See
optim
for details.
Returns a list containing the estimated intercept and
log(scale)
, the covariance matrix of these two parameters,
convergence status from optim
, and the heuristic residuals.