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 alcoholmetabolismeGastric: de activiteitssnelheid van maag-alcoholdehydrogenaseSex: het geslacht van de drinker (MaleofFemale)
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
Oefeninstructies
- Schrijf een formule die
Metaboluitdrukt als functie vanGastricenSexzonder interacties.- Wijs de formule toe aan de variabele
fmla_adden print die.
- Wijs de formule toe aan de variabele
- Schrijf een formule die
Metaboluitdrukt als functie van de interactie tussenGastricenSex.- Voeg
Gastrictoe als hoofdeffect, maar nietSex. - Wijs de formule toe aan de variabele
fmla_interactionen print die.
- Voeg
- Fit een lineair model met alleen hoofdeffecten:
model_addop de gegevens. - Fit een lineair model met de interactie:
model_interactionop 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
___
___