R Under development (unstable) (2026-07-04 r90207) -- "Unsuffered Consequences"
Copyright (C) 2026 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu

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.

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.

> library("lme4")
Loading required package: Matrix
> L <- load(system.file("testdata","crabs_randdata2.Rda",package="lme4"))
> ## randdata0: simulated data, in form suitable for plotting
> ## randdata: simulated data, in form suitable for analysis
> 
> ## fr  ## alive/dead formula
> ## fr2 ## proportion alive formula (use with weights=initial.snail.density)
> 
> ## FIXME: there are still bigger differences than I'd like between the approaches
> ## (mostly in the random-effects correlation).  It's not clear who's right;
> ## lme4 thinks its parameters are better, but ??  Could be explored further.
> 
> if (FALSE) {
+ ## library(ggplot2)  ## commented to avoid triggering Suggests: requirement
+ library(grid)
+ zmargin <- theme(panel.margin=unit(0,"lines"))
+ theme_set(theme_bw())
+ g1 <- ggplot(randdata0,aes(x=snail.size,y=surv,colour=snail.size,fill=snail.size))+
+     geom_hline(yintercept=1,colour="black")+
+     stat_sum(aes(size=factor(..n..)),alpha=0.6)+
+     facet_grid(.~ttt)+zmargin+
+     geom_boxplot(fill=NA,outlier.colour=NULL,outlier.shape=3)+  ## set outliers to same colour as points
+     ## (hard to see which are outliers, but it doesn't really matter in this case)
+     scale_size_discrete("# obs",range=c(2,5))
+ }
> 
> if (.Platform$OS.type != "windows") {
+ t1 <- system.time(glmer1 <- glmer(fr2,weights=initial.snail.density,
+                                   family ="binomial", data=randdata))
+ t1B <- system.time(glmer1B <- glmer(fr,family ="binomial", data=randdata))
+ 
+ res1 <- c(fixef(glmer1),c(VarCorr(glmer1)$plot))
+ res1B <- c(fixef(glmer1B),c(VarCorr(glmer1B)$plot))
+ p1 <- unlist(getME(glmer1,c("theta","beta")))
+ stopifnot(all.equal(res1,res1B))
+ 
+ dfun <- update(glmer1,devFunOnly=TRUE)
+ stopifnot(all.equal(dfun(p1),c(-2*logLik(glmer1))))
+ ##
+ ## library(lme4.0)  ## version 0.999999.2 results
+ ## t1_lme4.0 <- system.time(glmer1X <-
+ ##                           glmer(fr2,weights=initial.snail.density,
+ ##                                 family ="binomial", data=randdata))
+ ## dput(c(fixef(glmer1X),c(VarCorr(glmer1X)$plot)))
+ ## p1X <- c(getME(glmer1X,"theta"),getME(glmer1X,"beta"))
+ p1X <- c(0.681301656652347, -1.14775239687404, 0.436143018123226, 
+ 2.77730476938968, 0.609023583738824, -1.60055813739844, 2.0324468778545, 
+ 0.624173873057839, -1.7908793509579, -2.44540201631615, -1.42365990002708, 
+ -2.26780929006268, 0.700928084600075, -1.26220238391029, 0.369024582097804, 
+ 3.44325347343035, 2.26400391093108)
+ stopifnot(all.equal(unname(p1),p1X,tolerance=0.03))
+ dfun(p1X)
+ dfun(p1)
+ ## ~ 1.8 seconds elapsed time
+ lme4.0_res <- 
+     structure(c(2.77730476938968, 0.609023583738824, -1.60055813739844, 
+                 2.0324468778545, 0.624173873057839, -1.7908793509579, -2.44540201631615, 
+                 -1.42365990002708, -2.26780929006268, 0.700928084600075, -1.26220238391029, 
+                 0.369024582097804, 3.44325347343035, 2.26400391093108, 0.464171947357232, 
+                 -0.532754465140956, -0.532754465140956, 0.801690946568518),
+           .Names = c("(Intercept)", 
+           "crab.speciesS", "crab.speciesW", "crab.sizeS", "crab.sizeM", 
+           "snail.sizeS", "crab.speciesS:crab.sizeS", "crab.speciesS:crab.sizeM", 
+           "crab.speciesS:snail.sizeS", "crab.speciesW:snail.sizeS",
+           "crab.sizeS:snail.sizeS", "crab.sizeM:snail.sizeS",
+           "crab.speciesS:crab.sizeS:snail.sizeS", 
+           "crab.speciesS:crab.sizeM:snail.sizeS", "", "", "", ""))
+ 
+ stopifnot(all.equal(res1,lme4.0_res,tolerance=0.015))
+ 
+ 
+ ## library("glmmADMB")
+ ## prop/weights formulation: ~ 7 seconds
+ ## t1_glmmadmb <- system.time(glmer1B <- glmmadmb(fr,family ="binomial",
+ ##                                               corStruct="full",data=randdata))
+ ## dput(c(fixef(glmer1B),c(VarCorr(glmer1B)$plot)))
+ glmmADMB_res <- structure(c(2.7773101267224, 0.609026276823218, -1.60055704634712, 
+ 2.03244174458562, 0.624171008585953, -1.79088398816641, -2.44540300134182, 
+ -1.42366043619683, -2.26780858382505, 0.700927141726545, -1.26219964572264, 
+ 0.369029052442189, 3.44326297908383, 2.26403738918967, 0.46417, 
+ -0.53253, -0.53253, 0.80169), .Names = c("(Intercept)", "crab.speciesS", 
+ "crab.speciesW", "crab.sizeS", "crab.sizeM", "snail.sizeS", "crab.speciesS:crab.sizeS", 
+ "crab.speciesS:crab.sizeM", "crab.speciesS:snail.sizeS", "crab.speciesW:snail.sizeS", 
+ "crab.sizeS:snail.sizeS", "crab.sizeM:snail.sizeS", "crab.speciesS:crab.sizeS:snail.sizeS", 
+ "crab.speciesS:crab.sizeM:snail.sizeS", "", "", "", ""))
+ 
+ stopifnot(all.equal(res1B,glmmADMB_res,tolerance=0.015))
+ } ## skip on windows (for speed)
=================================================================
==4048105==ERROR: AddressSanitizer: heap-use-after-free on address 0x7e16d0c3e130 at pc 0x7bd6ba1d0d9d bp 0x7ffe9fab01e0 sp 0x7ffe9fab01d8
READ of size 8 at 0x7e16d0c3e130 thread T0
    #0 0x7bd6ba1d0d9c in Eigen::internal::product_evaluator<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, 8, Eigen::DiagonalShape, Eigen::DenseShape, double, double>::coeff(long, long) const /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925
    #1 0x7bd6ba1d0d9c in Eigen::internal::generic_dense_assignment_kernel<Eigen::internal::evaluator<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > >, Eigen::internal::evaluator<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >, Eigen::internal::assign_op<double, double>, 0>::assignCoeff(long, long) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:654
    #2 0x7bd6ba1d0d9c in Eigen::internal::generic_dense_assignment_kernel<Eigen::internal::evaluator<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > >, Eigen::internal::evaluator<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >, Eigen::internal::assign_op<double, double>, 0>::assignCoeffByOuterInner(long, long) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:668
    #3 0x7bd6ba1d0d9c in Eigen::internal::dense_assignment_loop<Eigen::internal::generic_dense_assignment_kernel<Eigen::internal::evaluator<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > >, Eigen::internal::evaluator<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >, Eigen::internal::assign_op<double, double>, 0>, 0, 0>::run(Eigen::internal::generic_dense_assignment_kernel<Eigen::internal::evaluator<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > >, Eigen::internal::evaluator<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >, Eigen::internal::assign_op<double, double>, 0>&) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:347
    #4 0x7bd6ba1d0d9c in void Eigen::internal::call_dense_assignment_loop<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, Eigen::internal::assign_op<double, double> >(Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >&, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> const&, Eigen::internal::assign_op<double, double> const&) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:785
    #5 0x7bd6ba1d0d9c in Eigen::internal::Assignment<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, Eigen::internal::assign_op<double, double>, Eigen::internal::Dense2Dense, void>::run(Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >&, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> const&, Eigen::internal::assign_op<double, double> const&) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:954
    #6 0x7bd6ba1d0d9c in void Eigen::internal::call_assignment_no_alias<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, Eigen::internal::assign_op<double, double> >(Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >&, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> const&, Eigen::internal::assign_op<double, double> const&) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890
    #7 0x7bd6ba1d0d9c in void Eigen::internal::call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, Eigen::internal::assign_op<double, double> >(Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >&, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> const&, Eigen::internal::assign_op<double, double> const&, Eigen::internal::enable_if<!Eigen::internal::evaluator_assume_aliasing<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, Eigen::internal::evaluator_traits<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >::Shape>::value, void*>::type) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:858
    #8 0x7bd6ba1d0d9c in void Eigen::internal::call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >(Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >&, Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> const&) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:836
    #9 0x7bd6ba1d0d9c in Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >& Eigen::MatrixBase<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > >::operator=<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> >(Eigen::DenseBase<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1> > const&) /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/Assign.h:66
    #10 0x7bd6ba1d0d9c in lme4::merPredD::updateXwts(Eigen::Array<double, -1, 1, 0, -1, 1> const&) /data/gannet/ripley/R/packages/tests-gcc-SAN/lme4/src/predModule.cpp:221
    #11 0x7bd6ba046416 in internal_glmerWrkIter /data/gannet/ripley/R/packages/tests-gcc-SAN/lme4/src/external.cpp:274
    #12 0x7bd6ba046416 in pwrssUpdate /data/gannet/ripley/R/packages/tests-gcc-SAN/lme4/src/external.cpp:330
    #13 0x7bd6ba05e3d4 in glmerLaplace /data/gannet/ripley/R/packages/tests-gcc-SAN/lme4/src/external.cpp:386
    #14 0x00000075005e in R_doDotCall /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:770
    #15 0x0000008ed70a in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8700
    #16 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #17 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #18 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #19 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #20 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #21 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #22 0x000000894a2e in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585
    #23 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #24 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001
    #25 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #26 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #27 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #28 0x00000088038b in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #29 0x0000008eb805 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8121
    #30 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #31 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #32 0x000000930ad1 in forcePromise.part.0.lto_priv.0 /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976
    #33 0x00000086405c in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1185
    #34 0x0000008952f4 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683
    #35 0x00000086388c in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1249
    #36 0x000000894a2e in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585
    #37 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #38 0x0000008952f4 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683
    #39 0x00000086388c in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1249
    #40 0x0000008952f4 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683
    #41 0x00000086388c in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1249
    #42 0x000000930ad1 in forcePromise.part.0.lto_priv.0 /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976
    #43 0x000000863f37 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1191
    #44 0x0000008d6d02 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150
    #45 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #46 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #47 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #48 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #49 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #50 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #51 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001
    #52 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #53 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #54 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #55 0x00000088038b in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #56 0x0000008eb805 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8121
    #57 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #58 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #59 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #60 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #61 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #62 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #63 0x0000006237b5 in do_docall /data/gannet/ripley/R/svn/R-devel/src/main/coerce.c:2788
    #64 0x0000008d6d02 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150
    #65 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #66 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #67 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #68 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #69 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #70 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #71 0x000000894a2e in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585
    #72 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #73 0x000000930ad1 in forcePromise.part.0.lto_priv.0 /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976
    #74 0x0000008a9bb1 in forcePromise /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:956
    #75 0x0000008a9bb1 in getvar /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:5867
    #76 0x0000008e6b7c in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7880
    #77 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #78 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #79 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #80 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #81 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #82 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #83 0x000000894a2e in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585
    #84 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #85 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001
    #86 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #87 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #88 0x000000a0b688 in Rf_ReplIteration /data/gannet/ripley/R/svn/R-devel/src/main/main.c:264
    #89 0x000000a0b688 in R_ReplConsole /data/gannet/ripley/R/svn/R-devel/src/main/main.c:317
    #90 0x000000a1a97a in run_Rmainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1235
    #91 0x000000a1aa12 in Rf_mainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1242
    #92 0x00000041306f in main /data/gannet/ripley/R/svn/R-devel/src/main/Rmain.c:29
    #93 0x7fd6d4012574 in __libc_start_call_main (/lib64/libc.so.6+0x3574) (BuildId: 92b5376d35bb29c098175948cf3e7cbcae3aeae1)
    #94 0x7fd6d4012627 in __libc_start_main@GLIBC_2.2.5 (/lib64/libc.so.6+0x3627) (BuildId: 92b5376d35bb29c098175948cf3e7cbcae3aeae1)
    #95 0x000000413a14 in _start (/data/gannet/ripley/R/gcc-SAN3/bin/exec/R+0x413a14) (BuildId: f060bf4ee83b802127121439cefaeb137652c0ab)

