Diskusi 2 Rancangan Percobaan [STAT4431]

Nama : Tengku Mashitah Crisanty NIM : 054983862

# ================================
# 1. Input Data
# ================================
perlakuan <- factor(rep(c("Tanpa_Emulsi", "Dengan_Emulsi"), each = 5))
respon <- c(6, 7, 6, 8, 8,
            9, 8, 11, 9, 9)

data <- data.frame(perlakuan, respon)
print(data)

##        perlakuan respon
## 1   Tanpa_Emulsi      6
## 2   Tanpa_Emulsi      7
## 3   Tanpa_Emulsi      6
## 4   Tanpa_Emulsi      8
## 5   Tanpa_Emulsi      8
## 6  Dengan_Emulsi      9
## 7  Dengan_Emulsi      8
## 8  Dengan_Emulsi     11
## 9  Dengan_Emulsi      9
## 10 Dengan_Emulsi      9

# ================================
# 2. Analisis Ragam (ANOVA)
# ================================
model <- aov(respon ~ perlakuan, data = data)
summary(model)

##             Df Sum Sq Mean Sq F value Pr(>F)  
## perlakuan    1   12.1    12.1      11 0.0106 *
## Residuals    8    8.8     1.1                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

# ================================
# 3. Rata-rata per perlakuan
# ================================
tapply(data$respon, data$perlakuan, mean)

## Dengan_Emulsi  Tanpa_Emulsi 
##           9.2           7.0

# ================================
# 4. Uji lanjut (jika diperlukan, misalnya LSD)
# ================================
# install.packages("agricolae") # jika belum terinstal
library(agricolae)
uji_lsd <- LSD.test(model, "perlakuan", p.adj="none")
print(uji_lsd)

## $statistics
##   MSerror Df Mean       CV  t.value     LSD
##       1.1  8  8.1 12.94826 2.306004 1.52963
## 
## $parameters
##         test p.ajusted    name.t ntr alpha
##   Fisher-LSD      none perlakuan   2  0.05
## 
## $means
##               respon      std r        se      LCL       UCL Min Max Q25 Q50
## Dengan_Emulsi    9.2 1.095445 5 0.4690416 8.118388 10.281612   8  11   9   9
## Tanpa_Emulsi     7.0 1.000000 5 0.4690416 5.918388  8.081612   6   8   6   7
##               Q75
## Dengan_Emulsi   9
## Tanpa_Emulsi    8
## 
## $comparison
## NULL
## 
## $groups
##               respon groups
## Dengan_Emulsi    9.2      a
## Tanpa_Emulsi     7.0      b
## 
## attr(,"class")
## [1] "group"

Karena p-value (0.0104) < 0.05 → tolak H0.atau Karena:

𝐹hitung=11.00>𝐹tabel(0.05;1,8)=5.32. Fhitung =11.00>Ftabel(0.05;1,8)=5.32 Maka H₀ ditolak dan H₁ diterima.Artinya ada pengaruh signifikan pemberian emulsi terhadap kekuatan regang beton. Kesimpulan sama: emulsi berpengaruh signifikan terhadap kekuatan regang beton.