==3511189== Memcheck, a memory error detector
==3511189== Copyright (C) 2002-2024, and GNU GPL'd, by Julian Seward et al.
==3511189== Using Valgrind-3.24.0 and LibVEX; rerun with -h for copyright info
==3511189== Command: /data/blackswan/ripley/R/R-devel-vg/bin/exec/R -f glmer-1.R --restore --save --no-readline --vanilla
==3511189== 

R Under development (unstable) (2026-07-02 r90204) -- "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 (lme4:::testLevel() > 1 || .Platform$OS.type!="windows") withAutoprint({
+ 
+     ## generalized linear mixed model
+     stopifnot(suppressPackageStartupMessages(require(lme4)))
+     options(show.signif.stars = FALSE)
+ 
+     source(system.file("test-tools-1.R", package = "Matrix"), keep.source = FALSE)
+     ##
+     ##' Check that coefficient +- "2" * SD  contains true value
+     ##'
+     ##' @title Check that confidence interval for coefficients contains true value
+     ##' @param fm fitted model, e.g., from  lm(), lmer(), glmer(), ..
+     ##' @param true.coef numeric vector of true (fixed effect) coefficients
+     ##' @param conf.level confidence level for confidence interval
+     ##' @param sd.factor the "2", i.e. default 1.96 factor for the confidence interval
+     ##' @return TRUE or a string of "error"
+     ##' @author Martin Maechler
+     chkFixed <- function(fm, true.coef, conf.level = 0.95,
+                          sd.factor = qnorm((1+conf.level)/2))
+     {
+         stopifnot(is.matrix(cf <- coefficients(summary(fm))), ncol(cf) >= 2)
+         cc <- cf[,1]
+         sd <- cf[,2]
+         if(any(out1 <- true.coef < cc - sd.factor*sd))
+             return(sprintf("true coefficient[j], j=%s, is smaller than lower confidence limit",
+                            paste(which(out1), collapse=", ")))
+         if(any(out2 <- true.coef > cc + sd.factor*sd))
+             return(sprintf("true coefficient[j], j=%s, is larger than upper confidence limit",
+                            paste(which(out2), collapse=", ")))
+         ## else, return
+         TRUE
+     }
+ 
+ 
+     ## TODO: (1) move these to ./glmer-ex.R [DONE]
+     ## ----  (2) "rationalize" with ../man/cbpp.Rd
+                                         #m1e <- glmer1(cbind(incidence, size - incidence) ~ period + (1 | herd),
+                                         #              family = binomial, data = cbpp, doFit = FALSE)
+     ## now
+                                         #bobyqa(m1e, control = list(iprint = 2L))
+ 
+     m1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
+                 family = binomial, data = cbpp)
+     m1. <- update(m1, start = getME(m1, c("theta", "fixef")))
+     dm1 <- drop1(m1)
+     stopifnot(all.equal(drop1(m1.), dm1, tol = 1e-10))# Lnx(F28) 64b: 4e-12
+     ## response as a vector of probabilities and usage of argument "weights"
+     m1p <- glmer(incidence / size ~ period + (1 | herd), weights = size,
+                  family = binomial, data = cbpp)
+     ## Confirm that these are equivalent:
+     stopifnot(all.equal(fixef(m1), fixef(m1p)),
+               all.equal(ranef(m1), ranef(m1p)),
+               TRUE)
+     ## for(m in c(m1, m1p)) {
+     ##     cat("-------\\n\\nCall: ",
+     ##         paste(format(getCall(m)), collapse="\\n"), "\\n")
+     ##     print(logLik(m)); cat("AIC:", AIC(m), "\\n") ; cat("BIC:", BIC(m),"\\n")
+     ## }
+     stopifnot(all.equal(logLik(m1), logLik(m1p)),
+               all.equal(AIC(m1),    AIC(m1p)),
+               all.equal(BIC(m1),    BIC(m1p)))
+ 
+ 
+     ## changed tolPwrss to 1e-7 to match other default
+     m1b <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
+                  family = binomial, data = cbpp, verbose = 2L,
+                  control =
+                      glmerControl(optimizer="bobyqa", tolPwrss=1e-7,
+                                   optCtrl=list(rhobeg=0.2, rhoend=2e-7)))
+ 
+     ## using nAGQ=9L provides a better evaluation of the deviance
+     m.9 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
+                  family = binomial, data = cbpp, nAGQ = 9)
+ 
+     ## check with nAGQ = 25
+     m2 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
+                 family = binomial, data = cbpp, nAGQ = 25)
+ 
+     ## loosened tolerance on parameters
+     stopifnot(is((cm2 <- coef(m2)), "coef.mer"),
+               dim(cm2$herd) == c(15,4),
+               all.equal(fixef(m2),
+ ### lme4a [from an Ubuntu 11.10 amd64 system]
+                         c(-1.39922533406847, -0.991407294757321,
+                           -1.12782184600404, -1.57946627431248),
+                         ##c(-1.3766013, -1.0058773,
+                         ##  -1.1430128, -1.5922817),
+                         tolerance = 5.e-4,
+                         check.attributes=FALSE),
+               all.equal(c(-2*logLik(m2)), 100.010030538022, tolerance=1e-9),
+               all.equal(deviance(m2), 73.373, tolerance=1e-5)
+               ## with bobyqa first (AGQ=0), then
+               ##all.equal(deviance(m2), 101.119749563, tolerance=1e-9)
+               )
+ 
+     ## 32-bit Ubuntu 10.04:
+     coef_m1_lme4.0 <- structure(c(-1.39853505102576,
+                                   -0.992334712470269, -1.12867541092127,
+                                   -1.58037389566025),
+                                 .Names = c("(Intercept)", "period2", "period3",
+                                            "period4"))
+ 
+     ## library(glmmADMB)
+     ## mg <- glmmadmb(cbind(incidence, size - incidence) ~ period + (1 | herd),
+     ##                family = "binomial", data = cbpp)
+     coef_m1_glmmadmb <- structure(c(-1.39853810064827, -0.99233330126975, -1.12867317840779,
+                                     -1.58031150854503), .Names = c("(Intercept)", "period2", "period3",
+                                                                    "period4"))
+ 
+     ## library(glmmML)
+     ## mm <- glmmML(cbind(incidence, size - incidence) ~ period,
+     ##              cluster=herd,
+     ##             family = "binomial", data = cbpp)
+     coef_m1_glmmML <- structure(c(-1.39853234657711, -0.992336901732793, -1.12867036466201,
+                                   -1.58030977686564), .Names = c("(Intercept)", "period2", "period3",
+                                                                  "period4"))
+ 
+     ## lme4[r 1636], 64-bit ubuntu 11.10:
+     ## c(-1.3788385, -1.0589543,
+     ##                      -1.1936382, -1.6306271),
+ 
+     stopifnot(is((cm1 <- coef(m1b)), "coef.mer"),
+               dim(cm1$herd) == c(15,4),
+               all.equal(fixef(m1b),fixef(m1),tolerance=4e-5),
+               is.all.equal4(fixef(m1b),
+                             coef_m1_glmmadmb,
+                             coef_m1_lme4.0,
+                             coef_m1_glmmML,
+                             tol = 5e-4)
+               )
+ 
+ 
+     ## Deviance for the new algorithm is lower, eventually we should change the previous test
+     ##stopifnot(deviance(m1) <= deviance(m1e))
+ 
+     showProc.time() #
+ 
+     if (require('MASS', quietly = TRUE)) {
+         bacteria$wk2 <- bacteria$week > 2
+         contrasts(bacteria$trt) <-
+             structure(contr.sdif(3),
+                       dimnames = list(NULL, c("diag", "encourage")))
+         print(fm5 <- glmer(y ~ trt + wk2 + (1|ID),
+                            data=bacteria, family=binomial))
+         showProc.time() #
+ 
+         stopifnot(
+             all.equal(logLik(fm5),
+                       ## was	  -96.127838
+                       structure(-96.13069, nobs = 220L, nall = 220L,
+                                 df = 5L, REML = FALSE,
+                                 class = "logLik"),
+                       tolerance = 5e-4, check.attributes = FALSE)
+            ,
+             all.equal(fixef(fm5),
+                       ## was		 2.834218798		 -1.367099481
+                       c("(Intercept)"= 2.831609490, "trtdiag"= -1.366722631,
+                         ## now	 0.5842291915,		 -1.599148773
+                         "trtencourage"=0.5840147802, "wk2TRUE"=-1.598591346),
+                       tolerance = 1e-4 )
+         )
+     }
+ 
+     ## Failure to specify a random effects term - used to give an obscure message
+     ## Ensure *NON*-translated message; works on Linux,... :
+     if(.Platform$OS.type == "unix") {
+         Sys.setlocale("LC_MESSAGES", "C")
+         tc <- tryCatch(
+             m2 <- glmer(incidence / size ~ period, weights = size,
+                         family = binomial, data = cbpp)
+           , error = function(.) .)
+         stopifnot(inherits(tc, "error"),
+                   identical(tc$message,
+                             "No random effects terms specified in formula"))
+     }
+ 
+ 
+     ## glmer - Modeling overdispersion as "mixture" aka
+     ## ----- - *ONE* random effect *PER OBSERVATION" -- example inspired by Ben Bolker:
+ 
+     ##' <description>
+     ##'
+     ##' <details>
+     ##' @title
+     ##' @param ng number of groups
+     ##' @param nr number of "runs", i.e., observations per groups
+     ##' @param sd standard deviations of group and "Individual" random effects,
+     ##'    (\sigma_f, \sigma_I)
+     ##' @param b  true beta (fixed effects)
+     ##' @return a data frame (to be used in glmer()) with columns
+     ##'    (x, f, obs, eta0, eta, mu, y), where y ~ Pois(lambda(x)),
+     ##'                                   log(lambda(x_i)) = b_1 + b_2 * x + G_{f(i)} + I_i
+     ##'    and G_k ~ N(0, \sigma_f);  I_i ~ N(0, \sigma_I)
+     ##' @author Ben Bolker and Martin Maechler
+     rPoisGLMMi <- function(ng, nr, sd=c(f = 1, ind = 0.5), b=c(1,2))
+     {
+         stopifnot(nr >= 1, ng >= 1,
+                   is.numeric(sd), names(sd) %in% c("f","ind"), sd >= 0)
+         ntot <- nr*ng
+         b.reff <- rnorm(ng,  sd= sd[["f"]])
+         b.rind <- rnorm(ntot,sd= sd[["ind"]])
+         x <- runif(ntot)
+         within(data.frame(x,
+                           f = factor(rep(LETTERS[1:ng], each=nr)),
+                           obs = 1:ntot,
+                           eta0 = cbind(1, x) %*% b),
+         {
+             eta <- eta0 + b.reff[f] + b.rind[obs]
+             mu <- exp(eta)
+             y <- rpois(ntot, lambda=mu)
+         })
+     }
+ 
+     set.seed(1)
+     dd <- rPoisGLMMi(12, 20)
+     m0  <- glmer(y~x + (1|f),           family="poisson", data=dd)
+     m1 <- glmer(y~x + (1|f) + (1|obs), family="poisson", data=dd)
+     stopifnot(isTRUE(chkFixed(m0, true.coef = c(1,2))),
+               isTRUE(chkFixed(m1, true.coef = c(1,2))))
+     (a01 <- anova(m0, m1))
+ 
+     stopifnot(all.equal(a01$Chisq[2], 554.334056, tolerance=1e-5),
+               all.equal(a01$logLik, c(-1073.77193, -796.604902), tolerance=1e-6),
+               a01$ npar == 3:4,
+               na.omit(a01$ Df) == 1)
+ 
+     if(lme4:::testLevel() > 1) {
+         nsim <- 10
+         set.seed(2)
+         system.time(
+             simR <- lapply(1:nsim,  function(i) {
+                 cat(i,"", if(i %% 20 == 0)"\n")
+                 dd <- rPoisGLMMi(10 + rpois(1, lambda=3),
+                                  16 + rpois(1, lambda=5))
+                 m0 <- glmer(y~x + (1|f),           family="poisson", data=dd)
+                 m1 <- glmer(y~x + (1|f) + (1|obs), family="poisson", data=dd)
+                 a01 <- anova(m0, m1)
+                 stopifnot(a01$ npar == 3:4,
+                           na.omit(a01$ Df) == 1)
+                 list(chk0 = chkFixed(m0, true.coef = c(1,2)),
+                      chk1 = chkFixed(m1, true.coef = c(1,2)),
+                      chisq= a01$Chisq[2],
+                      lLik = a01$logLik)
+             }))
+ 
+         ## m0 is the wrong model, so we don't expect much here:
+         table(unlist(lapply(simR, `[[`, "chk0")))
+ 
+ 
+         ## If the fixed effect estimates were unbiased and the standard errors correct,
+         ## and N(0,sigma^2) instead of t_{nu} good enough for the fixed effects,
+         ## the confidence interval should contain the true coef in ~95 out of 100:
+         table(unlist(lapply(simR, `[[`, "chk1")))
+ 
+         ## The tests are all highly significantly in favor of  m1 :
+         summary(chi2s <- sapply(simR, `[[`, "chisq"))
+         ##  Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
+         ## 158.9   439.0   611.4   698.2   864.3  2268.0
+         stopifnot(chi2s > qchisq(0.9999, df = 1))
+     }
+ 
+     showProc.time()
+ })  ## skip if windows and testLevel<1
> stopifnot(suppressPackageStartupMessages(require(lme4)))
> options(show.signif.stars = FALSE)
> source(system.file("test-tools-1.R", package = "Matrix"), keep.source = FALSE)
Loading required package: tools
> chkFixed <- function(fm, true.coef, conf.level = 0.95, sd.factor = qnorm((1 + 
+     conf.level)/2)) {
+     stopifnot(is.matrix(cf <- coefficients(summary(fm))), ncol(cf) >= 2)
+     cc <- cf[, 1]
+     sd <- cf[, 2]
+     if (any(out1 <- true.coef < cc - sd.factor * sd)) 
+         return(sprintf("true coefficient[j], j=%s, is smaller than lower confidence limit", 
+             paste(which(out1), collapse = ", ")))
+     if (any(out2 <- true.coef > cc + sd.factor * sd)) 
+         return(sprintf("true coefficient[j], j=%s, is larger than upper confidence limit", 
+             paste(which(out2), collapse = ", ")))
+     TRUE
+ }
> m1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), 
+     family = binomial, data = cbpp)
==3511189== Invalid read of size 8
==3511189==    at 0x24BDA545: coeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925)
==3511189==    by 0x24BDA545: assignCoeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:654)
==3511189==    by 0x24BDA545: assignCoeffByOuterInner (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:668)
==3511189==    by 0x24BDA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:347)
==3511189==    by 0x24BDA545: call_dense_assignment_loop<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:785)
==3511189==    by 0x24BDA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:954)
==3511189==    by 0x24BDA545: call_assignment_no_alias<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDA545: call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:858)
==3511189==    by 0x24BDA545: call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:836)
==3511189==    by 0x24BDA545: operator=<Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Assign.h:66)
==3511189==    by 0x24BDA545: lme4::merPredD::updateXwts(Eigen::Array<double, -1, 1, 0, -1, 1> const&) (packages/tests-vg/lme4/src/predModule.cpp:221)
==3511189==    by 0x24BB2A2F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:274)
==3511189==    by 0x24BB2A2F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Address 0x1c57cd20 is 48 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4D78D5: promiseArgsEx (svn/R-devel/src/main/eval.c:3789)
==3511189==    by 0x4F4712: Rf_promiseArgs (svn/R-devel/src/main/eval.c:3854)
==3511189==    by 0x4F4712: Rf_eval (svn/R-devel/src/main/eval.c:1276)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F760A: R_execMethod (svn/R-devel/src/main/eval.c:2567)
==3511189==    by 0x48BAA10: R_dispatchGeneric (svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD31F4: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD31F4: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD31F4: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD31F4: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3511189==    by 0x24BB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==  Address 0x1c57cd20 is 48 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4D78D5: promiseArgsEx (svn/R-devel/src/main/eval.c:3789)
==3511189==    by 0x4F4712: Rf_promiseArgs (svn/R-devel/src/main/eval.c:3854)
==3511189==    by 0x4F4712: Rf_eval (svn/R-devel/src/main/eval.c:1276)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F760A: R_execMethod (svn/R-devel/src/main/eval.c:2567)
==3511189==    by 0x48BAA10: R_dispatchGeneric (svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD3204: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD3204: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD3204: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD3204: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3511189==    by 0x24BB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==  Address 0x1c57cd30 is 64 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4D78D5: promiseArgsEx (svn/R-devel/src/main/eval.c:3789)
==3511189==    by 0x4F4712: Rf_promiseArgs (svn/R-devel/src/main/eval.c:3854)
==3511189==    by 0x4F4712: Rf_eval (svn/R-devel/src/main/eval.c:1276)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F760A: R_execMethod (svn/R-devel/src/main/eval.c:2567)
==3511189==    by 0x48BAA10: R_dispatchGeneric (svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD3211: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD3211: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD3211: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD3211: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3511189==    by 0x24BB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==  Address 0x1c57cd40 is 80 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4D78D5: promiseArgsEx (svn/R-devel/src/main/eval.c:3789)
==3511189==    by 0x4F4712: Rf_promiseArgs (svn/R-devel/src/main/eval.c:3854)
==3511189==    by 0x4F4712: Rf_eval (svn/R-devel/src/main/eval.c:1276)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F760A: R_execMethod (svn/R-devel/src/main/eval.c:2567)
==3511189==    by 0x48BAA10: R_dispatchGeneric (svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD321E: pmul<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237)
==3511189==    by 0x24BD321E: pmadd<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959)
==3511189==    by 0x24BD321E: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95)
==3511189==    by 0x24BD321E: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3511189==    by 0x24BB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==  Address 0x1c57cd50 is 96 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4D78D5: promiseArgsEx (svn/R-devel/src/main/eval.c:3789)
==3511189==    by 0x4F4712: Rf_promiseArgs (svn/R-devel/src/main/eval.c:3854)
==3511189==    by 0x4F4712: Rf_eval (svn/R-devel/src/main/eval.c:1276)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F760A: R_execMethod (svn/R-devel/src/main/eval.c:2567)
==3511189==    by 0x48BAA10: R_dispatchGeneric (svn/R-devel/src/library/methods/src/methods_list_dispatch.c:1154)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
> m1. <- update(m1, start = getME(m1, c("theta", "fixef")))
> dm1 <- drop1(m1)
> stopifnot(all.equal(drop1(m1.), dm1, tol = 1e-10))
> m1p <- glmer(incidence/size ~ period + (1 | herd), weights = size, family = binomial, 
+     data = cbpp)
> stopifnot(all.equal(fixef(m1), fixef(m1p)), all.equal(ranef(m1), ranef(m1p)), 
+     TRUE)
> stopifnot(all.equal(logLik(m1), logLik(m1p)), all.equal(AIC(m1), AIC(m1p)), 
+     all.equal(BIC(m1), BIC(m1p)))
> m1b <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), 
+     family = binomial, data = cbpp, verbose = 2L, control = glmerControl(optimizer = "bobyqa", 
+         tolPwrss = 1e-07, optCtrl = list(rhobeg = 0.2, rhoend = 2e-07)))
npt = 3 , n =  1 
rhobeg =  0.2 , rhoend =  2e-07 
start par. =  1 fn =  186.7231 
rho:    0.020 eval:   4 fn:      184.166 par:0.600000 
rho:   0.0020 eval:   7 fn:      184.110 par:0.649419 
rho:  0.00020 eval:  10 fn:      184.109 par:0.641956 
rho:  2.0e-05 eval:  12 fn:      184.109 par:0.641847 
rho:  2.0e-06 eval:  13 fn:      184.109 par:0.641847 
rho:  2.0e-07 eval:  15 fn:      184.109 par:0.641839 
At return
eval:  18 fn:      184.10869 par: 0.641839
npt = 7 , n =  5 
rhobeg =  0.2 , rhoend =  2e-07 
start par. =  0.6418386 -1.360476 -0.9761732 -1.111073 -1.559676 fn =  184.1086 
rho:    0.020 eval:   8 fn:      184.109 par:0.641839 -1.36048 -0.976173 -1.11107 -1.55968 
rho:   0.0020 eval:  15 fn:      184.056 par:0.641943 -1.40262 -0.981786 -1.13822 -1.57895 
rho:  0.00020 eval:  30 fn:      184.053 par:0.642120 -1.39845 -0.991559 -1.12809 -1.58002 
rho:  2.0e-05 eval:  37 fn:      184.053 par:0.642074 -1.39838 -0.991818 -1.12817 -1.57968 
rho:  2.0e-06 eval:  49 fn:      184.053 par:0.642064 -1.39834 -0.991914 -1.12821 -1.57974 
rho:  2.0e-07 eval:  57 fn:      184.053 par:0.642064 -1.39833 -0.991924 -1.12821 -1.57975 
At return
eval:  69 fn:      184.05313 par: 0.642064 -1.39833 -0.991924 -1.12821 -1.57975
> m.9 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), 
+     family = binomial, data = cbpp, nAGQ = 9)
==3511189== Invalid read of size 8
==3511189==    at 0x24BDA545: coeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925)
==3511189==    by 0x24BDA545: assignCoeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:654)
==3511189==    by 0x24BDA545: assignCoeffByOuterInner (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:668)
==3511189==    by 0x24BDA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:347)
==3511189==    by 0x24BDA545: call_dense_assignment_loop<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:785)
==3511189==    by 0x24BDA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:954)
==3511189==    by 0x24BDA545: call_assignment_no_alias<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDA545: call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:858)
==3511189==    by 0x24BDA545: call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:836)
==3511189==    by 0x24BDA545: operator=<Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Assign.h:66)
==3511189==    by 0x24BDA545: lme4::merPredD::updateXwts(Eigen::Array<double, -1, 1, 0, -1, 1> const&) (packages/tests-vg/lme4/src/predModule.cpp:221)
==3511189==    by 0x24BB2A2F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:274)
==3511189==    by 0x24BB2A2F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB937D: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:428)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Address 0x1a7b9a10 is 48 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD31F4: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD31F4: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD31F4: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD31F4: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB9674: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:451)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a10 is 48 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD3204: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD3204: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD3204: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD3204: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB9674: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:451)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a20 is 64 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD3211: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD3211: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD3211: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD3211: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB9674: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:451)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a30 is 80 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD321E: pmul<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237)
==3511189==    by 0x24BD321E: pmadd<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959)
==3511189==    by 0x24BD321E: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95)
==3511189==    by 0x24BD321E: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB9674: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:451)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a40 is 96 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD31F4: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD31F4: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD31F4: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD31F4: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB987D: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:457)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a10 is 48 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD3204: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD3204: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD3204: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD3204: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB987D: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:457)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a20 is 64 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD3211: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD3211: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD3211: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD3211: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB987D: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:457)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a30 is 80 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD321E: pmul<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237)
==3511189==    by 0x24BD321E: pmadd<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959)
==3511189==    by 0x24BD321E: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95)
==3511189==    by 0x24BD321E: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB987D: glmerAGQ (packages/tests-vg/lme4/src/external.cpp:457)
==3511189==    by 0x4A91A7: R_doDotCall (svn/R-devel/src/main/dotcode.c:780)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3511189==  Address 0x1a7b9a40 is 96 bytes inside a block of size 1,840 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488)
==3511189==    by 0x4C6E51: Rf_defineVar (svn/R-devel/src/main/envir.c:1996)
==3511189==    by 0x4EF2D0: bcEval_loop (svn/R-devel/src/main/eval.c:7914)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
> m2 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), 
+     family = binomial, data = cbpp, nAGQ = 25)
> stopifnot(is((cm2 <- coef(m2)), "coef.mer"), dim(cm2$herd) == c(15, 4), 
+     all.equal(fixef(m2), c(-1.39922533406847, -0.991407294757321, -1.12782184600404, 
+         -1.57946627431248), tolerance = 5e-04, check.attributes = FALSE), all.equal(c(-2 * 
+         logLik(m2)), 100.010030538022, tolerance = 1e-09), all.equal(deviance(m2), 
+         73.373, tolerance = 1e-05))
> coef_m1_lme4.0 <- structure(c(-1.39853505102576, -0.992334712470269, -1.12867541092127, 
+     -1.58037389566025), .Names = c("(Intercept)", "period2", "period3", "period4"))
> coef_m1_glmmadmb <- structure(c(-1.39853810064827, -0.99233330126975, 
+     -1.12867317840779, -1.58031150854503), .Names = c("(Intercept)", "period2", "period3", 
+     "period4"))
> coef_m1_glmmML <- structure(c(-1.39853234657711, -0.992336901732793, -1.12867036466201, 
+     -1.58030977686564), .Names = c("(Intercept)", "period2", "period3", "period4"))
> stopifnot(is((cm1 <- coef(m1b)), "coef.mer"), dim(cm1$herd) == c(15, 4), 
+     all.equal(fixef(m1b), fixef(m1), tolerance = 4e-05), is.all.equal4(fixef(m1b), 
+         coef_m1_glmmadmb, coef_m1_lme4.0, coef_m1_glmmML, tol = 5e-04))
> showProc.time()
Time (user system elapsed): 57.521 0.247 58.152 
> if (require("MASS", quietly = TRUE)) {
+     bacteria$wk2 <- bacteria$week > 2
+     contrasts(bacteria$trt) <- structure(contr.sdif(3), dimnames = list(NULL, c("diag", 
+         "encourage")))
+     print(fm5 <- glmer(y ~ trt + wk2 + (1 | ID), data = bacteria, family = binomial))
+     showProc.time()
+     stopifnot(all.equal(logLik(fm5), structure(-96.13069, nobs = 220L, nall = 220L, 
+         df = 5L, REML = FALSE, class = "logLik"), tolerance = 5e-04, check.attributes = FALSE), 
+         all.equal(fixef(fm5), c(`(Intercept)` = 2.83160949, trtdiag = -1.366722631, 
+             trtencourage = 0.5840147802, wk2TRUE = -1.598591346), tolerance = 1e-04))
+ }
==3511189== Invalid read of size 16
==3511189==    at 0x24BD360C: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3511189==    by 0x24BD360C: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3511189==    by 0x24BD360C: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3511189==    by 0x24BD360C: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:201)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3511189==    by 0x24BB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==  Address 0x1b193300 is 1,776 bytes inside a block of size 7,088 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3511189==    by 0x52D9A2: Rf_matchArgs_NR (svn/R-devel/src/main/match.c:196)
==3511189==    by 0x52E1B8: Rf_matchArgs_RC (svn/R-devel/src/main/match.c:392)
==3511189==    by 0x4D7080: make_applyClosure_env (svn/R-devel/src/main/eval.c:2261)
==3511189==    by 0x4F7251: applyClosure_core (svn/R-devel/src/main/eval.c:2304)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
==3511189== Invalid read of size 16
==3511189==    at 0x24BD361C: pmul<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237)
==3511189==    by 0x24BD361C: pmadd<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959)
==3511189==    by 0x24BD361C: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95)
==3511189==    by 0x24BD361C: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:201)
==3511189==    by 0x24BE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3511189==    by 0x24BE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3511189==    by 0x24BE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3511189==    by 0x24BDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3511189==    by 0x24BDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3511189==    by 0x24BDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3511189==    by 0x24BDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3511189==    by 0x24BDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3511189==    by 0x24BDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3511189==    by 0x24BDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3511189==    by 0x24BDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3511189==    by 0x24BB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3511189==    by 0x24BB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3511189==    by 0x24BB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3511189==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3511189==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==  Address 0x1b193310 is 1,792 bytes inside a block of size 7,088 free'd
==3511189==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3511189==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3511189==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3511189==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3511189==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3511189==    by 0x52D9A2: Rf_matchArgs_NR (svn/R-devel/src/main/match.c:196)
==3511189==    by 0x52E1B8: Rf_matchArgs_RC (svn/R-devel/src/main/match.c:392)
==3511189==    by 0x4D7080: make_applyClosure_env (svn/R-devel/src/main/eval.c:2261)
==3511189==    by 0x4F7251: applyClosure_core (svn/R-devel/src/main/eval.c:2304)
==3511189==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3511189==    by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3511189==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3511189==  Block was alloc'd at
==3511189==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3511189==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3511189==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3511189==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3511189==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3511189==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3511189==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3511189==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3511189==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3511189== 
Generalized linear mixed model fit by maximum likelihood (Laplace
  Approximation) [glmerMod]
 Family: binomial  ( logit )