0x7e16d0c3e130 is located 48 bytes inside of 7664-byte region [0x7e16d0c3e100,0x7e16d0c3fef0)
freed by thread T0 here:
    #0 0x7fd6d56e5beb in free.part.0 (/lib64/libasan.so.8+0xe5beb) (BuildId: d3cb6206dff19da52969c009f4cd93611901c478)
    #1 0x000000a1bf9f in ReleaseLargeFreeVectors /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:1167
    #2 0x000000a49e0d in RunGenCollect /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:1951
    #3 0x000000a7244d in R_gc_internal /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:3237
    #4 0x000000a872bb in CONS_NR /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:2519
    #5 0x0000008e18ea in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8059
    #6 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #7 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #8 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #9 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #10 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #11 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #12 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #13 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001
    #14 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #15 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #16 0x000000882ecd in R_execMethod /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2567
    #17 0x7fd6d0d28e3a in R_dispatchGeneric /data/gannet/ripley/R/svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154
    #18 0x000000ab5d56 in do_standardGeneric /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:1348
    #19 0x0000008f0311 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8085
    #20 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #21 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #22 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #23 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #24 0x000000ab91dc in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #25 0x000000ab91dc in R_possible_dispatch /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:1687
    #26 0x00000089b0be in tryDispatch /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:6090
    #27 0x0000008ce28f in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8320
    #28 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #29 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #30 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #31 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314

