goscar(A, y, tp, s1, s2, RmaxIter = 100,
RaMaxIter = 1000, Rrho = 5, Rtau = 0.15,
Rwt = rep(1, length(tp)), Rtol = 0.001,
RaTol = 0.001, x0 = rep(0, ncol(A)))A The data matrix of size $n \times p$, each row corresponds to one sample.y The response vector of length n.tp The edges vector of length 2*g (eg. (1,2,3,4) means an edge between 1 and 2, and an edge between 3 and 4, g=2 is the number of edges).s1 The $l_1$ regularization parameter, $s1 >=0$.s2 Tge grouping penatly parameter, $s2 >=0$.RmaxIter The maximum number of iterations in DC programming (default 100).RaMaxIter The maximum number of iterations in ADMM (default 1000).Rrho The dual update length ofor ADMM (default 5).Rtau The tuning parameter for non-convex penalty (default 0.15).Rwt The weight and signs of edges (default rep(1,g)).Rtol The tolerance for convergence in DC programming (default 1e-3).RaTol The tolerance for convergence in ADMM (default 1e-3).x0 The returned weight vector (default rep(0,p)).A<-matrix(rnorm(25),5,5)
y<-rnorm(5)
tp<-c(1,2,2,3,3,4,4,5)
goscar(A,y,tp,0,0)Run the code above in your browser using DataLab