getME function crashes
Johannes, on Linux-64, with the same version of R and the same version of all loaded packages, your example code runs fine. Apart from Windows-32, the only difference I've spotted is that a module named "tools_3.0.2" does not appear in my sessionInfo(), but that might be platform dependant. Hugo On Thu, 7 Nov 2013 17:31:41 +0100
Johannes Radinger <johannesradinger at gmail.com> wrote:
Hi, Actually I tried to use the sim function from the package "arm" to get posterior simulations of beta from a merMod. But somehow this function crashed and I could trace the problem back to the function getME(model,"L") which is the cause for the crash. For example here an example dataset for a binomial model (glmer) from the sim function in the 'arm'-package: set.seed (1) J <- 15 n <- J*(J+1)/2 group <- rep (1:J, 1:J) mu.a <- 5 sigma.a <- 2 a <- rnorm (J, mu.a, sigma.a) b <- -3 x <- rnorm (n, 2, 1) y <- rbinom (n, 1, invlogit (a[group] + b*x)) mod <- glmer (y ~ x + (1 | group), family=binomial(link="logit")) When I now run getME(mod, "L") my R session completely crashes with following error: *** caught segfault *** address (nil), cause 'memory not mapped' Traceback: 1: .Call(merPredDL, ptr()) 2: PR$L() 3: getME(mod, "L") Possible actions: 1: abort (with core dump, if enabled) 2: normal R exit 3: exit R without saving workspace 4: exit R saving workspace I am running following setup: R version 3.0.2 (2013-09-25) Platform: i686-pc-linux-gnu (32-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] arm_1.6-09 lme4_1.0-5 lattice_0.20-24 Matrix_1.1-0 MASS_7.3-29 loaded via a namespace (and not attached): [1] abind_1.4-0 coda_0.16-1 grid_3.0.2 minqa_1.2.1 nlme_3.1-111 splines_3.0.2 tools_3.0.2 Can anyone reproduce this problem? /Johannes [[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
--