3

我使用这个glmnet包来运行多项套索回归。当使用family="multinomial具有p变量和n样本的数据集pmax=x时，如果是奇数，则会发生分段错误x（如果不是pmax>p。在这种情况下，它很可能被忽略，因为它没有影响）。一个例子：

n=100
p=20
require(glmnet)
D= as.data.frame(replicate(p, rnorm(n)))
D[,p] = as.factor(round(rnorm(n)))

lasso  <- glmnet(data.matrix(D[, -p]), D[, p], standardize=T, family="multinomial")         ## works
lasso  <- glmnet(data.matrix(D[, -p]), D[, p], standardize=T, family="multinomial", pmax=7) ## works, because it is odd
lasso  <- glmnet(data.matrix(D[, -p]), D[, p], standardize=T, family="multinomial", pmax=24 ## works, because pmax>p
lasso  <- glmnet(data.matrix(D[, -p]), D[, p], standardize=T, family="multinomial", pmax=10)## crashes

和错误消息：

 *** caught segfault ***
address 0x22de58a8, cause 'memory not mapped'

Traceback:
 1: .Fortran("lognet", parm = alpha, nobs, nvars, nc, as.double(x),     y, offset, jd, vp, cl, ne, nx, nlam, flmin, ulam, thresh,     isd, intr, maxit, kopt, lmu = integer(1), a0 = double(nlam *         nc), ca = double(nx * nlam * nc), ia = integer(nx), nin = integer(nlam),     nulldev = double(1), dev = double(nlam), alm = double(nlam),     nlp = integer(1), jerr = integer(1), PACKAGE = "glmnet")
 2: lognet(x, is.sparse, ix, jx, y, weights, offset, alpha, nobs,     nvars, jd, vp, cl, ne, nx, nlam, flmin, ulam, thresh, isd,     intr, vnames, maxit, kopt, family)
 3: glmnet(data.matrix(D[, -p]), D[, p], standardize = T, family = "multinomial",     pmax = 10)

我的第一个问题是：为什么？这有数学上的原因吗？（我想是这样...）

第二个是：没有比分段错误更好的解决方案吗？......像警告一样？或者只是使用pmax<-pmax-1

编辑：好的，它似乎有点复杂。有时，只有当我第二次为 pmax 执行相同编号的相同命令时，才会发生分段错误。

另外我发现了这个错误

*** glibc detected *** /usr/lib64/R/bin/exec/R: double free or corruption (out): 0x0000000005c41720 ***
======= Backtrace: =========
....

对于 pmax 的偶数和奇数...

现在它对我来说更像是一个错误......或者？

编辑 2：我在 linux 环境（64 位）中使用 glmnet 1.9-5 运行 R 2.15.2 我还在使用 ubuntu 64 位和 R 3.0.2 的另一台 PC 上遇到分段错误

4

2 回答 2

0

这是我在 R 3.0.2、64k、Windows 7、glmnet1.9-5 下得到的

lasso  <- glmnet(data.matrix(D[, -p]), D[, p], standardize=T, family="multinomial", pmax=10)
Warning message:
from glmnet Fortran code (error code -10005); Number of nonzero coefficients along the path exceeds pmax=10 at 5th lambda value; solutions for larger lambdas returned

您没有说明您的设置，因此我无法评论您为什么不捕获错误，但此消息应该提供相当多的洞察力:-) 为何大 pmax 会造成麻烦。

编辑：澄清：我没有遇到段错误。

于 2014-01-10T12:38:58.810 回答
-1

我知道这个答案为时已晚，但以防万一有人发现同样的问题：

我遇到了同样的问题，并按照This Question中的步骤解决了它。基本上，对于您拥有的任何版本的 MATLAB，请尝试

mex -v -setup

它在下降时会给出要使用的编译器列表。安装更新的版本，但不是在 MATLAB 版本之后（对于 2016a，我做了 XE2016 和 VS2015）。然后

mex glmnetMex.F glmnet.f

将完成这项工作。

在 Linux 机器上（使用 2014a），我试过

mex -largeArrayDims glmnetMex.F GLMnet.f

使用gcc编译器，它也可以工作。

希望这可能会有所帮助。我花了好几个星期才弄清楚这一点。

于 2017-04-09T19:46:21.073 回答