Nothing
`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
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.