R Under development (unstable) (2026-01-24 r89326 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2026 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(MEMSS)
Loading required package: lme4
Loading required package: Matrix

Attaching package: 'MEMSS'

The following objects are masked from 'package:datasets':

    CO2, Orange, Theoph

> options(show.signif.stars = FALSE, digits = 3)
> old_optimizer = lmerControl(optimizer = "bobyqa")  # use an older default
> 
> m1 <- lmer(pixel ~ day + I(day^2) + (1|Dog:Side) + (day|Dog),
+            Pixel, control = old_optimizer)
> print(m1, corr = FALSE)
Linear mixed model fit by REML ['lmerMod']
Formula: pixel ~ day + I(day^2) + (1 | Dog:Side) + (day | Dog)
   Data: Pixel
REML criterion at convergence: 825
Random effects:
 Groups   Name        Std.Dev. Corr 
 Dog:Side (Intercept) 16.82         
 Dog      (Intercept) 28.37         
          day          1.84    -0.55
 Residual              8.99         
Number of obs: 102, groups:  Dog:Side, 20; Dog, 10
Fixed Effects:
(Intercept)          day     I(day^2)  
   1073.339        6.130       -0.367  
> unclass(ranef(m1))
$`Dog:Side`
    (Intercept)
A:L       -5.97
A:R      -11.77
B:L       -4.05
B:R       -9.80
C:L       -7.97
C:R        3.73
D:L        2.11
D:R        2.85
E:L        3.31
E:R        7.11
F:L       10.91
F:R       -9.90
G:L       -5.91
G:R       21.12
H:L        8.13
H:R       -7.42
I:L       11.17
I:R      -17.05
J:L       31.61
J:R      -22.20

$Dog
  (Intercept)     day
A      -24.71 -1.1954
B      -23.58 -0.4324
C      -27.08  2.1948
D      -16.66  3.0960
E       25.30 -0.5613
F       10.82 -1.0370
G       49.35 -2.2745
H       -7.05  0.9903
I       -5.75 -0.6811
J       19.37 -0.0994

> 
> proc.time()
   user  system elapsed 
   1.39    0.12    1.51