Memodelkan sebuah interaksi
Dalam latihan ini, Anda akan menggunakan interaksi untuk memodelkan pengaruh jenis kelamin dan aktivitas lambung terhadap metabolisme alkohol.
Data frame alcohol telah dimuat sebelumnya, dan memiliki kolom:
Metabol: laju metabolisme alkoholGastric: laju aktivitas enzim alkohol dehidrogenase di lambungSex: jenis kelamin peminum (MaleatauFemale)
Dalam video, kita memasangkan tiga model pada data alcohol:
- satu dengan hanya term aditif (efek utama):
Metabol ~ Gastric + Sex - dua model, masing-masing dengan interaksi antara aktivitas lambung dan jenis kelamin
Anda melihat bahwa salah satu model dengan term interaksi memiliki R-squared yang lebih baik daripada model aditif, yang menunjukkan bahwa penggunaan term interaksi memberikan kecocokan yang lebih baik. Dalam latihan ini, Anda akan membandingkan R-squared dari salah satu model interaksi dengan model yang hanya berisi efek utama.
Ingat bahwa operator : menyatakan interaksi antara dua variabel. Operator * menyatakan interaksi antara dua variabel tersebut, ditambah efek utamanya.
x*y = x + y + x:y
Latihan ini adalah bagian dari kursus
Supervised Learning di R: Regresi
Petunjuk latihan
- Tulis sebuah formula yang menyatakan
Metabolsebagai fungsi dariGastricdanSextanpa interaksi.- Tetapkan formula tersebut ke variabel
fmla_adddan cetak.
- Tetapkan formula tersebut ke variabel
- Tulis sebuah formula yang menyatakan
Metabolsebagai fungsi dari interaksi antaraGastricdanSex.- Tambahkan
Gastricsebagai efek utama, tetapi bukanSex. - Tetapkan formula tersebut ke variabel
fmla_interactiondan cetak.
- Tambahkan
- Pasangkan model linear hanya dengan efek utama:
model_addpada data. - Pasangkan model linear dengan interaksi:
model_interactionpada data. - Panggil
summary()pada kedua model. Mana yang memiliki R-squared lebih baik?
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# 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
___
___