Formula: y ~ trt + wk2 + (1 | ID)
   Data: bacteria
      AIC       BIC    logLik -2*log(L)  df.resid 
 202.2614  219.2296  -96.1307  192.2614       215 
Random effects:
 Groups Name        Std.Dev.
 ID     (Intercept) 1.242   
Number of obs: 220, groups:  ID, 50
Fixed Effects:
 (Intercept)       trtdiag  trtencourage       wk2TRUE  
       2.832        -1.367         0.584        -1.598  
Time (user system elapsed): 12.701 0.158 12.941 
> if (.Platform$OS.type == "unix") {
+     Sys.setlocale("LC_MESSAGES", "C")
+     tc <- tryCatch(m2 <- glmer(incidence/size ~ period, weights = size, family = binomial, 
+         data = cbpp), error = function(.) .)
+     stopifnot(inherits(tc, "error"), identical(tc$message, "No random effects terms specified in formula"))
+ }
> rPoisGLMMi <- function(ng, nr, sd = c(f = 1, ind = 0.5), b = c(1, 2)) {
+     stopifnot(nr >= 1, ng >= 1, is.numeric(sd), names(sd) %in% c("f", "ind"), sd >= 
+         0)
+     ntot <- nr * ng
+     b.reff <- rnorm(ng, sd = sd[["f"]])
+     b.rind <- rnorm(ntot, sd = sd[["ind"]])
+     x <- runif(ntot)
+     within(data.frame(x, f = factor(rep(LETTERS[1:ng], each = nr)), obs = 1:ntot, 
+         eta0 = cbind(1, x) %*% b), {
+         eta <- eta0 + b.reff[f] + b.rind[obs]
+         mu <- exp(eta)
+         y <- rpois(ntot, lambda = mu)
+     })
+ }
> set.seed(1)
> dd <- rPoisGLMMi(12, 20)
> m0 <- glmer(y ~ x + (1 | f), family = "poisson", data = dd)
> m1 <- glmer(y ~ x + (1 | f) + (1 | obs), family = "poisson", data = dd)
> stopifnot(isTRUE(chkFixed(m0, true.coef = c(1, 2))), isTRUE(chkFixed(m1, 
+     true.coef = c(1, 2))))
> (a01 <- anova(m0, m1))
Data: dd
Models:
m0: y ~ x + (1 | f)
m1: y ~ x + (1 | f) + (1 | obs)
   npar    AIC    BIC   logLik -2*log(L)  Chisq Df Pr(>Chisq)
