Dalam pertemuan ini, anda akan mempelajari cara menggunakan Rstudio untuk melakukan Two-Way ANOVA. Two-Way ANOVA bermanfaat untuk mencari faktor yang paling berpengaruh terhadap hasil eksperimen. Dengan mengetahui faktor tersebut, anda dapat mendesain eksperimen dengan protokol atau prosedur yang lebih baik demi tercapainya hasil yang maksimal. ## Definisi data Sebelum anda melangkah lebih jauh, anda perlu mendefinisikan semua data yang anda punya. Pada slide yang diberikan sehari lalu, dicontohkan pengerjaan Two-Way ANOVA secara manual menggunakan data percobaan ligan pengelat A, B, C, dan D. Logam pengelat digunakan untuk mengekstrak ion logam dalam dari sebuah larutan. Banyaknya konsentrasi ion logam yang diekstrak, diduga dipengaruhi oleh jenis pengelat. Konsentrasi ion logam yang diekstrak didefinisikan di bawah ini:
c_A = c(84,79,83)
c_B = c(80,77,78)
c_C = c(83,80,80)
c_D = c(79,79,78)Setelah berhasil mendefinisikan konsentrasi ion logam, semua konsentrasi ion logam digabungkan menjadi satu sebagai variabel konsentrasi. Karena percobaan dihasilkan pada hari yang berbeda untuk masing-masing ion logam, maka hari atau dengan kata lain, waktu percobaan juga menjadi faktor penentu keberhasilan eksperimen. Maka, variabel hari juga didefinisikan. Variabel agen juga perlu didefinisikan untuk menjadi memetakan variabel konsentrasi yang telah dibuat. Terakhir, semua variabel dimasukkan ke dalam sebuah tabel dengan nama dataku seperti ditunjukkan di bawah:
konsentrasi = c(c_A,c_B,c_C,c_D)
hari = rep(c("1","2","3"),4)
agen = rep(c("A","B","C","D"),each=3)
dataku = data.table(hari,agen,konsentrasi)
dataku## hari agen konsentrasi
## 1: 1 A 84
## 2: 2 A 79
## 3: 3 A 83
## 4: 1 B 80
## 5: 2 B 77
## 6: 3 B 78
## 7: 1 C 83
## 8: 2 C 80
## 9: 3 C 80
## 10: 1 D 79
## 11: 2 D 79
## 12: 3 D 78Perhitungan ANOVA dengan mudah dilakukan menggunakan syntax seperti ditunjukkan di bawah:
aov2 = aov(konsentrasi~agen+hari, data=dataku)
summary(aov2)## Df Sum Sq Mean Sq F value Pr(>F)
## agen 3 28.667 9.556 5.831 0.0328 *
## hari 2 15.500 7.750 4.729 0.0585 .
## Residuals 6 9.833 1.639
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Dari uji ANOVA di atas, diperoleh bahwa faktor jenis logam pengelat (yang diberikan tanda bintang) memberikan pengaruh yang signifikan terhadap hasil eksperimen. Jika dibandingkan dengan waktu dilakukannya percobaan, ternyata kualitas logam pengelat jauh lebih penting. Tampak bahwa nilai \(Pr(>F)\) untuk agen pengelat lebih kecil dibandingkan pengaruh hari yang mencapai \(0.0585\). Jika nilai \(Pr(>F) < 0.05\), maka dikatakan bahwa faktor tersebut berpengaruh secara signifikan terhadap keberlangsungan eksperimen.