[Lme4-authors] Crash in the latest release
Thanks, posted to https://github.com/lme4/lme4/issues/412 for further discussion ... On Wed, Feb 22, 2017 at 10:03 AM, Therneau, Terry M., Ph.D.
<therneau at mayo.edu> wrote:
I found this by accident yesterday. The program that crashes is the first two lines of the example from the help page for nlmer. That example hasn't changed in a long time, so I assumed that it is an R-devel issue. It could also be a long latent nlmer bug. The second run with valgrind is puzzling. Terry T.
library(lmer) sessionInfo()
R Under development (unstable) (2017-02-21 r72241) Platform: i686-pc-linux-gnu (32-bit) Running under: Ubuntu 14.04.1 LTS Matrix products: default BLAS: /usr/local/src/R-devel/lib/libRblas.so LAPACK: /usr/local/src/R-devel/lib/libRlapack.so locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=C [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] lme4_1.1-12 Matrix_1.2-8 loaded via a namespace (and not attached): [1] minqa_1.2.4 MASS_7.3-46 compiler_3.4.0 tools_3.4.0 [5] Rcpp_0.12.9 splines_3.4.0 nlme_3.1-131 grid_3.4.0 [9] nloptr_1.0.4 lattice_0.20-34
startvec <- c(Asym = 200, xmid = 725, scal = 350) (nm1 <- nlmer(circumference ~ SSlogis(age, Asym, xmid, scal) ~ Asym|Tree,
+ Orange, start = startvec)) *** Error in `/usr/local/src/R-devel/bin/exec/R': malloc(): memory corruption: 0x0e3df988 *** Aborted (core dumped) --------------------- The valgrind version is interesting too: it complains and then finishes! tmt-local2437% R -d valgrind ==31901== Memcheck, a memory error detector ==31901== Copyright (C) 2002-2013, and GNU GPL'd, by Julian Seward et al. ==31901== Using Valgrind-3.10.0.SVN and LibVEX; rerun with -h for copyright info ==31901== Command: /usr/local/src/R-devel/bin/exec/R ==31901== R Under development (unstable) (2017-02-21 r72241) -- "Unsuffered Consequences" Copyright (C) 2017 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. [Previously saved workspace restored]
library(lme4)
Loading required package: Matrix
startvec <- c(Asym = 200, xmid = 725, scal = 350) (nm1 <- nlmer(circumference ~ SSlogis(age, Asym, xmid, scal) ~ Asym|Tree,
+ Orange, start = startvec)) ==31901== Invalid write of size 4 ==31901== at 0xA13DA20: cholmod_start (cholmod_common.c:184) ==31901== by 0xD8B3065: cholmod_start (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0xD8B728D: Eigen::CholmodBase<Eigen::SparseMatrix<double, 0, int>, 1, Eigen::CholmodDecomposition<Eigen::SparseMatrix<double, 0, int>, 1>
::CholmodBase() (in /home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0xD8B413D: Eigen::CholmodDecomposition<Eigen::SparseMatrix<double, 0, int>, 1>::CholmodDecomposition() (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0xD8B39B8: lme4::lme4CholmodDecomposition<Eigen::SparseMatrix<double, 0, int>, 1>::lme4CholmodDecomposition() (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0xD8AEBE4: lme4::merPredD::merPredD(SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*) (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0xD8727F1: merPredDCreate (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0x80C0B9E: R_doDotCall (dotcode.c:648) ==31901== by 0x80C159F: do_dotcall (dotcode.c:1252) ==31901== by 0x80FFABA: Rf_eval (eval.c:728) ==31901== by 0x810308A: do_set (eval.c:2579) ==31901== by 0x80FF907: Rf_eval (eval.c:700) ==31901== Address 0xe310828 is 0 bytes after a block of size 2,408 alloc'd ==31901== at 0x402A6DC: operator new(unsigned int) (in /usr/lib/valgrind/vgpreload_memcheck-x86-linux.so) ==31901== by 0xD872770: merPredDCreate (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0x80C0B9E: R_doDotCall (dotcode.c:648) ==31901== by 0x80C159F: do_dotcall (dotcode.c:1252) ==31901== by 0x80FFABA: Rf_eval (eval.c:728) ==31901== by 0x810308A: do_set (eval.c:2579) ==31901== by 0x80FF907: Rf_eval (eval.c:700) ==31901== by 0x8101F20: do_begin (eval.c:2187) ==31901== by 0x80FF907: Rf_eval (eval.c:700) ==31901== by 0x810141B: R_execClosure (eval.c:1614) ==31901== by 0x80FF711: Rf_eval (eval.c:747) ==31901== by 0x8101F20: do_begin (eval.c:2187) ==31901== ==31901== Invalid write of size 4 ==31901== at 0xA13DA31: cholmod_start (cholmod_common.c:184) ==31901== by 0xD8B3065: cholmod_start (in /home/therneau/Rlib/lme4/libs/lme4.so) ==31901== by 0xD8B728D: Eigen::CholmodBase<Eigen::SparseMatrix<double, 0, int>, 1, Eigen::CholmodDecomposition<Eigen::SparseMatrix<double, 0, int>, 1>
::CholmodBase() (in /home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0xD8B413D:
Eigen::CholmodDecomposition<Eigen::SparseMatrix<double, 0, int>,
1>::CholmodDecomposition() (in /home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0xD8B39B8:
lme4::lme4CholmodDecomposition<Eigen::SparseMatrix<double, 0, int>,
1>::lme4CholmodDecomposition() (in /home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0xD8AEBE4: lme4::merPredD::merPredD(SEXPREC*, SEXPREC*,
SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*,
SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*, SEXPREC*,
SEXPREC*) (in /home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0xD8727F1: merPredDCreate (in
/home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0x80C0B9E: R_doDotCall (dotcode.c:648)
==31901== by 0x80C159F: do_dotcall (dotcode.c:1252)
==31901== by 0x80FFABA: Rf_eval (eval.c:728)
==31901== by 0x810308A: do_set (eval.c:2579)
==31901== by 0x80FF907: Rf_eval (eval.c:700)
==31901== Address 0xe31082c is 4 bytes after a block of size 2,408 alloc'd
==31901== at 0x402A6DC: operator new(unsigned int) (in
/usr/lib/valgrind/vgpreload_memcheck-x86-linux.so)
==31901== by 0xD872770: merPredDCreate (in
/home/therneau/Rlib/lme4/libs/lme4.so)
==31901== by 0x80C0B9E: R_doDotCall (dotcode.c:648)
==31901== by 0x80C159F: do_dotcall (dotcode.c:1252)
==31901== by 0x80FFABA: Rf_eval (eval.c:728)
==31901== by 0x810308A: do_set (eval.c:2579)
==31901== by 0x80FF907: Rf_eval (eval.c:700)
==31901== by 0x8101F20: do_begin (eval.c:2187)
==31901== by 0x80FF907: Rf_eval (eval.c:700)
==31901== by 0x810141B: R_execClosure (eval.c:1614)
==31901== by 0x80FF711: Rf_eval (eval.c:747)
==31901== by 0x8101F20: do_begin (eval.c:2187)
==31901==
Nonlinear mixed model fit by maximum likelihood ['nlmerMod']
Formula: circumference ~ SSlogis(age, Asym, xmid, scal) ~ Asym | Tree
Data: Orange
AIC BIC logLik deviance df.resid
273.1438 280.9205 -131.5719 263.1438 30
Random effects:
Groups Name Std.Dev.
Tree Asym 31.646
Residual 7.843
Number of obs: 35, groups: Tree, 5
Fixed Effects:
Asym xmid scal
192.1 727.9 348.1
q()
==31901== ==31901== HEAP SUMMARY: ==31901== in use at exit: 94,537,359 bytes in 46,733 blocks ==31901== total heap usage: 194,149 allocs, 147,416 frees, 248,586,365 bytes allocated ==31901== ==31901== LEAK SUMMARY: ==31901== definitely lost: 0 bytes in 0 blocks ==31901== indirectly lost: 0 bytes in 0 blocks ==31901== possibly lost: 0 bytes in 0 blocks ==31901== still reachable: 94,537,359 bytes in 46,733 blocks ==31901== suppressed: 0 bytes in 0 blocks ==31901== Rerun with --leak-check=full to see details of leaked memory ==31901== ==31901== For counts of detected and suppressed errors, rerun with: -v ==31901== ERROR SUMMARY: 3 errors from 2 contexts (suppressed: 0 from 0)
_______________________________________________ Lme4-authors mailing list Lme4-authors at lists.r-forge.r-project.org https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/lme4-authors