Lasso and Glmnet library

15 次查看（过去 30 天）
jason beckell
jason beckell 2018-4-22
评论： Le Zheng 2020-12-30
Dear all,
I have a question about Lasso with a high-dimensional problem, both with the standard Lasso procedure in matlab and with the wrapper of the Glmnet Library. As to the standard procedure in Matlab (which solves the problem by means of ADMM algorithm), I can always get a solution which is sometimes less accurate than that produced by the wrapper of Glmnet library in Matlab. But, when I use the wrapper, I sometimes get this exception:
Error using ==
Matrix dimensions must agree.
Error in cvglmnet (line 358)
idmin = options.lambda==CVerr.lambda_min;
Even if I explicitly declare the sequence of lambdas to be used, I still might get this exception. I have also seen that, in this case, the structure returned from the following lines code of function elnet (lines 23-24):
task = 11;
[lmu,a0,ca,ia,nin,rsq,alm,nlp,jerr] = glmnetMex(task,parm,x,y-offset,jd,vp,ne,nx,nlam,flmin,ulam,thresh,isd,weights,ka,cl,intr,maxit);
alm does not contain the sequence of lambdas required, but a vector of NaNs. What this might be due to? And what can I do to solve this problem?
Many thank for your attention and my best regards,
Jason.
  1 个评论
Le Zheng
Le Zheng 2020-12-30
Hi Jason, I'm facing the same problem, too. Did you find a solution?

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Albena Vassileva
Albena Vassileva 2019-7-9
编辑：Albena Vassileva 2019-7-30
EDIT: Figured it out. Or rather, "patched it".. The input should be in double precision. Still not sure exactly why it happened...
But I had it in single and was getting this error too. I can only assume it went wrong somewhere around the initialization of the lambda search, and that's why it was sometimes working.
  1 个评论
Le Zheng
Le Zheng 2020-12-29
Hi Albena, I'm facing the same problem, too. But the error still comes out when I use cmd double to convert my data to double precision. Do you have any clues?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by