# Load library
library(tidyverse)

# Membaca data
data <- read.csv("data_penelitian1.csv")

# Cek data
head(data)

# Pastikan kolom faktor dikonversi ke faktor (kategori)
data$Aplikasi <- factor(data$Aplikasi)
data$Komposisi <- factor(data$Komposisi)

# Model ANOVA: Komposisi tersarang dalam Aplikasi
model <- aov(Jlh_daun ~ Aplikasi + Aplikasi:Komposisi, data = data)

# Ringkasan hasil ANOVA
summary(model)

# Panggil library agricolae
library(agricolae)

# Interaksi perlakuan
data$Interaksi <- interaction(data$Aplikasi, data$Komposisi)

# ANOVA dasar
anova_model <- aov(Jlh_daun ~ Interaksi, data = data)

# Uji DMRT (Duncan)
hasil_duncan <- duncan.test(anova_model, "Interaksi", group = TRUE)

# Lihat hasil huruf
hasil_duncan$groups
