Membandingkan hasil panen bawang merah berdasarkan 3 faktor, diantaranya Luas Panen Habis, Produksi Habis, dan Luas Panen Belum Habis ( data diambil dari tahun 2018 ). Untuk mengetahui perbandingan hasil panen dari 3 faktor tersebut saya memilih 6 kecamatan secara acak untuk percobaan yang masing-masing diujikan dengan 3 faktor tersebut. Dengan menggunakan level signifikan 0,05, dilakukan pengujian apakah rata-rata hasil panen dari 3 faktor tersebut sama.
Proses Pengerjaan:
library(readxl)
data = read_xlsx ("statis.xlsx")
data
#Gabungkan baris data menjadi vektor dengan menggunakan fungsi t(matrix transpose) agar nilai urut
r = c(t(as.matrix(data)))
r
#Membuat variabel baru untuk kategori, level treatment, dan data treatment
f = c("LP Bongkar Habis", "Produksi Bongkar habis", "LP Belum Habis")
k = 3 #treatment levels
n = 6 #data tiap treatment
#Membuat factor level sesuai anggota variabel r dengan menggunakan fungsi gl
tm = gl(k, 1, n*k, factor(f))
tm
#Gunakan fungsi aov untuk menguji anova, kemudian print out hasil akhirnya menggunakan summary.
av = aov(r~tm)
summary(av)