# ==============================
# 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$Tinggi)

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

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

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

# Lihat hasil kelompok Duncan
hasil_duncan_tinggi$groups

