Nothing
logitpath <- function(x, y, nlam, flmin, ulam, isd, intr, eps, dfmax, pmax, jd,
pf, pf2, maxit, lam2, nobs, nvars, vnames) {
################################################################################
## data setup
y <- as.factor(y)
y <- c(-1, 1)[as.numeric(y)]
if (!all(y %in% c(-1, 1)))
stop("y should be a factor with two levels")
storage.mode(x) <- "double"
################################################################################
## call Fortran core
fit <- .Fortran("loglassoNET", lam2, nobs, nvars, x, as.double(y), jd, pf,
pf2, dfmax, pmax, nlam, flmin, ulam, eps, isd, intr, maxit,
nalam = integer(1), b0 = double(nlam),
beta = double(pmax * nlam), ibeta = integer(pmax),
nbeta = integer(nlam), alam = double(nlam),
npass = integer(1), jerr = integer(1), PACKAGE = "gcdnet")
################################################################################
## output
outlist <- getoutput(fit, maxit, pmax, nvars, vnames)
outlist <- c(outlist, list(npasses = fit$npass, jerr = fit$jerr))
class(outlist) <- c("logitpath")
outlist
}
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.