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.

> if (.Platform$OS.type != "windows") {
+ 
+ library(lme4)
+ library(testthat)
+ 
+ source(system.file("testdata/lme-tst-funs.R", package="lme4", mustWork=TRUE))
+ ##-> gSim(), a general simulation function ...
+ 
+ ## hand-coded Pearson residuals {for  sumFun() }
+ mypresid <- function(x) {
+     mu <- fitted(x)
+     (getME(x,"y") - mu) * sqrt(weights(x)) / sqrt(x@resp$family$variance(mu))
+ }
+ 
+ ## should be equal (up to numerical error) to weights(.,type="working")
+ workingWeights <- function(mod) mod@resp$weights*(mod@resp$muEta()^2)/mod@resp$variance()
+ 
+ ##' Sum of weighted residuals, 4 ways; the last three are identical
+ sumFun <- function(m) {
+     wrss1 <- m@devcomp$cmp["wrss"]
+     wrss2 <- sum(residuals(m,type="pearson")^2)
+     wrss3 <- sum(m@resp$wtres^2)
+     ## compare to hand-fitted Pearson resids ...
+     wrss4 <- sum(mypresid(m)^2)
+     c(wrss1,wrss2,wrss3,wrss4)
+ }
+ ## The relative "error"/differences of the weights w[] entries
+ rel.diff <- function(w) abs(1 - w[-1]/w[1])
+ 
+ set.seed(101)
+ 
+ ## GAMMA
+ g0 <-  glmer(y~x+(1|block),data=gSim(),family=Gamma)
+ expect_true(all(rel.diff(sumFun(g0)) < 1e-13))
+ expect_equal(weights(g0, type = "working"), workingWeights(g0),
+              tolerance = 1e-4)  ## FIXME:  why is such a high tolerance required?
+ 
+ ## BERNOULLI
+ g1 <-  glmer(y~x+(1|block),data=gSim(family=binomial(),nbinom=1),
+              family=binomial)
+ expect_true(all(rel.diff(sumFun(g1)) < 1e-13))
+ expect_equal(weights(g1, type = "working"), workingWeights(g1),
+              tolerance = 1e-5)  ## FIXME:  why is such a high tolerance required?
+ 
+ 
+ ## POISSON
+ (n <- nrow(d.P <- gSim(family=poisson())))
+ g2 <-  glmer(y ~ x + (1|block), data = d.P, family=poisson)
+ g2W <- glmer(y ~ x + (1|block), data = d.P, family=poisson, weights = rep(2,n))
+ expect_true(all(rel.diff(sumFun(g2 )) < 1e-13))
+ expect_true(all(rel.diff(sumFun(g2W)) < 1e-13))
+ ## correct
+ expect_equal(weights(g2, type = "working"), workingWeights(g2),
+              tolerance = 1e-5)  ## FIXME:  why is such a high tolerance required?
+ expect_equal(weights(g2W, type = "working"), workingWeights(g2W),
+              tolerance = 1e-5)  ## FIXME:  why is such a high tolerance required?
+ 
+ 
+ ## non-Bernoulli BINOMIAL
+ g3 <- glmer(y ~ x + (1|block), data= gSim(family=binomial(), nbinom=10),
+             family=binomial)
+ expect_true(all(rel.diff(sumFun(g3)) < 1e-13))
+ expect_equal(weights(g3, type = "working"), workingWeights(g3),
+              tolerance = 1e-4)  ## FIXME:  why is such a high tolerance required?
+ 
+ 
+ 
+ d.b.2 <- gSim(nperblk = 2, family=binomial())
+ g.b.2 <- glmer(y ~ x + (1|block), data=d.b.2, family=binomial)
+ 
+ expect_true(all(rel.diff(sumFun(g.b.2 )) < 1e-13))
+ 
+ 
+ ## Many blocks of only 2 observations each - (but nicely balanced)
+ ## Want this "as" https://github.com/lme4/lme4/issues/47
+ ## (but it "FAILS" survival already):
+ ##
+ ## n2 = n/2 :
+ n2 <- 2048
+ if(FALSE)
+ n2 <-  100 # for building/testing
+ set.seed(47)
+ dB2 <- gSim(n2, nperblk = 2, x= rep(0:1, each= n2), family=binomial())
+ ##                       --  --     ---  --------
+ gB2 <- glmer(y ~ x + (1|block), data=dB2, family=binomial)
+ expect_true(all(rel.diff(sumFun(gB2)) < 1e-13))
+ 
+ ## NB: Finite sample bias of  \hat\sigma_1 and  \hat\beta_1 ("Intercept")
+ ##     tend to zero only slowly for  n2 -> Inf,  e.g., for
+ ## n2 = 2048,  b1 ~= 4.3 (instead of 4);  s1 ~= 1.3 (instead of 1)
+ 
+ ## FAILS -----
+ ## library(survival)
+ ## (gSurv.B2 <- clogit(y ~ x + strata(block), data=dB2))
+ ## ## --> Error in Surv(rep(1, 200L), y) : Time and status are different lengths
+ ## summary(gSurv.B2)
+ ## (SE.surf <- sqrt(diag(vcov(gSurv.B2))))
+ 
+ 
+ 
+ g3 <-  glmer(y ~ x + (1|block),data=gSim(family=binomial(),nbinom=10),
+              family=binomial)
+ expect_equal(var(sumFun(g3)),0)
+ 
+ ## check dispersion parameter
+ ## (lowered tolerance to pass checks on my machine -- SCW)
+ expect_equal(sigma(g0)^2, 0.4888248, tolerance=1e-4)
+ 
+ } ## skip on windows (for speed)
Loading required package: Matrix
=================================================================
==555415==ERROR: AddressSanitizer: heap-use-after-free on address 0x7ea16c5ac430 at pc 0x7bc157bc493c bp 0x7ffdf7b1e410 sp 0x7ffdf7b1e408
READ of size 16 at 0x7ea16c5ac430 thread T0
    #0 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925:46
    #1 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:654:58
    #2 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:668:5
    #3 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:347:16
    #4 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:785:3
    #5 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:954:5
    #6 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890:3
    #7 0x7bc157bc493b 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<!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>>::value, void*>::type) /data/gannet/ripley/R/test-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:858:3
    #8 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:836:3
    #9 0x7bc157bc493b 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-clang/RcppEigen/include/Eigen/src/Core/Assign.h:66:3
    #10 0x7bc157bc493b in Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0>>& 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-clang/RcppEigen/include/Eigen/src/Core/Map.h:162:5
    #11 0x7bc157bc493b in lme4::merPredD::updateXwts(Eigen::Array<double, -1, 1, 0, -1, 1> const&) /data/gannet/ripley/R/packages/tests-clang-ASAN/lme4/src/predModule.cpp:221:30
    #12 0x7bc157b52cd2 in internal_glmerWrkIter(lme4::merPredD*, lme4::glmResp*, bool) /data/gannet/ripley/R/packages/tests-clang-ASAN/lme4/src/external.cpp:274:13
    #13 0x7bc157b52cd2 in pwrssUpdate(lme4::glmResp*, lme4::merPredD*, bool, double, int, int) /data/gannet/ripley/R/packages/tests-clang-ASAN/lme4/src/external.cpp:330:20
    #14 0x7bc157b51bc8 in glmerLaplace /data/gannet/ripley/R/packages/tests-clang-ASAN/lme4/src/external.cpp:386:9
    #15 0x55bfce7d0b06 in R_doDotCall /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:770:11
    #16 0x55bfce894d2d in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8700:21
    #17 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #18 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #19 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #20 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #21 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #22 0x55bfce87bd28 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278:12
    #23 0x55bfce8d4bc5 in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585:8
    #24 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #25 0x55bfce8d2269 in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001:10
    #26 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #27 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #28 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #29 0x55bfce89ceba in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #30 0x55bfce89ceba in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8121:16
    #31 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #32 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #33 0x55bfce87d60a in forcePromise /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976:13
    #34 0x55bfce87b4f5 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1185:6
    #35 0x55bfce87d9a2 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683:12
    #36 0x55bfce87b9d1 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1249:20
    #37 0x55bfce8d4bc5 in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585:8
    #38 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #39 0x55bfce87d9a2 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683:12
    #40 0x55bfce87b9d1 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1249:20
    #41 0x55bfce87d9a2 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683:12
    #42 0x55bfce87b9d1 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1249:20
    #43 0x55bfce87d60a in forcePromise /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976:13
    #44 0x55bfce87b2b7 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1191:2
    #45 0x55bfce88977a in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150:14
    #46 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #47 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #48 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #49 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #50 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #51 0x55bfce87bd28 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278:12
    #52 0x55bfce8d2269 in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001:10
    #53 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #54 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #55 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #56 0x55bfce89ceba in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #57 0x55bfce89ceba in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8121:16
    #58 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #59 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #60 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #61 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #62 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #63 0x55bfce87bd28 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278:12
    #64 0x55bfce731abc in do_docall /data/gannet/ripley/R/svn/R-devel/src/main/coerce.c:2788:12
    #65 0x55bfce88977a in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150:14
    #66 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #67 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #68 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #69 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #70 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #71 0x55bfce87bd28 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278:12
    #72 0x55bfce8d4bc5 in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585:8
    #73 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #74 0x55bfce8d2269 in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001:10
    #75 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #76 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #77 0x55bfce94cbe1 in Rf_ReplIteration /data/gannet/ripley/R/svn/R-devel/src/main/main.c:264:23
    #78 0x55bfce94f410 in R_ReplConsole /data/gannet/ripley/R/svn/R-devel/src/main/main.c:317:11
    #79 0x55bfce94f410 in run_Rmainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1235:5
    #80 0x55bfce94f4a2 in Rf_mainloop /data/gannet/ripley/R/svn/R-devel/src/main/main.c:1242:5
    #81 0x55bfce650d3c in main /data/gannet/ripley/R/svn/R-devel/src/main/Rmain.c:29:5
    #82 0x7fc16f212574 in __libc_start_call_main (/lib64/libc.so.6+0x3574) (BuildId: 92b5376d35bb29c098175948cf3e7cbcae3aeae1)
    #83 0x7fc16f212627 in __libc_start_main@GLIBC_2.2.5 (/lib64/libc.so.6+0x3627) (BuildId: 92b5376d35bb29c098175948cf3e7cbcae3aeae1)
    #84 0x55bfce5687d4 in _start (/data/gannet/ripley/R/clang-ASAN/bin/exec/R+0x17d4)

