From bdaa10785eeca7961a0c29fb57713db40e988081 Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Tue, 26 Aug 2025 15:14:20 +0200
Subject: [PATCH 1/6] tuesday afternoon

---
 sesssions/2025-08-26-tuesday-pm.jl | 8 ++++++++
 1 file changed, 8 insertions(+)
 create mode 100644 sesssions/2025-08-26-tuesday-pm.jl

diff --git a/sesssions/2025-08-26-tuesday-pm.jl b/sesssions/2025-08-26-tuesday-pm.jl
new file mode 100644
index 0000000..712befc
--- /dev/null
+++ b/sesssions/2025-08-26-tuesday-pm.jl
@@ -0,0 +1,8 @@
+using CairoMakie
+using MixedModels
+using MixedModelsMakie
+using SMLP2025: datasets
+
+kb07 = dataset("kb07")
+
+fm_max = lmm(@formula(rt_))

From d51e434fac3e62dff9a44f0714cbdbb6120b9362 Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Tue, 26 Aug 2025 15:19:02 +0200
Subject: [PATCH 2/6] moar

---
 sesssions/2025-08-26-tuesday-pm.jl | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/sesssions/2025-08-26-tuesday-pm.jl b/sesssions/2025-08-26-tuesday-pm.jl
index 712befc..df3558f 100644
--- a/sesssions/2025-08-26-tuesday-pm.jl
+++ b/sesssions/2025-08-26-tuesday-pm.jl
@@ -1,8 +1,14 @@
 using CairoMakie
 using MixedModels
 using MixedModelsMakie
-using SMLP2025: datasets
+using SMLP2025: dataset
 
 kb07 = dataset("kb07")
+contrasts = Dict(:spkr => HelmertCoding(),
+                 :prec => HelmertCoding(),
+                 :load => HelmertCoding())
 
-fm_max = lmm(@formula(rt_))
+fm_max = lmm(@formula(rt_trunc ~ 1 + spkr * load * prec +
+                                (1 + spkr * load * prec | subj) +
+                                (1 + spkr * load * prec | item)),
+             kb07; contrasts=contrasts)

From b5fed7b9d4da9191124959e7a99abf974b7b6df3 Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Tue, 26 Aug 2025 15:28:13 +0200
Subject: [PATCH 3/6] wip

---
 sesssions/2025-08-26-tuesday-pm.jl | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/sesssions/2025-08-26-tuesday-pm.jl b/sesssions/2025-08-26-tuesday-pm.jl
index df3558f..c653045 100644
--- a/sesssions/2025-08-26-tuesday-pm.jl
+++ b/sesssions/2025-08-26-tuesday-pm.jl
@@ -3,6 +3,11 @@ using MixedModels
 using MixedModelsMakie
 using SMLP2025: dataset
 
+fm_sleep = lmm(@formula(reaction ~ 1 + days + (1+days|subj)),
+               dataset("sleepstudy"))
+
+shrinkageplot(fm_sleep)
+
 kb07 = dataset("kb07")
 contrasts = Dict(:spkr => HelmertCoding(),
                  :prec => HelmertCoding(),

From 707055b4c7a70855b4cd33bf1b7315d03199308e Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Tue, 26 Aug 2025 15:55:18 +0200
Subject: [PATCH 4/6] wip

---
 sesssions/2025-08-26-tuesday-pm.jl | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/sesssions/2025-08-26-tuesday-pm.jl b/sesssions/2025-08-26-tuesday-pm.jl
index c653045..ba85496 100644
--- a/sesssions/2025-08-26-tuesday-pm.jl
+++ b/sesssions/2025-08-26-tuesday-pm.jl
@@ -8,6 +8,10 @@ fm_sleep = lmm(@formula(reaction ~ 1 + days + (1+days|subj)),
 
 shrinkageplot(fm_sleep)
 
+shrinkageplot!(Figure(;size=(600,600)), fm_sleep; ellipse=true)
+
+caterpillar(fm_sleep)
+
 kb07 = dataset("kb07")
 contrasts = Dict(:spkr => HelmertCoding(),
                  :prec => HelmertCoding(),
@@ -17,3 +21,9 @@ fm_max = lmm(@formula(rt_trunc ~ 1 + spkr * load * prec +
                                 (1 + spkr * load * prec | subj) +
                                 (1 + spkr * load * prec | item)),
              kb07; contrasts=contrasts)
+
+shrinkageplot(fm_max, :subj)
+
+fm_max.rePCA[:subj]
+
+fm_max.PCA[:subj]

From bf1d377b48d82624bfcbe9524e3061a533781b11 Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Tue, 26 Aug 2025 16:38:18 +0200
Subject: [PATCH 5/6] saving

---
 sesssions/2025-08-26-tuesday-pm.jl | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/sesssions/2025-08-26-tuesday-pm.jl b/sesssions/2025-08-26-tuesday-pm.jl
index ba85496..a245bfd 100644
--- a/sesssions/2025-08-26-tuesday-pm.jl
+++ b/sesssions/2025-08-26-tuesday-pm.jl
@@ -8,8 +8,19 @@ fm_sleep = lmm(@formula(reaction ~ 1 + days + (1+days|subj)),
 
 shrinkageplot(fm_sleep)
 
+# notice the ! -- we're modifying an existing figure
+# (admittedly one that we created in the same line)
+# the advantage to this here is that we can specify the
+# figure's size manually
 shrinkageplot!(Figure(;size=(600,600)), fm_sleep; ellipse=true)
 
+# how do we save a figure?
+f = shrinkageplot!(Figure(;size=(600,600)), fm_sleep; ellipse=true)
+
+save("figure.png", f)
+save("figure.pdf", f)
+save("figure.svg", f)
+
 caterpillar(fm_sleep)
 
 kb07 = dataset("kb07")

From 8d116ca82369da189df1fb32eb4b914991076cf9 Mon Sep 17 00:00:00 2001
From: Phillip Alday <me@phillipalday.com>
Date: Tue, 26 Aug 2025 16:52:39 +0200
Subject: [PATCH 6/6] show me the methods

---
 sesssions/2025-08-26-tuesday-pm.jl | 20 ++++++++++++++++++++
 1 file changed, 20 insertions(+)

diff --git a/sesssions/2025-08-26-tuesday-pm.jl b/sesssions/2025-08-26-tuesday-pm.jl
index a245bfd..ccb7bc5 100644
--- a/sesssions/2025-08-26-tuesday-pm.jl
+++ b/sesssions/2025-08-26-tuesday-pm.jl
@@ -38,3 +38,23 @@ shrinkageplot(fm_max, :subj)
 fm_max.rePCA[:subj]
 
 fm_max.PCA[:subj]
+
+
+# one final thing....
+# Julia makes it easy to show() things in different formats
+
+show(stdout, MIME("text/plain"), fm_max)
+
+show(stdout, MIME("text/markdown"), fm_max)
+
+show(stdout, MIME("text/markdown"), VarCorr(fm_max))
+
+show(stdout, MIME("text/latex"), fm_max)
+
+show(stdout, MIME("text/xelatex"), fm_max)
+
+show(stdout, MIME("text/html"), fm_max)
+
+open("mytable.md", "w") do io
+    show(io,  MIME("text/markdown"), fm_max)
+end