Моделювання взаємодії
У цій вправі ви використаєте взаємодії, щоб змоделювати вплив статі та шлункової активності на метаболізм алкоголю.
Датафрейм alcohol уже завантажено. Він містить стовпці:
Metabol: швидкість метаболізму алкоголюGastric: швидкість активності шлункової алкогольдегідрогеназиSex: стать того, хто вживає алкоголь (MaleабоFemale)
У відео ми підігнали три моделі до даних alcohol:
- одну лише з додатковими (основними) ефектами:
Metabol ~ Gastric + Sex - дві моделі з взаємодіями між шлунковою активністю та статтю
Ви бачили, що одна з моделей із термінами взаємодії мала кращий R-squared, ніж адитивна модель. Це свідчить, що додавання взаємодій дає кращу відповідність. У цій вправі ви порівняєте R-squared однієї з моделей із взаємодією з моделлю лише з головними ефектами.
Нагадаємо: оператор : позначає взаємодію між двома змінними. Оператор * позначає взаємодію між двома змінними плюс головні ефекти.
x*y = x + y + x:y
Ця вправа є частиною курсу
Кероване навчання в R: регресія
Інструкції до вправи
- Напишіть формулу, яка виражає
Metabolяк функцію відGastricтаSexбез взаємодій.- Присвойте формулу змінній
fmla_addі виведіть її.
- Присвойте формулу змінній
- Напишіть формулу, яка виражає
Metabolяк функцію взаємодії міжGastricіSex.- Додайте
Gastricяк головний ефект, але неSex. - Присвойте формулу змінній
fmla_interactionі виведіть її.
- Додайте
- Навчіть лінійну модель лише з головними ефектами:
model_addна цих даних. - Навчіть лінійну модель з взаємодією:
model_interactionна цих даних. - Викличте
summary()для обох моделей. У якої R-squared кращий?
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# 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
___
___