0x7ea16c5ac430 is located 48 bytes inside of 41648-byte region [0x7ea16c5ac400,0x7ea16c5b66b0)
freed by thread T0 here:
    #0 0x55bfce60d916 in free /home/runner/work/llvm-project/llvm-project/compiler-rt/lib/asan/asan_malloc_linux.cpp:51:3
    #1 0x55bfce96a31b in ReleaseLargeFreeVectors /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:1167:7
    #2 0x55bfce96a31b in RunGenCollect /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:1951:5
    #3 0x55bfce95d18b in R_gc_internal /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:3237:19
    #4 0x55bfce95a195 in Rf_cons /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:2488:2
    #5 0x55bfce830b7b in Rf_defineVar /data/gannet/ripley/R/svn/R-devel/src/main/envir.c:1996:21
    #6 0x55bfce8cd50e in R_execMethod /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2550:5
    #7 0x7fc16ccbdf40 in R_dispatchGeneric /data/gannet/ripley/R/svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154:19
    #8 0x55bfce992853 in do_standardGeneric /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:1348:13
    #9 0x55bfce897a1f in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8085:12
    #10 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #11 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #12 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #13 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #14 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #15 0x55bfce993e97 in R_possible_dispatch /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:1687:10
    #16 0x55bfce8dbd49 in Rf_DispatchOrEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:4191:14
    #17 0x55bfceaafa0a in R_DispatchOrEvalSP /data/gannet/ripley/R/svn/R-devel/src/main/subset.c:671:16
    #18 0x55bfceabcb2f in do_subset3 /data/gannet/ripley/R/svn/R-devel/src/main/subset.c:1264:8
    #19 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #20 0x55bfce87b541 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1219:19
    #21 0x55bfce8d2269 in do_begin /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3001:10
    #22 0x55bfce87b840 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1230:12
    #23 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #24 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #25 0x55bfce89ceba in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #26 0x55bfce89ceba in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8121:16
    #27 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #28 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #29 0x55bfce87d60a in forcePromise /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:976:13
    #30 0x55bfce87b4f5 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1185:6
    #31 0x55bfce87d9a2 in Rf_evalList /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3683:12

