Aan de slagGa gratis aan de slag

Een interactie modelleren

In deze oefening gebruik je interacties om het effect van geslacht en maagactiviteit op alcoholmetabolisme te modelleren.

De alcohol-dataframe is vooraf geladen en bevat de kolommen:

  • Metabol: de snelheid van alcoholmetabolisme
  • Gastric: de activiteitssnelheid van maag-alcoholdehydrogenase
  • Sex: het geslacht van de drinker (Male of Female)

In de video hebben we drie modellen op de alcohol-gegevens gefit:

  • één met alleen additieve (hoofdeffect) termen: Metabol ~ Gastric + Sex
  • twee modellen, elk met interacties tussen maagactiviteit en geslacht

Je zag dat een van de modellen met interactietermen een hogere R-squared had dan het additieve model, wat suggereert dat interactietermen een betere fit geven. In deze oefening vergelijk je de R-squared van een van de interactiemodellen met het model met alleen hoofdeffecten.

Onthoud dat de operator : de interactie tussen twee variabelen aangeeft. De operator * geeft de interactie tussen de twee variabelen aan, plus de hoofdeffecten.

x*y = x + y + x:y

Deze oefening maakt deel uit van de cursus

Supervised Learning in R: Regressie

Cursus bekijken

Oefeninstructies

  • Schrijf een formule die Metabol uitdrukt als functie van Gastric en Sex zonder interacties.
    • Wijs de formule toe aan de variabele fmla_add en print die.
  • Schrijf een formule die Metabol uitdrukt als functie van de interactie tussen Gastric en Sex.
    • Voeg Gastric toe als hoofdeffect, maar niet Sex.
    • Wijs de formule toe aan de variabele fmla_interaction en print die.
  • Fit een lineair model met alleen hoofdeffecten: model_add op de gegevens.
  • Fit een lineair model met de interactie: model_interaction op de gegevens.
  • Roep summary() aan op beide modellen. Welke heeft een betere R-squared?

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# alcohol is available
summary(alcohol)

# Create the formula with main effects only
(fmla_add <- ___ )

# Create the formula with interactions
(fmla_interaction <- ___ )

# Fit the main effects only model
model_add <- ___

# Fit the interaction model
model_interaction <- ___

# Call summary on both models and compare
___
___
Code bewerken en uitvoeren