[Rd] Hard crash of lme4 in R-devel

[Therneau, Terry M., Ph.D.]

I ran across this by accident when working up an example.    It uses a data set from the survival package, but nothing else from there.  Fails on the Intel machine shown below, and on a virtual linux instance on a newer Mac.

Terry

> library(survival)
> library(lme4)
Loading required package: Matrix
> sessionInfo()
R Under development (unstable) (2024-06-14 r86747)
Platform: aarch64-unknown-linux-gnu
Running under: Ubuntu 22.04.4 LTS

Matrix products: default
BLAS:   /usr/local/src/R-devel/lib/libRblas.so
LAPACK: /usr/lib/aarch64-linux-gnu/lapack/liblapack.so.3.10.0

locale:
[1] LC_CTYPE=C.UTF-8       LC_NUMERIC=C           LC_TIME=C.UTF-8
 [4] LC_COLLATE=C           LC_MONETARY=C.UTF-8    LC_MESSAGES=C.UTF-8
 [7] LC_PAPER=C.UTF-8       LC_NAME=C              LC_ADDRESS=C
[10] LC_TELEPHONE=C         LC_MEASUREMENT=C.UTF-8 LC_IDENTIFICATION=C

time zone: America/Chicago
tzcode source: system (glibc)

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] lme4_1.1-35.1  Matrix_1.7-0   survival_3.6-6

loaded via a namespace (and not attached):
[1] minqa_1.2.5    MASS_7.3-60.2  compiler_4.5.0 Rcpp_1.0.11    splines_4.5.0
 [6] nlme_3.1-164   grid_4.5.0     nloptr_2.0.3   boot_1.3-30    lattice_0.22-6
>
> # Create a data set of all PBC subjects with 4 or more values
> count <- table(pbcseq$id)
> data4 <- subset(pbcseq, id %in% names(count)[count>3])
> data4$year <- data4$day/365.25
>
> fit <- lmer(log(bili) ~ year + (1+year | id), data4)

-- more lines than I can capture --

= NULL, verbose = 0L,     control = list(optimizer = "nloptwrap", restart_edge = TRUE,         boundary.tol = 1e-05, calc.derivs = TRUE, use.last.params = FALSE,         checkControl = list(check.nobs.vs.rankZ = "ignore", check.nobs.vs.nlev = "stop",             check.nlev.gtreq.5 = "ignore", check.nlev.gtr.1 = "stop",             check.nobs.vs.nRE = "stop", check.rankX = "message+drop.cols",             check.scaleX = "warning", check.formula.LHS = "stop"),         checkConv = list(check.conv.grad = list(action = "warning",             tol = 0.002, relTol = NULL), check.conv.singular = list(            action = "message", tol = 1e-04), check.conv.hess = list(            action = "warning", tol = 1e-06)), optCtrl = list()))
10: do.call(mkLmerDevfun, c(lmod, list(start = start, verbose = verbose,     control = control)))
11: lmer(log(bili) ~ year + (1 + year | id), data4)

Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace
Selection:


	[[alternative HTML version deleted]]