Nothing
glm_profiling <-
function(X, Y, W, etastart, O, fam, glmcontrol, est, OriginalDeviance, DispersionParameter, K, b){
z <- numeric(length(b))
pmi <- matrix(nrow=length(b), ncol=length(est))
obj <- list()
Kest <- as.vector(K %*% est)
signs <- rep(1, length(b))
signs[Kest > b] <- -1
bdist <- abs(b-Kest)
nbi <- order(bdist[signs < 0])
pbi <- order(bdist[signs > 0]) + max(nbi)
bind <- c(nbi,pbi)
oeta <- logical(length=length(b))
oeta[max(nbi + 1)] <- TRUE
etavec <- etastart
for (i in bind){
fm <- try(orglm.fit(x=X, y=Y, weights=W, etastart=etavec, offset=O, family=fam, control=glmcontrol, constr=K, rhs=b[i], nec=1), silent=TRUE)
if (class(fm)[1] != "try-error"){
zz <- signs[i]*sqrt(abs(fm$deviance - OriginalDeviance)/DispersionParameter)
z[i] <- zz
pmi[i,] <- fm$coefficients
obj[[i]] <- fm
if (fm$converged == TRUE & oeta[i] == FALSE) etavec <- X %*% fm$coefficients
if (oeta[i] == TRUE) etavec <- etastart
} else {
z[i] <- NA
obj[[i]] <- NA
}
}
list(stats=cbind(b, z), param=pmi, object=obj)
}
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.