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.

> ## require(agridat)
> ## dat <- gotway.hessianfly
> 
> if (.Platform$OS.type != "windows") {
+ ## don't actually use gotway_hessianfly_fit or gotway_hessianfly_prof,
+ ## so we should be OK even with R< 3.0.1
+ load(system.file("testdata","gotway_hessianfly.rda",package="lme4"))
+ # Block random.  See Glimmix manual, output 1.18.
+ # Note: (Different parameterization)
+ 
+ ## require("lme4.0")
+ ## fit2 <- glmer(cbind(y, n-y) ~ gen + (1|block), data=dat, family=binomial)
+ ## params <- list(fixef=fixef(fit2),theta=getME(fit2,"theta"))
+ ## detach("package:lme4.0")
+ lme4.0fit <- structure(list(fixef = structure(c(1.50345713031203, -0.193853259383803, 
+ -0.540808391060274, -1.43419379979154, -0.203701042949808, -0.978322555343941, 
+ -0.604078624475678, -1.67742449813309, -1.39842466673692, -0.681709344788684, 
+ -1.46295367186169, -1.45908310198959, -3.55285756517073, -2.50731975980307, 
+ -2.08716296677356, -2.96974270029992), .Names = c("(Intercept)", 
+ "genG02", "genG03", "genG04", "genG05", "genG06", "genG07", "genG08", 
+ "genG09", "genG10", "genG11", "genG12", "genG13", "genG14", "genG15", 
+ "genG16")), theta = structure(0.0319087494293615, .Names = "block.(Intercept)")), .Names = c("fixef", 
+ "theta"))
+ 
+ ## start doesn't work because we don't get there
+ library(lme4)
+ m1 <- glmer(cbind(y, n-y) ~ gen + (1|block), data=gotway.hessianfly,
+             family=binomial)
+ m1B <- update(m1,control=glmerControl(optimizer="bobyqa"))
+ max(abs(m1@optinfo$derivs$gradient))  ## 0.0012
+ max(abs(m1B@optinfo$derivs$gradient)) ## 2.03e-5
+ abs(m1@optinfo$derivs$gradient)/abs(m1B@optinfo$derivs$gradient)
+ ## bobyqa gets gradients *at least* 1.64* lower
+ 
+ lme4fit <- list(fixef=fixef(m1),theta=getME(m1,"theta"))
+ 
+ ## hack around slight naming differences
+ lme4fit$theta <- unname(lme4fit$theta)
+ lme4.0fit$theta <- unname(lme4.0fit$theta)
+ ## difference in theta on x86_64-w64-mingw32 (64-bit) with r-devel is  0.000469576
+ stopifnot(all.equal(lme4fit, lme4.0fit, tolerance = 5e-4))
+ 
+ ## Fun stuff: visualize and alternative model
+ 
+ ## library(ggplot2)
+ ## dat$prop <- dat$y/dat$n
+ ## theme_set(theme_bw())
+ ## ggplot(dat,aes(x=gen,y=prop,colour=block))+geom_point(aes(size=n))+
+ ##     geom_line(aes(group=block,colour=block))+
+ ##     geom_smooth(family=binomial,aes(weight=n,colour=block,group=block),method="glm",
+ ##                 alpha=0.1)
+ 
+ ## dat$obs <- factor(seq(nrow(dat)))
+ ## m2 <- glmer(cbind(y, n-y) ~ block+ (1|gen) + (1|obs), data=dat, family=binomial)
+ 
+ }  ## not on windows/CRAN
Loading required package: Matrix
=================================================================
==3819436==ERROR: AddressSanitizer: heap-use-after-free on address 0x7df7a44c0930 at pc 0x7ba7801d0d9d bp 0x7ffce8b7e8f0 sp 0x7ffce8b7e8e8
READ of size 8 at 0x7df7a44c0930 thread T0
    #0 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba7801d0d9c 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 0x7ba780046416 in internal_glmerWrkIter /data/gannet/ripley/R/packages/tests-gcc-SAN/lme4/src/external.cpp:274
    #12 0x7ba780046416 in pwrssUpdate /data/gannet/ripley/R/packages/tests-gcc-SAN/lme4/src/external.cpp:330
    #13 0x7ba78005e3d4 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 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001
    #74 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #75 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #76 0x000000a0b688 in Rf_ReplIteration /data/gannet/ripley/R/svn/R-devel/src/main/main.c:264
    #77 0x000000a0b688 in R_ReplConsole /data/gannet/ripley/R/svn/R-devel/src/main/main.c:317
    #78 0x000000a1a97a in run_Rmainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1235
    #79 0x000000a1aa12 in Rf_mainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1242
    #80 0x00000041306f in main /data/gannet/ripley/R/svn/R-devel/src/main/Rmain.c:29
    #81 0x7fa79a012574 in __libc_start_call_main (/lib64/libc.so.6+0x3574) (BuildId: 92b5376d35bb29c098175948cf3e7cbcae3aeae1)
    #82 0x7fa79a012627 in __libc_start_main@GLIBC_2.2.5 (/lib64/libc.so.6+0x3627) (BuildId: 92b5376d35bb29c098175948cf3e7cbcae3aeae1)
    #83 0x000000413a14 in _start (/data/gannet/ripley/R/gcc-SAN3/bin/exec/R+0x413a14) (BuildId: f060bf4ee83b802127121439cefaeb137652c0ab)

0x7df7a44c0930 is located 48 bytes inside of 8240-byte region [0x7df7a44c0900,0x7df7a44c2930)
freed by thread T0 here:
    #0 0x7fa79b8e5beb 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 0x000000a782fb in Rf_cons /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:2488
    #5 0x00000080993d in Rf_defineVar /data/gannet/ripley/R/svn/R-devel/src/main/envir.c:1996
    #6 0x000000882d02 in R_execMethod /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2544
    #7 0x7fa79712be3a in R_dispatchGeneric /data/gannet/ripley/R/svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154
    #8 0x000000ab5d56 in do_standardGeneric /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:1348
    #9 0x0000008f0311 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8085
    #10 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #11 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #12 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #13 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #14 0x000000ab91dc in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #15 0x000000ab91dc in R_possible_dispatch /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:1687
    #16 0x000000899880 in Rf_DispatchOrEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:4191
    #17 0x000000d0bd6b in R_DispatchOrEvalSP /data/gannet/ripley/R/svn/R-devel/src/main/subassign.c:1564
    #18 0x000000d6ba7f in do_subset3 /data/gannet/ripley/R/svn/R-devel/src/main/subset.c:1264
    #19 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #20 0x000000862f2e in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1219
    #21 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001
    #22 0x0000008634c6 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230
    #23 0x00000087ad2a in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398
    #24 0x00000087ecdd in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314
    #25 0x00000088038b in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333
    #26 0x0000008eb805 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8121
    #27 0x0000008c3f63 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533
    #28 0x0000008629ea in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167
    #29 0x000000930ad1 in forcePromise.part.0.lto_priv.0 /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976
    #30 0x00000086405c in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1185

previously allocated by thread T0 here:
    #0 0x7fa79b8e6f2b 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 0x7ba78894b36e 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 0x00000087729b in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001

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:
  0x7df7a44c0680: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7df7a44c0700: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7df7a44c0780: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7df7a44c0800: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7df7a44c0880: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
=>0x7df7a44c0900: fd fd fd fd fd fd[fd]fd fd fd fd fd fd fd fd fd
  0x7df7a44c0980: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7df7a44c0a00: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7df7a44c0a80: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7df7a44c0b00: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7df7a44c0b80: 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
==3819436==ABORTING