ANOVA

Deskripsi Data

Bawang merah merupakan salah saru komoditas sayuran unggulan yang sejak lama telah diusahakan oleh petani secara intensif.Permintaan akan bawang merah terus meningkat sejalan peningkatan jumlah penduduk.Pada data dinas pertanian dan ketahanan pangan kab.brebes terdapat laporan hasil tanaman sayuran (bawang merah). Dari data tersebut terdapat hasil laporan selama satu tahun (2018). Dari data tersebut terdapat 18 kecamatan. Masing-masing kecamatan terdapat tiga faktor.
Faktor pertama yaitu Luas panen dibongkar habis.
Faktor kedua yaitu Produksi dibongkar habis.
Faktor yang ketiga Luas panen belum habis.
Dari ketiga faktor tersebut dijadikan variabel.
Variabel pertama LPDH (Luas Panen Dibongkar Habis), variabel kedua PrBH (Produksi Dibongkar Habis), variabel LPBH (Luas Panen Belum Habis). Dari 18 kecamatan tersebut. Saya mengambil beberapa sampel kecamatan (saya ambil 5 kecamatan secara acak yaitu Kec. Bumiayu, Kec.Larangan, Kec.rebes, Kec.Tonjong, Kec.Kersana). Saya melakukan perhitungan rata-rata perkecamatan tiap variabel/faktor.

Dengan menggunakan level signifikan 0,05 lakukan pengujian apakah rata-rata dari ketiga faktor tersebut sama.

library(readxl)
data = read_xlsx ("dinas.xlsx")
data
## # A tibble: 5 x 3
##    LPDH  PrDH  LPBH
##   <dbl> <dbl> <dbl>
## 1     0     0     0
## 2   583 57371   105
## 3   334 33003    97
## 4     0     0     0
## 5    98  9381    78

#Gabungkan baris data menjadi vektor. gunakan fungsi t(matrix tranpose) agar nilai urut

r = c(t(as.matrix(data)))
r
##  [1]     0     0     0   583 57371   105   334 33003    97     0     0     0
## [13]    98  9381    78

#Membuat variabel baru

f = c("LPDH","PrDH", "LPBH")
k = 3
n = 5
tm = gl(k, 1, n*k, factor(f))
tm
##  [1] LPDH PrDH LPBH LPDH PrDH LPBH LPDH PrDH LPBH LPDH PrDH LPBH LPDH PrDH LPBH
## Levels: LPDH PrDH LPBH

#Gunakan aov untuk uji anova

av = aov(r~tm)
summary(av)
##             Df    Sum Sq   Mean Sq F value Pr(>F)  
## tm           2 1.310e+09 654846782    3.17 0.0785 .
## Residuals   12 2.479e+09 206557221                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Karena p value(0.0785) lebih besar dari level signifikan(0,05) maka hipotesis diterima (dengan setara), sehingga rata-rata dari 3 faktor tersebut sama.