Cara membuat plot interaksi di r

ANOVA dua arah digunakan untuk menentukan apakah terdapat perbedaan antara rata-rata tiga atau lebih kelompok independen yang dibagi berdasarkan dua faktor.

Kami menggunakan ANOVA dua arah ketika kami ingin mengetahui apakah dua faktor spesifik mempengaruhi variabel respon tertentu.

Namun terkadang terdapat efek interaksi antara kedua faktor tersebut, yang dapat mempengaruhi cara kita menafsirkan hubungan antara faktor dan variabel respon.

Misalnya, kita mungkin ingin mengetahui apakah faktor (1) olahraga dan (2) jenis kelamin memengaruhi variabel respons penurunan berat badan . Meskipun kedua faktor tersebut mungkin mempengaruhi penurunan berat badan, ada kemungkinan juga keduanya berinteraksi satu sama lain.

Misalnya, olahraga mungkin menyebabkan penurunan berat badan pada tingkat yang berbeda pada pria dan wanita. Dalam hal ini terdapat efek interaksi antara olahraga dan gender.

Cara paling sederhana untuk mendeteksi dan memahami pengaruh interaksi antara dua faktor adalah dengan menggunakan grafik interaksi .

Ini adalah jenis plot yang menampilkan nilai pas dari variabel respons pada sumbu y dan nilai faktor pertama pada sumbu x. Sedangkan garis-garis pada grafik mewakili nilai faktor kepentingan kedua.

Tutorial ini menjelaskan cara membuat dan menafsirkan plot interaksi di R.

Contoh: plot interaksi di R

Katakanlah peneliti ingin menentukan apakah intensitas olahraga dan jenis kelamin berdampak pada penurunan berat badan. Untuk menguji hal ini, mereka merekrut 30 pria dan 30 wanita untuk berpartisipasi dalam percobaan di mana mereka secara acak menugaskan 10 dari mereka untuk mengikuti program tanpa olahraga, olahraga ringan, atau olahraga intens selama sebulan.

Gunakan langkah-langkah berikut untuk membuat bingkai data di R, lakukan ANOVA dua arah, dan buat plot interaksi untuk memvisualisasikan efek interaksi antara olahraga dan gender.

Langkah 1: Buat datanya.

Kode berikut menunjukkan cara membuat bingkai data di R:

 #make this example reproducible
set.seed(10)

#create data frame
data <- data.frame(gender = rep (c("Male", "Female"), each = 30 ),
                   exercise = rep (c("None", "Light", "Intense"), each = 10 , times = 2 ),
                   weight_loss = c(runif(10, -3, 3), runif(10, 0, 5), runif(10, 5, 9),
                                   runif(10, -4, 2), runif(10, 0, 3), runif(10, 3, 8)))

#view first six rows of data frame
head(data)

  gender exercise weight_loss
1 Male None 0.04486922
2 Male None -1.15938896
3 Male None -0.43855400
4 Male None 1.15861249
5 Male None -2.48918419
6 Male None -1.64738030

Langkah 2: Pasangkan model ANOVA dua arah.

Kode berikut menunjukkan cara memasukkan ANOVA dua arah ke data:

 #fit the two-way ANOVA model
model <- aov(weight_loss ~ gender * exercise, data = data)

#view the model output
summary(model)

# Df Sum Sq Mean Sq F value Pr(>F)    
#gender 1 15.8 15.80 11.197 0.0015 ** 
#exercise 2 505.6 252.78 179.087 <2e-16 ***
#gender:exercise 2 13.0 6.51 4.615 0.0141 *  
#Residuals 54 76.2 1.41                   
#---
#Significant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Perhatikan bahwa nilai p ( 0,0141 ) untuk interaksi antara olahraga dan gender signifikan secara statistik, yang menunjukkan bahwa terdapat efek interaksi yang signifikan antara kedua faktor tersebut.

Langkah 3: Buat plot interaksi.

Kode berikut menunjukkan cara membuat plot interaksi untuk olahraga dan gender:

 interaction.plot(x.factor = data$exercise, #x-axis variable
                 trace.factor = data$gender, #variable for lines
                 response = data$weight_loss, #y-axis variable
                 fun = median, #metric to plot
                 ylab = "Weight Loss",
                 xlab = "Exercise Intensity",
                 col = c("pink", "blue"),
                 lty = 1, #line type
                 lwd = 2, #linewidth
                 trace.label = "Gender")

Secara umum, jika dua garis plot interaksi sejajar maka tidak ada efek interaksi. Namun, jika garis tersebut bersilangan, kemungkinan besar terjadi efek interaksi.

Dari grafik tersebut terlihat adanya persilangan garis antara laki-laki dan perempuan, yang menunjukkan kemungkinan adanya pengaruh interaksi antara variabel intensitas olahraga dan gender.

Hal ini sesuai dengan fakta bahwa nilai p pada hasil tabel ANOVA signifikan secara statistik untuk istilah interaksi dalam model ANOVA.

Sumber daya tambahan

Bagaimana melakukan ANOVA satu arah di R
Cara melakukan ANOVA dua arah di R