m0    3 2153.5 2164.0 -1073.77    2147.5                     
m1    4 1601.2 1615.1  -796.61    1593.2 554.33  1  < 2.2e-16
> stopifnot(all.equal(a01$Chisq[2], 554.334056, tolerance = 1e-05), all.equal(a01$logLik, 
+     c(-1073.77193, -796.604902), tolerance = 1e-06), a01$npar == 3:4, na.omit(a01$Df) == 
+     1)
> if (lme4:::testLevel() > 1) {
+     nsim <- 10
+     set.seed(2)
+     system.time(simR <- lapply(1:nsim, function(i) {
+         cat(i, "", if (i%%20 == 0) 
+             "\n")
+         dd <- rPoisGLMMi(10 + rpois(1, lambda = 3), 16 + rpois(1, lambda = 5))
+         m0 <- glmer(y ~ x + (1 | f), family = "poisson", data = dd)
+         m1 <- glmer(y ~ x + (1 | f) + (1 | obs), family = "poisson", data = dd)
+         a01 <- anova(m0, m1)
+         stopifnot(a01$npar == 3:4, na.omit(a01$Df) == 1)
+         list(chk0 = chkFixed(m0, true.coef = c(1, 2)), chk1 = chkFixed(m1, true.coef = c(1, 
+             2)), chisq = a01$Chisq[2], lLik = a01$logLik)
+     }))
+     table(unlist(lapply(simR, `[[`, "chk0")))
+     table(unlist(lapply(simR, `[[`, "chk1")))
+     summary(chi2s <- sapply(simR, `[[`, "chisq"))
+     stopifnot(chi2s > qchisq(0.9999, df = 1))
+ }
> showProc.time()
Time (user system elapsed): 16.427 0.1 16.63 
> 
> proc.time()
   user  system elapsed 
