This function fits a multi-response pliable lasso model over a path of regularization values.
admm_MADMMplasso_cpp(
beta0,
theta0,
beta,
beta_hat,
theta,
rho1,
X,
Z,
max_it,
W_hat,
XtY,
y,
N,
e_abs,
e_rel,
alpha,
lambda,
alph,
svd_w_tu,
svd_w_tv,
svd_w_d,
C,
CW,
gg,
my_print = TRUE
)predicted values for the ADMM part
a vector of length ncol(y) of estimated beta_0 coefficients
matrix of the initial theta_0 coefficients ncol(Z) by ncol(y)
a matrix of the initial beta coefficients ncol(X) by ncol(y)
a matrix of the initial beta and theta coefficients (ncol(X)+ncol(X) by ncol(Z)) by ncol(y)
an array of initial theta coefficients ncol(X) by ncol(Z) by ncol(y)
the Lagrange variable for the ADMM which is usually included as rho in the MADMMplasso call.
n by p matrix of predictors
n by nz matrix of modifying variables. The elements of z may represent quantitative or categorical variables, or a mixture of the two. Categorical variables should be coded by 0-1 dummy variables: for a k-level variable, one can use either k or k-1 dummy variables.
maximum number of iterations in loop for one lambda during the ADMM optimization. This is usually included in the MADMMplasso call
N by (p+(p by nz)) of the main and interaction predictors. This generated internally when MADMMplasso is called or by using the function generate_my_w.
a matrix formed by multiplying the transpose of X by y.
N by D matrix of responses. The X and Z variables are centered in the function. We recommend that X and Z also be standardized before the call
nrow(X)
absolute error for the ADMM. This is included int the call of MADMMplasso.
relative error for the ADMM. This is included int the call of MADMMplasso.
mixing parameter, usually obtained from the MADMMplasso call. When the goal is to include more interactions, alpha should be very small and vice versa.
a vector lambda_3 values for the ADMM call with length ncol(y). This is usually calculated in the MADMMplasso call. In our current setting, we use the same the lambda_3 value for all responses.
an overrelaxation parameter in [1, 1.8], usually obtained from the MADMMplasso call.
the transpose of the U matrix from the SVD of W_hat
the transpose of the V matrix from the SVD of W_hat
the D matrix from the SVD of W_hat
the trained tree
weights for the trained tree The easy way to obtain this is by using the function (tree_parms) which gives a default clustering. However, user decide on a specific structure and then input a tree that follows such structure.
penalty terms for the tree structure for lambda_1 and lambda_2 for the ADMM call.
Should information form each ADMM iteration be printed along the way? Default TRUE. This prints the dual and primal residuals