R/separation.detection.R

`separation.detection` <-
function (fit, nsteps = 30) 
{
    fit.class <- class(fit)[1]
    if (fit.class != "glm") 
        stop("Only objects of class 'glm' are accepted.")
    eps <- .Machine$double.eps
    betasNames <- names(betas <- coef(fit))
    noNA <- !is.na(betas)
    stdErrors <- matrix(0, nsteps, length(betas))
    for (i in 1:nsteps) {
        suppressWarnings(temp.fit <- update(fit, control = glm.control(maxit = i, 
            epsilon = eps)))
        stdErrors[i, noNA] <- summary(temp.fit)$coef[betasNames[noNA], "Std. Error"]
    }
    res <- sweep(stdErrors, 2, stdErrors[1, ], "/")
    colnames(res) <- names(coef(fit))
    res
}

Try the brglm package in your browser

Any scripts or data that you put into this service are public.

brglm documentation built on April 6, 2019, 3 a.m.