# Input data sesuai dengan tabel yang diberikan
data <- data.frame(
  Urutan = factor(rep(1:4, each=4)),    
  Operator = factor(rep(1:4, times=4)), 
  Metode = factor(c("C", "B", "D", "A", "B", "C", "A", "D", "A", "D", "B", "C", "D", "A", "C", "B")),
  Tempat = factor(c("β", "y", "δ", "α", "α", "δ", "y", "β", "δ", "α", "β", "y", "y", "β", "α", "δ")),
  Waktu = c(11, 10, 14, 8, 8, 12, 10, 12, 9, 11, 7, 15, 9, 8, 18, 6)    
)

# Lakukan analisis variansi (ANOVA) Graeco-Latin
anova_graeco_latin <- aov(Waktu ~ Urutan + Operator + Metode + Tempat, data=data)

# Lihat hasil ANOVA
summary(anova_graeco_latin)
##             Df Sum Sq Mean Sq F value Pr(>F)
## Urutan       3    0.5    0.17   0.018  0.996
## Operator     3   19.0    6.33   0.691  0.616
## Metode       3   95.5   31.83   3.473  0.167
## Tempat       3    7.5    2.50   0.273  0.843
## Residuals    3   27.5    9.17
# Hasil dalam bentuk tabel
anova_table <- summary(anova_graeco_latin)
anova_table
##             Df Sum Sq Mean Sq F value Pr(>F)
## Urutan       3    0.5    0.17   0.018  0.996
## Operator     3   19.0    6.33   0.691  0.616
## Metode       3   95.5   31.83   3.473  0.167
## Tempat       3    7.5    2.50   0.273  0.843
## Residuals    3   27.5    9.17