# ==============================
# Analisis Parameter Jumlah Anakan
# ==============================

# Load Data
data <- read.csv("data_penelitian1.csv")

# Konversi ke faktor
data$Aplikasi <- factor(data$Aplikasi)
data$Komposisi <- factor(data$Komposisi)

# Cek normalitas
shapiro.test(data$Jlh_anakan)

# Cek homogenitas
library(car)
leveneTest(Jlh_anakan ~ interaction(Aplikasi, Komposisi), data = data)

# ANOVA tersarang: Komposisi dalam Aplikasi
model_anakan <- aov(Jlh_anakan ~ Aplikasi + Aplikasi:Komposisi, data = data)
summary(model_anakan)

# Uji lanjut Duncan
library(agricolae)
data$Interaksi <- interaction(data$Aplikasi, data$Komposisi)
hasil_duncan_anakan <- duncan.test(aov(Jlh_anakan ~ Interaksi, data = data), "Interaksi", group = TRUE)

# Lihat hasil kelompok Duncan
hasil_duncan_anakan$groups