previously allocated by thread T0 here:
    #0 0x55bfce60dbb4 in malloc /home/runner/work/llvm-project/llvm-project/compiler-rt/lib/asan/asan_malloc_linux.cpp:67:3
    #1 0x55bfce957c7c in Rf_allocVector3 /data/gannet/ripley/R/svn/R-devel/src/main/memory.c:2894:7
    #2 0x55bfce6b032b in Rf_allocVector /data/gannet/ripley/R/svn/R-devel/src/include/Rinlinedfuns.h:609:12
    #3 0x55bfce6b032b in Rf_allocMatrix /data/gannet/ripley/R/svn/R-devel/src/main/array.c:236:17
    #4 0x7bc15e6d80ce in modelmatrix /data/gannet/ripley/R/svn/R-devel/src/library/stats/src/model.c:676:17
    #5 0x55bfce7c5365 in do_External /data/gannet/ripley/R/svn/R-devel/src/main/dotcode.c:573:11
    #6 0x55bfce88977a in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150:14
    #7 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #8 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #9 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #10 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #11 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #12 0x55bfce989e14 in dispatchMethod /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:473:16
    #13 0x55bfce988a98 in Rf_usemethod /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:523:9
    #14 0x55bfce98b073 in do_usemethod /data/gannet/ripley/R/svn/R-devel/src/main/objects.c:579:9
    #15 0x55bfce88ef93 in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8170:15
    #16 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #17 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #18 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #19 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #20 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #21 0x55bfce87bd28 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278:12
    #22 0x55bfce8d9653 in do_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3975:13
    #23 0x55bfce88977a in bcEval_loop /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:8150:14
    #24 0x55bfce87cd74 in bcEval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:7533:16
    #25 0x55bfce87b301 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1167:8
    #26 0x55bfce8c3a9c in R_execClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2398:39
    #27 0x55bfce8c2c58 in applyClosure_core /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2314:16
    #28 0x55bfce87e3be in Rf_applyClosure /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:2333:16
    #29 0x55bfce87bd28 in Rf_eval /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:1278:12
    #30 0x55bfce8d4bc5 in do_set /data/gannet/ripley/R/svn/R-devel/src/main/eval.c:3585:8

SUMMARY: AddressSanitizer: heap-use-after-free /data/gannet/ripley/R/test-clang/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925:46 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:
  0x7ea16c5ac180: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7ea16c5ac200: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7ea16c5ac280: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7ea16c5ac300: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
  0x7ea16c5ac380: fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa fa
=>0x7ea16c5ac400: fd fd fd fd fd fd[fd]fd fd fd fd fd fd fd fd fd
  0x7ea16c5ac480: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7ea16c5ac500: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7ea16c5ac580: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7ea16c5ac600: fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd fd
  0x7ea16c5ac680: 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
==555415==ABORTING