R Under development (unstable) (2023-11-02 r85465 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. > library(lme4) Loading required package: Matrix > > ## setup > ## library(polytomous) > ## data(think) > ## think.polytomous.lmer1 <- polytomous(Lexeme ~ Agent + Patient + (1|Register), > ## data=think, heuristic="poisson.reformulation") > ## save("formula.poisson","data.poisson",file="polytomous_test.RData") > > load(system.file("testdata","polytomous_test.RData",package="lme4")) > > if (FALSE) { + ## infinite loop + glmer(formula.poisson,data=data.poisson,family=poisson,verbose=10) + ## Cholmod not positive definite -> infinite loop + glmer(formula.poisson,data=data.poisson,family=poisson, + verbose=10,control=glmerControl(optimizer="bobyqa")) + ## caught warning: maxfun < 10 * length(par)^2 is not recommended. -> infinite loop + } > ## works but sloooow .... > if (FALSE) { + try(g1 <- glmer(formula.poisson,data=data.poisson,family=poisson, + control=glmerControl(compDev=FALSE),verbose=1)) + ## runs for 2880 steps until: + ## Error in pp$updateDecomp() : Downdated VtV is not positive definite + } > > (testLevel <- lme4:::testLevel()) [1] 1 > if (testLevel > 2) { + glmer(formula.poisson,data=data.poisson,family=poisson, + control=glmerControl(compDev=FALSE,optimizer="bobyqa")) + ## caught warning: maxfun < 10 * length(par)^2 is not recommended. + ## but runs to completion + } > > > > > > proc.time() user system elapsed 1.64 0.34 1.92