R Under development (unstable) (2024-06-18 r86781 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 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 
   0.98    0.17    1.14