library(readxl)
## Warning: package 'readxl' was built under R version 4.3.2
suplement<-read_xlsx("D:/rancob/Tugas_Kelompok.xlsx")
suplement
## # A tibble: 24 × 2
##    Pabrik Kemampuan
##    <chr>      <dbl>
##  1 A           6.78
##  2 A           6.51
##  3 A           6.92
##  4 A           7.08
##  5 B           6.88
##  6 B           6.13
##  7 B           7.21
##  8 B           7.84
##  9 C           9.8 
## 10 C           9.72
## # ℹ 14 more rows

Kasus & Hipotesis

  1. Apakah terdapat perbedaan Bobot ayam dari sumplement pabrik A dengan pabrik B, C, D, E, F. \[ H_0 : 5\mu_a-\mu_b-\mu_c-\mu_d-\mu_e-\mu_f=0\\ H_1 : 5\mu_a-\mu_b-\mu_c-\mu_d-\mu_e-\mu_f≠0\\ \]
  2. Apakah terdapat perbedaan Bobot ayam dari suplement pabrik B dengan pabrik C, D, E, F. \[ H_0 : 4\mu_b-\mu_c-\mu_d-\mu_e-\mu_f=0\\ H_1 : 4\mu_b-\mu_c-\mu_d-\mu_e-\mu_f≠0\\ \]
  3. Apakah terdapat perrbedaan Bobot ayam dari suplement pabrik pabrik C dengan pabrik D, E, F. \[ H_0 : 3\mu_c-\mu_d-\mu_e-\mu_f=0\\ H_1 : 3\mu_c-\mu_d-\mu_e-\mu_f≠0\\ \]
  4. Apakah terdapat perbedaan Bobot ayam dari suplement pabrik D dengan pabrik E, F. \[ H_0 : 2\mu_d-\mu_e-\mu_f=0\\ H_1 : 2\mu_d-\mu_e-\mu_f≠0\\ \]
  5. Apakah terdapat perrbedaan Bobot ayam dari suplement pabrik E dengan pabrik F. \[ H_0 : \mu_e-\mu_f=0\\ H_1 : \mu_e-\mu_f≠0\\ \]

Buat Kontras

suplement$Pabrik<-as.factor(suplement$Pabrik)
levels(suplement$Pabrik)
## [1] "A" "B" "C" "D" "E" "F"
contrasts(suplement$Pabrik) <- cbind(c(5, -1, -1, -1, -1, -1), c(0,4,-1,-1,-1,-1), c(0,0,3,-1,-1,-1), c(0,0,0,2,-1,-1), c(0,0,0,0,1,-1))
contrasts(suplement$Pabrik)
##   [,1] [,2] [,3] [,4] [,5]
## A    5    0    0    0    0
## B   -1    4    0    0    0
## C   -1   -1    3    0    0
## D   -1   -1   -1    2    0
## E   -1   -1   -1   -1    1
## F   -1   -1   -1   -1   -1

Anova dengan Kontras Orthogonal

aovKontras<-aov(Kemampuan~Pabrik, suplement)
summary.aov(aovKontras, split = list (Pabrik = list('A vs BCDEF' = 1, 'B vs CDEF' = 2, 'C vs DEF' = 3, 'D vs EF' = 4, 'E vs F' = 5)))
##                      Df Sum Sq Mean Sq F value   Pr(>F)    
## Pabrik                5 25.654   5.131  11.326 4.72e-05 ***
##   Pabrik: A vs BCDEF  1  3.387   3.387   7.476   0.0136 *  
##   Pabrik: B vs CDEF   1  3.325   3.325   7.340   0.0144 *  
##   Pabrik: C vs DEF    1 16.579  16.579  36.597 1.02e-05 ***
##   Pabrik: D vs EF     1  0.019   0.019   0.041   0.8413    
##   Pabrik: E vs F      1  2.344   2.344   5.173   0.0354 *  
## Residuals            18  8.154   0.453                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Dapat disimpulkan bahwa rata-rata bobot ayam antara perbandingan hasil produksi suplemen dari pabrik A dengan pabrik B,C,D,E,F, pabrik B dengan pabrik C,D,E,F, pabrik C dengan pabrik D,E,F, dan pabrik E dengan F semuanya berbeda signifikan, kecuali rata-rata bobot ayam yang berasal dari hasil produksi suplemen pabrik D dengan pabrik E,F tidak berbeda signifikan.