Randomisatie-inferentie en t-inferentie vergelijken

Als aan de technische voorwaarden wordt voldaan (zie het volgende hoofdstuk), zouden de conclusies van de randomisatietoets en de toets met de t-verdeling overeenkomen. Ze geven niet exact hetzelfde antwoord, omdat ze op verschillende methoden zijn gebaseerd. Maar de p-waarden en betrouwbaarheidsintervallen zouden redelijk dicht bij elkaar moeten liggen.

Deze oefening maakt deel uit van de cursus

Inferentie voor lineaire regressie in R

Bekijk cursus

Oefeninstructies

Bereken de absolute waarde van de geobserveerde helling, obs_slope.
Voeg een kolom toe aan perm_slope met de absolute waarde van de helling in elke gerandomiseerde replicatie. De kolom met de helling heet stat.
Bereken in de aanroep van summarize() de p-waarde als de proportie absolute schattingen van de helling die minstens zo extreem zijn als de absolute geobserveerde schatting van de helling.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

# The slope in the observed data and each permutation replicate
obs_slope
perm_slope

# Calculate the absolute value of the observed slope
abs_obs_slope <- ___

# Find the p-value
perm_slope %>%
  # Add a column for the absolute value of stat
  ___ %>%
  summarize(
    # Calculate prop'n permuted values at least as extreme as observed
    p_value = ___
  )

Code bewerken en uitvoeren