previously allocated by thread T0 here:
    #0 0x7fd6d56e6f2b in malloc (/lib64/libasan.so.8+0xe6f2b) (BuildId: d3cb6206dff19da52969c009f4cd93611901c478)
    #1 0x000000a7a3e0 in Rf_allocVector3 /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:2894
    #2 0x0000005077b3 in Rf_allocVector /data/gannet/ripley/R/svn/R-devel/src/include/Rinlinedfuns.h:609
    #3 0x0000005077b3 in Rf_allocMatrix /data/gannet/ripley/R/svn/R-devel/src/main/array.c:236
    #4 0x7bd6c294b36e in modelmatrix /data/gannet/ripley/R/svn/R-devel/src/library/stats/src/model.c:676
    #5 0x00000073b52a in do_External /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:573
    #6 0x0000008d6d02 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150
    #7 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #8 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #9 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #10 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #11 0x000000a9f646 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #12 0x000000a9f646 in applyMethod /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:120
    #13 0x000000aa3480 in dispatchMethod /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:473
    #14 0x000000aa56c0 in Rf_usemethod.isra.0 /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:523
    #15 0x000000aa6026 in do_usemethod /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:579
    #16 0x0000008e3f11 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8170
    #17 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #18 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #19 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #20 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #21 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #22 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #23 0x00000089c5a1 in do_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3975
    #24 0x0000008d6d02 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150
    #25 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #26 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #27 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #28 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #29 0x000000863098 in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #30 0x000000863098 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278
    #31 0x000000894a2e in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585
    #32 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #33 0x000000930ad1 in forcePromise.part.0.lto_priv.0 /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976

