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]]

