有没有懂adaptive lasso的大牛,现身说法一下。我用了两个文件包分别做,但是得到的优化系数差别很大呢
一个是lqa
penalty<-adaptive.lasso(lambda=l,as.weights=w)
lm1<- lqa(y=y,x=X,penalty=penalty,weights=w,intercept=F)$coefficients
得到的结果如下:
1.544867e-06 1.740825e-06 1.481752e-06 1.445347e-06 1.263677e-06 1.109394e-06
再一个是用的最近刚出的一个包msgps
ft<-msgps(X,y,penalty="alasso",gamma=1,lambda=l)
coef(ft)
BIC
0.7672238
0.7402016
0.7280660
0.7181718
0.6699783