Nothing
picasso <- function(X,
Y,
lambda = NULL,
nlambda = 100,
lambda.min.ratio = 0.05,
family = "gaussian",
method = "l1",
type.gaussian = "naive",
gamma = 3,
df = NULL,
standardize = TRUE,
intercept = TRUE,
prec = 1e-7,
max.ite = 1e3,
verbose = FALSE)
{
if (family != "gaussian" && family != "binomial" && family != "poisson" && family != "sqrtlasso"){
cat(" Wrong \"family\" input. \n \"family\" should be
one of \"gaussian\", \"binomial\" and \"poisson\".\n",
family," is not supported in this version. \n")
return(NULL)
}
if (family == "gaussian") {
if (!is.matrix(Y))
Y = as.matrix(Y)
p = ncol(Y)
if (p == 1)
out = picasso.gaussian(X = X, Y = Y, lambda = lambda, nlambda = nlambda,
lambda.min.ratio = lambda.min.ratio,
method = method, type.gaussian = type.gaussian, gamma = gamma, df = df,
standardize = standardize, intercept= intercept,
prec = prec,
max.ite = max.ite, verbose = verbose)
}
if (family == "binomial") {
if(!is.matrix(Y))
Y = as.matrix(Y)
out = picasso.logit(X = X, Y = Y, lambda = lambda, nlambda = nlambda,
lambda.min.ratio = lambda.min.ratio,
method = method, gamma = gamma, standardize = standardize, intercept=intercept,
prec = prec, max.ite = max.ite, verbose = verbose)
}
if (family == "sqrtlasso"){
if(!is.matrix(Y))
Y = as.matrix(Y)
out = picasso.sqrtlasso(X = X, Y = Y, lambda = lambda, nlambda = nlambda,
lambda.min.ratio = lambda.min.ratio,
method = method, gamma = gamma, standardize = standardize, intercept=intercept,
prec = prec, max.ite = max.ite, verbose = verbose)
}
if(family=="poisson"){
out = picasso.poisson(X = X, Y=Y, lambda = lambda, nlambda = nlambda,
lambda.min.ratio = lambda.min.ratio,
method = method, gamma = gamma,
standardize = standardize,
intercept = intercept,
prec = prec, max.ite = max.ite,
verbose = verbose)
}
out$family = family
return(out)
}
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.