SUMMARY: AddressSanitizer: heap-use-after-free /data/gannet/ripley/R/test-dev/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925 in Eigen::internal::product_evaluator<Eigen::Product<Eigen::DiagonalWrapper<Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> > const>, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, 1>, 8, Eigen::DiagonalShape, Eigen::DenseShape, double, double>::coeff(long, long) const
Shadow bytes around the buggy address:
  0x7e16d0c3de80: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7e16d0c3df00: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7e16d0c3df80: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7e16d0c3e000: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7e16d0c3e080: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
=>0x7e16d0c3e100: fd fd fd fd fd fd[fd]fd fd fd fd fd fd fd fd fd
  0x7e16d0c3e180: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7e16d0c3e200: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7e16d0c3e280: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7e16d0c3e300: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7e16d0c3e380: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
Shadow byte legend (one shadow byte represents 8 application bytes):
  Addressable:           00
  Partially addressable: 01 02 03 04 05 06 07 
  Heap left redzone:       fa
  Freed heap region:       fd
  Stack left redzone:      f1
  Stack mid redzone:       f2
  Stack right redzone:     f3
  Stack after return:      f5
  Stack use after scope:   f8
  Global redzone:          f9
  Global init order:       f6
  Poisoned by user:        f7
  Container overflow:      fc
  Array cookie:            ac
  Intra object redzone:    bb
  ASan internal:           fe
  Left alloca redzone:     ca
  Right alloca redzone:    cb
==4048105==ABORTING