penguins-regr02

bayes

rope

string

Published

May 8, 2023

library(tidyverse)

── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr     1.1.3     ✔ readr     2.1.4
✔ forcats   1.0.0     ✔ stringr   1.5.0
✔ ggplot2   3.4.4     ✔ tibble    3.2.1
✔ lubridate 1.9.3     ✔ tidyr     1.3.0
✔ purrr     1.0.2     
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

Aufgabe

Beantworten Sie folgende Forschungsfrage:

Gibt es einen Zusammenhang von Schnabellänge und Gewicht (AV) bei Pinguinen?

Hinweise:

Nutzen Sie den Datensatz aus dem R-Paket palmerpenguins.
Verwenden Sie das Rope-Verfahren

Lösung

Wir rufen Stan:

library(rstanarm)

Loading required package: Rcpp

This is rstanarm version 2.26.1

- See https://mc-stan.org/rstanarm/articles/priors for changes to default priors!

- Default priors may change, so it's safest to specify priors, even if equivalent to the defaults.

- For execution on a local, multicore CPU with excess RAM we recommend calling

  options(mc.cores = parallel::detectCores())

library(easystats)

# Attaching packages: easystats 0.6.0 (red = needs update)
✔ bayestestR  0.13.1   ✔ correlation 0.8.4 
✔ datawizard  0.9.0    ✔ effectsize  0.8.6 
✔ insight     0.19.6   ✔ modelbased  0.8.6 
✔ performance 0.10.8   ✔ parameters  0.21.3
✔ report      0.5.7    ✖ see         0.8.0 

Restart the R-Session and update packages in red with `easystats::easystats_update()`.

library(tidyverse)
library(palmerpenguins)
data(penguins)
m1 <- stan_glm(body_mass_g ~ bill_length_mm, 
               seed = 42,
               refresh = 0,
               data = penguins)

plot(rope(m1))

Hier ist also keine klare Aussage zur Frage, ob der Effekt vernachlässigbar klein ist oder größer, möglich.

Categories:

lm
bayes
rope
string