Given A = {a_1,…,a_n}, the response y, and a set of edges E, this function aims to solves
min 1/2Axy^2 + λ_1x_1 + λ_2 ∑_{(i,j) in E}w_(i,j)x_ir_(i,j)x_j
where w_(i,j) is the weight of the edge (i,j), and r_(i,j) is the sign of the correlation between features a_i and a_j. The weight and sign can be specified in Rwt: w=Rwt, and r=sign(Rwt).
1 2 3 4 
A 

y 

tp 

s1 

s2 

RmaxIter 

RaMaxIter 

Rrho 

Rtau 

Rwt 

Rtol 

RaTol 

x0 

Returned value x0 is the solution to the optimizaiton problem.
Yiwen Sun
S.Kim and E.Xing. Statisctical estimation of correlated genome associations to a quantitative trait network. PLoS genetics, 5(8):e1000587, 2009
1 2 3 4 
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.