R Under development (unstable) (2023-11-01 r85459 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 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. > ## data set and formula extracted from ?prLogisticDelta example > ## (Thailand, clustered-data) in prLogistic package > load(system.file("testdata","prLogistic.RData",package="lme4")) > library(lme4) Loading required package: Matrix > > (testLevel <- lme4:::testLevel()) [1] 1 > if (testLevel > 2) { + print(system.time( + lme4_est <- glmer(rgi ~ sex + pped + (1|schoolid), + data = dataset, family=binomial) + )) + lme4_results <- list(sigma= sqrt(unname(unlist(VarCorr(lme4_est)))), + beta = fixef(lme4_est)) + + ## stored results from other pkgs + glmmML_est <- list(sigma = 1.25365353546143, + beta = c("(Intercept)" = -2.19478801858317, + "sex" = 0.548884468743364, "pped"= -0.623835613907385)) + lme4.0_est <- list(sigma = 1.25369539060849, + beta = c("(Intercept)" = -2.19474529099587, + "sex" = 0.548900267825802, "pped"= -0.623934772981894)) + + source(system.file("test-tools-1.R", package = "Matrix"))#-> assert.EQ() etc + assert.EQ.(lme4_results, glmmML_est, tol=3e-3) + assert.EQ.(lme4_results, lme4.0_est, tol=3e-3) + print(lme4_est) + } > > proc.time() user system elapsed 1.10 0.12 1.21