200.298   1.733 200.957 
==3511189== 
==3511189== HEAP SUMMARY:
==3511189==     in use at exit: 213,406,806 bytes in 45,521 blocks
==3511189==   total heap usage: 800,456 allocs, 754,935 frees, 747,094,536 bytes allocated
==3511189== 
==3511189== 416 bytes in 1 blocks are possibly lost in loss record 200 of 2,809
==3511189==    at 0x484B133: calloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:1675)
==3511189==    by 0x4011F63: UnknownInlinedFun (/usr/src/debug/glibc-2.39-38.fc40.x86_64/elf/../include/rtld-malloc.h:44)
==3511189==    by 0x4011F63: allocate_dtv (/usr/src/debug/glibc-2.39-38.fc40.x86_64/elf/../elf/dl-tls.c:395)
==3511189==    by 0x4012A61: _dl_allocate_tls (/usr/src/debug/glibc-2.39-38.fc40.x86_64/elf/../elf/dl-tls.c:673)
==3511189==    by 0x557CC03: allocate_stack (/usr/src/debug/glibc-2.39-38.fc40.x86_64/nptl/allocatestack.c:431)
==3511189==    by 0x557CC03: pthread_create@@GLIBC_2.34 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/nptl/pthread_create.c:660)
==3511189==    by 0x20968DB6: cli__start_thread (/tmp/RtmpOgXr6R/R.INSTALL3b00f1391d945e/cli/src/thread.c:54)
==3511189==    by 0x20968DF7: clic_start_thread (/tmp/RtmpOgXr6R/R.INSTALL3b00f1391d945e/cli/src/thread.c:78)
==3511189==    by 0x4A924D: R_doDotCall (svn/R-devel/src/main/dotcode.c:760)
==3511189==    by 0x4A97A3: do_dotcall (svn/R-devel/src/main/dotcode.c:1437)
==3511189==    by 0x4EEA31: bcEval_loop (svn/R-devel/src/main/eval.c:8085)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3511189==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3511189==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3511189==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3511189== 
==3511189== LEAK SUMMARY:
==3511189==    definitely lost: 0 bytes in 0 blocks
==3511189==    indirectly lost: 0 bytes in 0 blocks
==3511189==      possibly lost: 416 bytes in 1 blocks
==3511189==    still reachable: 213,406,390 bytes in 45,520 blocks
==3511189==         suppressed: 0 bytes in 0 blocks
==3511189== Reachable blocks (those to which a pointer was found) are not shown.
==3511189== To see them, rerun with: --leak-check=full --show-leak-kinds=all
==3511189== 
==3511189== For lists of detected and suppressed errors, rerun with: -s
==3511189== ERROR SUMMARY: 1597031 errors from 17 contexts (suppressed: 0 from 0)