Nothing
kfbed.gee.reg <- function(y, x, id, prior = NULL, reps = NULL, univ = NULL, alpha = 0.05, wei = NULL, K = 0:5, test = "testIndGEEReg", correl = "exchangeable", se = "jack") {
a <- MXM::fbed.gee.reg(target = y, dataset = x, id = id, prior = prior, reps = reps, ini = univ, threshold = alpha, wei = wei, K = max(K), test = test, correl = correl, se = se)
info <- a$info
k <- dim(info)[1]
mod <- list()
if ( k > 0 & info[1, 1] > 0 ) {
if (info[k, 1] == info[k - 1, 1]) k <- k - 1
sel <- info[1:k, 1]
for (i in 1:k) mod[[ i ]] <- a$res[1:sel[i], ] ## end if (backward)
names(mod) <- paste("K=", 0:(k-1), sep = "" )
} ## end if (k > 0)
list(res = a, mod = mod)
}
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.