Skip to content

Commit bed4901

Browse files
committed
synch with online class
1 parent 280d5d7 commit bed4901

File tree

7 files changed

+121
-34
lines changed

7 files changed

+121
-34
lines changed

_freeze/slides/data-wrangling-1/execute-results/html.json

Lines changed: 2 additions & 2 deletions
Large diffs are not rendered by default.

_freeze/slides/working-with-models/execute-results/html.json

Lines changed: 2 additions & 2 deletions
Large diffs are not rendered by default.

_site/search.json

Lines changed: 19 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -665,11 +665,25 @@
665665
"text": "Putting everything into one pipeline\n\n\nlm_summary <- penguins |>\n nest(data = -species) |>\n mutate(\n fit = map(data, ~lm(flipper_length_mm ~ body_mass_g, data = .x)),\n glance_out = map(fit, glance)\n ) |>\n select(species, glance_out) |>\n unnest(cols = glance_out)\n\nlm_summary\n\n# A tibble: 3 × 13\n species r.squared adj.r.squared sigma statistic p.value df logLik AIC\n <fct> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>\n1 Adelie 0.219 0.214 5.80 41.8 1.34e- 9 1 -479. 963.\n2 Gentoo 0.494 0.490 4.63 118. 1.33e-19 1 -362. 730.\n3 Chinstrap 0.412 0.403 5.51 46.2 3.75e- 9 1 -212. 429.\n# ℹ 4 more variables: BIC <dbl>, deviance <dbl>, df.residual <int>, nobs <int>"
666666
},
667667
{
668-
"objectID": "slides/working-with-models.html#make-label-data",
669-
"href": "slides/working-with-models.html#make-label-data",
668+
"objectID": "slides/working-with-models.html#turn-model-summary-into-label-data",
669+
"href": "slides/working-with-models.html#turn-model-summary-into-label-data",
670670
"title": "Working with models",
671-
"section": "Make label data",
672-
"text": "Make label data\n\n\nlabel_data <- lm_summary |>\n mutate(\n rsqr = signif(r.squared, 2), # round to 2 significant digits\n pval = signif(p.value, 2),\n label = glue(\"R^2 = {rsqr}, P = {pval}\"),\n body_mass_g = 6400, flipper_length_mm = 175 # label position in plot\n ) |>\n select(species, label, body_mass_g, flipper_length_mm)\n\nlabel_data\n\n# A tibble: 3 × 4\n species label body_mass_g flipper_length_mm\n <fct> <glue> <dbl> <dbl>\n1 Adelie R^2 = 0.22, P = 1.3e-09 6400 175\n2 Gentoo R^2 = 0.49, P = 1.3e-19 6400 175\n3 Chinstrap R^2 = 0.41, P = 3.7e-09 6400 175"
671+
"section": "Turn model summary into label data",
672+
"text": "Turn model summary into label data\n\n\n# Reminder: summary of the model fits\nlm_summary\n\n# A tibble: 3 × 13\n species r.squared adj.r.squared sigma statistic p.value df logLik AIC\n <fct> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>\n1 Adelie 0.219 0.214 5.80 41.8 1.34e- 9 1 -479. 963.\n2 Gentoo 0.494 0.490 4.63 118. 1.33e-19 1 -362. 730.\n3 Chinstrap 0.412 0.403 5.51 46.2 3.75e- 9 1 -212. 429.\n# ℹ 4 more variables: BIC <dbl>, deviance <dbl>, df.residual <int>, nobs <int>"
673+
},
674+
{
675+
"objectID": "slides/working-with-models.html#turn-model-summary-into-label-data-1",
676+
"href": "slides/working-with-models.html#turn-model-summary-into-label-data-1",
677+
"title": "Working with models",
678+
"section": "Turn model summary into label data",
679+
"text": "Turn model summary into label data\n\n\n# Turn R^2 and p value into labels\nlabel_data <- lm_summary |>\n mutate(\n rsqr = signif(r.squared, 2), # round to 2 significant digits\n pval = signif(p.value, 2),\n label = glue(\"R^2 = {rsqr}, P = {pval}\"),\n ) |>\n select(species, label)\n\nlabel_data\n\n# A tibble: 3 × 2\n species label \n <fct> <glue> \n1 Adelie R^2 = 0.22, P = 1.3e-09\n2 Gentoo R^2 = 0.49, P = 1.3e-19\n3 Chinstrap R^2 = 0.41, P = 3.7e-09"
680+
},
681+
{
682+
"objectID": "slides/working-with-models.html#turn-model-summary-into-label-data-2",
683+
"href": "slides/working-with-models.html#turn-model-summary-into-label-data-2",
684+
"title": "Working with models",
685+
"section": "Turn model summary into label data",
686+
"text": "Turn model summary into label data\n\n\n# Add label position in the plot\nlabel_data <- lm_summary |>\n mutate(\n rsqr = signif(r.squared, 2), # round to 2 significant digits\n pval = signif(p.value, 2),\n label = glue(\"R^2 = {rsqr}, P = {pval}\"),\n body_mass_g = 6400, flipper_length_mm = 175 # label position in plot\n ) |>\n select(species, label, body_mass_g, flipper_length_mm)\n\nlabel_data\n\n# A tibble: 3 × 4\n species label body_mass_g flipper_length_mm\n <fct> <glue> <dbl> <dbl>\n1 Adelie R^2 = 0.22, P = 1.3e-09 6400 175\n2 Gentoo R^2 = 0.49, P = 1.3e-19 6400 175\n3 Chinstrap R^2 = 0.41, P = 3.7e-09 6400 175"
673687
},
674688
{
675689
"objectID": "slides/working-with-models.html#and-plot",
@@ -4540,7 +4554,7 @@
45404554
"href": "slides/data-wrangling-1.html#which-to-use-native-pipe-or-old-school-pipe",
45414555
"title": "Data wrangling 1",
45424556
"section": "Which to use? Native pipe or old-school pipe?",
4543-
"text": "Which to use? Native pipe or old-school pipe?\n\n\nThe native pipe |> is the future. I’ll be using it throughout.\nThe magrittr %>% works fine. You can use it if you llike to."
4557+
"text": "Which to use? Native pipe or old-school pipe?\n\n\nThe native pipe |> is the future. I’ll be using it throughout.\nThe magrittr %>% works fine. You can use it if you like to."
45444558
},
45454559
{
45464560
"objectID": "slides/data-wrangling-1.html#pick-rows-from-a-table-filter",

_site/slides/data-wrangling-1.html

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -10,7 +10,7 @@
1010
<meta name="generator" content="quarto-1.6.39">
1111

1212
<meta name="author" content="Claus O. Wilke">
13-
<meta name="dcterms.date" content="2025-02-23">
13+
<meta name="dcterms.date" content="2025-06-18">
1414
<title>SDS 366 – Data wrangling 1</title>
1515
<meta name="apple-mobile-web-app-capable" content="yes">
1616
<meta name="apple-mobile-web-app-status-bar-style" content="black-translucent">
@@ -205,7 +205,7 @@ <h1 class="title">Data wrangling 1</h1>
205205
</div>
206206
</div>
207207

208-
<p class="date">2025-02-23</p>
208+
<p class="date">2025-06-18</p>
209209
</section>
210210
<section id="elementary-data-manipulations" class="slide level2">
211211
<h2>Elementary data manipulations</h2>
@@ -333,7 +333,7 @@ <h2>Which to use? Native pipe or old-school pipe?</h2>
333333
<div>
334334
<ul>
335335
<li class="fragment">The native pipe <code>|&gt;</code> is the future. I’ll be using it throughout.</li>
336-
<li class="fragment">The magrittr <code>%&gt;%</code> works fine. You can use it if you llike to.</li>
336+
<li class="fragment">The magrittr <code>%&gt;%</code> works fine. You can use it if you like to.</li>
337337
</ul>
338338
</div>
339339
<!-- Segment ends here -->

_site/slides/working-with-models.html

Lines changed: 64 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -10,7 +10,7 @@
1010
<meta name="generator" content="quarto-1.6.39">
1111

1212
<meta name="author" content="Claus O. Wilke">
13-
<meta name="dcterms.date" content="2025-03-06">
13+
<meta name="dcterms.date" content="2025-06-22">
1414
<title>SDS 366 – Working with models</title>
1515
<meta name="apple-mobile-web-app-capable" content="yes">
1616
<meta name="apple-mobile-web-app-status-bar-style" content="black-translucent">
@@ -205,7 +205,7 @@ <h1 class="title">Working with models</h1>
205205
</div>
206206
</div>
207207

208-
<p class="date">2025-03-06</p>
208+
<p class="date">2025-06-22</p>
209209
</section>
210210
<section id="how-do-we-obtain-information-about-model-fits" class="slide level2">
211211
<h2>How do we obtain information about model fits?</h2>
@@ -578,20 +578,64 @@ <h2>Putting everything into one pipeline</h2>
578578
<h1>How do we use this approach to annotate plots?</h1>
579579

580580
</section>
581-
<section id="make-label-data" class="slide level2">
582-
<h2>Make label data</h2>
581+
<section id="turn-model-summary-into-label-data" class="slide level2">
582+
<h2>Turn model summary into label data</h2>
583583
<div class="small-font">
584584
<div class="cell">
585-
<div class="sourceCode cell-code" id="cb32"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb32-1"><a></a>label_data <span class="ot">&lt;-</span> lm_summary <span class="sc">|&gt;</span></span>
586-
<span id="cb32-2"><a></a> <span class="fu">mutate</span>(</span>
587-
<span id="cb32-3"><a></a> <span class="at">rsqr =</span> <span class="fu">signif</span>(r.squared, <span class="dv">2</span>), <span class="co"># round to 2 significant digits</span></span>
588-
<span id="cb32-4"><a></a> <span class="at">pval =</span> <span class="fu">signif</span>(p.value, <span class="dv">2</span>),</span>
589-
<span id="cb32-5"><a></a> <span class="at">label =</span> <span class="fu">glue</span>(<span class="st">"R^2 = {rsqr}, P = {pval}"</span>),</span>
590-
<span id="cb32-6"><a></a> <span class="at">body_mass_g =</span> <span class="dv">6400</span>, <span class="at">flipper_length_mm =</span> <span class="dv">175</span> <span class="co"># label position in plot</span></span>
591-
<span id="cb32-7"><a></a> ) <span class="sc">|&gt;</span></span>
592-
<span id="cb32-8"><a></a> <span class="fu">select</span>(species, label, body_mass_g, flipper_length_mm)</span>
593-
<span id="cb32-9"><a></a></span>
594-
<span id="cb32-10"><a></a>label_data</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
585+
<div class="sourceCode cell-code" id="cb32"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb32-1"><a></a><span class="co"># Reminder: summary of the model fits</span></span>
586+
<span id="cb32-2"><a></a>lm_summary</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
587+
<div class="cell-output cell-output-stdout">
588+
<pre><code># A tibble: 3 × 13
589+
species r.squared adj.r.squared sigma statistic p.value df logLik AIC
590+
&lt;fct&gt; &lt;dbl&gt; &lt;dbl&gt; &lt;dbl&gt; &lt;dbl&gt; &lt;dbl&gt; &lt;dbl&gt; &lt;dbl&gt; &lt;dbl&gt;
591+
1 Adelie 0.219 0.214 5.80 41.8 1.34e- 9 1 -479. 963.
592+
2 Gentoo 0.494 0.490 4.63 118. 1.33e-19 1 -362. 730.
593+
3 Chinstrap 0.412 0.403 5.51 46.2 3.75e- 9 1 -212. 429.
594+
# ℹ 4 more variables: BIC &lt;dbl&gt;, deviance &lt;dbl&gt;, df.residual &lt;int&gt;, nobs &lt;int&gt;</code></pre>
595+
</div>
596+
</div>
597+
</div>
598+
</section>
599+
<section id="turn-model-summary-into-label-data-1" class="slide level2">
600+
<h2>Turn model summary into label data</h2>
601+
<div class="small-font">
602+
<div class="cell">
603+
<div class="sourceCode cell-code" id="cb34"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb34-1"><a></a><span class="co"># Turn R^2 and p value into labels</span></span>
604+
<span id="cb34-2"><a></a>label_data <span class="ot">&lt;-</span> lm_summary <span class="sc">|&gt;</span></span>
605+
<span id="cb34-3"><a></a> <span class="fu">mutate</span>(</span>
606+
<span id="cb34-4"><a></a> <span class="at">rsqr =</span> <span class="fu">signif</span>(r.squared, <span class="dv">2</span>), <span class="co"># round to 2 significant digits</span></span>
607+
<span id="cb34-5"><a></a> <span class="at">pval =</span> <span class="fu">signif</span>(p.value, <span class="dv">2</span>),</span>
608+
<span id="cb34-6"><a></a> <span class="at">label =</span> <span class="fu">glue</span>(<span class="st">"R^2 = {rsqr}, P = {pval}"</span>),</span>
609+
<span id="cb34-7"><a></a> ) <span class="sc">|&gt;</span></span>
610+
<span id="cb34-8"><a></a> <span class="fu">select</span>(species, label)</span>
611+
<span id="cb34-9"><a></a></span>
612+
<span id="cb34-10"><a></a>label_data</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
613+
<div class="cell-output cell-output-stdout">
614+
<pre><code># A tibble: 3 × 2
615+
species label
616+
&lt;fct&gt; &lt;glue&gt;
617+
1 Adelie R^2 = 0.22, P = 1.3e-09
618+
2 Gentoo R^2 = 0.49, P = 1.3e-19
619+
3 Chinstrap R^2 = 0.41, P = 3.7e-09</code></pre>
620+
</div>
621+
</div>
622+
</div>
623+
</section>
624+
<section id="turn-model-summary-into-label-data-2" class="slide level2">
625+
<h2>Turn model summary into label data</h2>
626+
<div class="small-font">
627+
<div class="cell">
628+
<div class="sourceCode cell-code" id="cb36"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb36-1"><a></a><span class="co"># Add label position in the plot</span></span>
629+
<span id="cb36-2"><a></a>label_data <span class="ot">&lt;-</span> lm_summary <span class="sc">|&gt;</span></span>
630+
<span id="cb36-3"><a></a> <span class="fu">mutate</span>(</span>
631+
<span id="cb36-4"><a></a> <span class="at">rsqr =</span> <span class="fu">signif</span>(r.squared, <span class="dv">2</span>), <span class="co"># round to 2 significant digits</span></span>
632+
<span id="cb36-5"><a></a> <span class="at">pval =</span> <span class="fu">signif</span>(p.value, <span class="dv">2</span>),</span>
633+
<span id="cb36-6"><a></a> <span class="at">label =</span> <span class="fu">glue</span>(<span class="st">"R^2 = {rsqr}, P = {pval}"</span>),</span>
634+
<span id="cb36-7"><a></a> <span class="at">body_mass_g =</span> <span class="dv">6400</span>, <span class="at">flipper_length_mm =</span> <span class="dv">175</span> <span class="co"># label position in plot</span></span>
635+
<span id="cb36-8"><a></a> ) <span class="sc">|&gt;</span></span>
636+
<span id="cb36-9"><a></a> <span class="fu">select</span>(species, label, body_mass_g, flipper_length_mm)</span>
637+
<span id="cb36-10"><a></a></span>
638+
<span id="cb36-11"><a></a>label_data</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
595639
<div class="cell-output cell-output-stdout">
596640
<pre><code># A tibble: 3 × 4
597641
species label body_mass_g flipper_length_mm
@@ -607,12 +651,12 @@ <h2>Make label data</h2>
607651
<h2>And plot</h2>
608652
<div class="small-font">
609653
<div class="cell">
610-
<div class="sourceCode cell-code" id="cb34"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb34-1"><a></a><span class="fu">ggplot</span>(penguins, <span class="fu">aes</span>(body_mass_g, flipper_length_mm)) <span class="sc">+</span> <span class="fu">geom_point</span>() <span class="sc">+</span></span>
611-
<span id="cb34-2"><a></a> <span class="fu">geom_text</span>(</span>
612-
<span id="cb34-3"><a></a> <span class="at">data =</span> label_data, <span class="fu">aes</span>(<span class="at">label =</span> label),</span>
613-
<span id="cb34-4"><a></a> <span class="at">size =</span> <span class="dv">10</span>, <span class="at">size.unit =</span> <span class="st">"pt"</span>, <span class="at">hjust =</span> <span class="dv">1</span> <span class="co"># 10pt, right-justified</span></span>
614-
<span id="cb34-5"><a></a> ) <span class="sc">+</span></span>
615-
<span id="cb34-6"><a></a> <span class="fu">geom_smooth</span>(<span class="at">method =</span> <span class="st">"lm"</span>, <span class="at">se =</span> <span class="cn">FALSE</span>) <span class="sc">+</span> <span class="fu">facet_wrap</span>(<span class="fu">vars</span>(species))</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
654+
<div class="sourceCode cell-code" id="cb38"><pre class="sourceCode numberSource r number-lines code-with-copy"><code class="sourceCode r"><span id="cb38-1"><a></a><span class="fu">ggplot</span>(penguins, <span class="fu">aes</span>(body_mass_g, flipper_length_mm)) <span class="sc">+</span> <span class="fu">geom_point</span>() <span class="sc">+</span></span>
655+
<span id="cb38-2"><a></a> <span class="fu">geom_text</span>(</span>
656+
<span id="cb38-3"><a></a> <span class="at">data =</span> label_data, <span class="fu">aes</span>(<span class="at">label =</span> label),</span>
657+
<span id="cb38-4"><a></a> <span class="at">size =</span> <span class="dv">10</span>, <span class="at">size.unit =</span> <span class="st">"pt"</span>, <span class="at">hjust =</span> <span class="dv">1</span> <span class="co"># 10pt, right-justified</span></span>
658+
<span id="cb38-5"><a></a> ) <span class="sc">+</span></span>
659+
<span id="cb38-6"><a></a> <span class="fu">geom_smooth</span>(<span class="at">method =</span> <span class="st">"lm"</span>, <span class="at">se =</span> <span class="cn">FALSE</span>) <span class="sc">+</span> <span class="fu">facet_wrap</span>(<span class="fu">vars</span>(species))</span></code><button title="Copy to Clipboard" class="code-copy-button"><i class="bi"></i></button></pre></div>
616660
<div class="cell-output-display">
617661
<div class="quarto-figure quarto-figure-center">
618662
<figure>

slides/data-wrangling-1.qmd

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -110,7 +110,7 @@ penguins %>%
110110

111111
::: {.incremental}
112112
- The native pipe `|>` is the future. I'll be using it throughout.
113-
- The magrittr `%>%` works fine. You can use it if you llike to.
113+
- The magrittr `%>%` works fine. You can use it if you like to.
114114
:::
115115

116116
<!-- Segment ends here -->

slides/working-with-models.qmd

Lines changed: 30 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -246,11 +246,40 @@ lm_summary
246246

247247
# How do we use this approach to annotate plots?
248248

249-
## Make label data
249+
## Turn model summary into label data
250250

251251
::: {.small-font}
252252
```{r}
253253
#| echo: true
254+
# Reminder: summary of the model fits
255+
lm_summary
256+
```
257+
:::
258+
259+
## Turn model summary into label data
260+
261+
::: {.small-font}
262+
```{r}
263+
#| echo: true
264+
# Turn R^2 and p value into labels
265+
label_data <- lm_summary |>
266+
mutate(
267+
rsqr = signif(r.squared, 2), # round to 2 significant digits
268+
pval = signif(p.value, 2),
269+
label = glue("R^2 = {rsqr}, P = {pval}"),
270+
) |>
271+
select(species, label)
272+
273+
label_data
274+
```
275+
:::
276+
277+
## Turn model summary into label data
278+
279+
::: {.small-font}
280+
```{r}
281+
#| echo: true
282+
# Add label position in the plot
254283
label_data <- lm_summary |>
255284
mutate(
256285
rsqr = signif(r.squared, 2), # round to 2 significant digits

0 commit comments

Comments
 (0)