Modellare un'interazione
In questo esercizio userai le interazioni per modellare l'effetto del sesso e dell'attività gastrica sul metabolismo dell'alcol.
Il data frame alcohol è già caricato e contiene le colonne:
Metabol: la velocità di metabolismo dell'alcolGastric: il tasso di attività della alcol deidrogenasi gastricaSex: il sesso della persona che beve (MaleoFemale)
Nel video abbiamo adattato tre modelli ai dati alcohol:
- uno con soli termini additivi (effetti principali):
Metabol ~ Gastric + Sex - due modelli, ciascuno con interazioni tra attività gastrica e sesso
Hai visto che uno dei modelli con termini di interazione aveva un R-squared migliore rispetto al modello additivo, suggerendo che usare termini di interazione offre un adattamento migliore. In questo esercizio confronterai l'R-squared di uno dei modelli con interazione con quello del modello con soli effetti principali.
Ricorda che l'operatore : indica l'interazione tra due variabili. L'operatore * indica l'interazione tra le due variabili, più gli effetti principali.
x*y = x + y + x:y
Questo esercizio fa parte del corso
Supervised Learning in R: Regression
Istruzioni dell'esercizio
- Scrivi una formula che esprima
Metabolcome funzione diGastriceSexsenza interazioni.- Assegna la formula alla variabile
fmla_adde stampala.
- Assegna la formula alla variabile
- Scrivi una formula che esprima
Metabolcome funzione dell'interazione traGastriceSex.- Aggiungi
Gastriccome effetto principale, ma nonSex. - Assegna la formula alla variabile
fmla_interactione stampala.
- Aggiungi
- Adatta un modello lineare con soli effetti principali:
model_addai dati. - Adatta un modello lineare con l'interazione:
model_interactionai dati. - Chiama
summary()su entrambi i modelli. Quale ha un R-squared migliore?
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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
___
___