Memasukkan Data
# Data hasil panen masing-masing kelompok (dalam kg per petak)
bio_subur <- c(4.0, 5.5, 4.2, 5.3, 4.5, 5.0, 4.6, 5.1)
nutri_prima <- c(5.2, 4.5, 5.5, 4.8, 5.0, 5.3, 4.6, 5.4, 4.7, 5.0)
alpha <- 0.10
# Tampilkan data
data.frame(
Kelompok = c(rep("BioSubur", length(bio_subur)), rep("NutriPrima", length(nutri_prima))),
Hasil_Panen = c(bio_subur, nutri_prima)
)
## Kelompok Hasil_Panen
## 1 BioSubur 4.0
## 2 BioSubur 5.5
## 3 BioSubur 4.2
## 4 BioSubur 5.3
## 5 BioSubur 4.5
## 6 BioSubur 5.0
## 7 BioSubur 4.6
## 8 BioSubur 5.1
## 9 NutriPrima 5.2
## 10 NutriPrima 4.5
## 11 NutriPrima 5.5
## 12 NutriPrima 4.8
## 13 NutriPrima 5.0
## 14 NutriPrima 5.3
## 15 NutriPrima 4.6
## 16 NutriPrima 5.4
## 17 NutriPrima 4.7
## 18 NutriPrima 5.0
Interpretasi Hasil
cat("--- Hasil Uji & Interpretasi Lengkap ---\n\n")
## --- Hasil Uji & Interpretasi Lengkap ---
cat("Data Sampel BioSubur:", paste(bio_subur, collapse = ", "), "\n")
## Data Sampel BioSubur: 4, 5.5, 4.2, 5.3, 4.5, 5, 4.6, 5.1
cat("Data Sampel NutriPrima:", paste(nutri_prima, collapse = ", "), "\n")
## Data Sampel NutriPrima: 5.2, 4.5, 5.5, 4.8, 5, 5.3, 4.6, 5.4, 4.7, 5
cat(sprintf("Rata-rata BioSubur (x̄₁): %.2f\n", mean(bio_subur)))
## Rata-rata BioSubur (x̄₁): 4.78
cat(sprintf("Rata-rata NutriPrima (x̄₂): %.2f\n", mean(nutri_prima)))
## Rata-rata NutriPrima (x̄₂): 5.00
cat("--------------------------------------------\n")
## --------------------------------------------
cat("Hipotesis Nol (H0): μ_BioSubur ≥ μ_NutriPrima\n")
## Hipotesis Nol (H0): μ_BioSubur ≥ μ_NutriPrima
cat("Hipotesis Alternatif (H1): μ_BioSubur < μ_NutriPrima\n")
## Hipotesis Alternatif (H1): μ_BioSubur < μ_NutriPrima
cat("--------------------------------------------\n")
## --------------------------------------------
cat(sprintf("Tingkat Signifikansi (α): %.3f\n", alpha))
## Tingkat Signifikansi (α): 0.100
cat(sprintf("T-Hitung (t-statistic): %.4f\n", t_hitung))
## T-Hitung (t-statistic): -1.0310
cat(sprintf("T-Kritis (t-critical): %.4f\n", t_kritis))
## T-Kritis (t-critical): -1.3597
cat(sprintf("P-value: %.4f\n", p_value))
## P-value: 0.1619
cat("--------------------------------------------\n\n")
## --------------------------------------------
# Keputusan berdasarkan P-Value
cat("Keputusan berdasarkan P-Value:\n")
## Keputusan berdasarkan P-Value:
if (p_value <= alpha) {
cat(sprintf("Karena P-value (%.4f) <= alpha (%.3f), maka Hipotesis Nol DITOLAK.\n", p_value, alpha))
} else {
cat(sprintf("Karena P-value (%.4f) > alpha (%.3f), maka Hipotesis Nol GAGAL DITOLAK.\n", p_value, alpha))
}
## Karena P-value (0.1619) > alpha (0.100), maka Hipotesis Nol GAGAL DITOLAK.
# Keputusan berdasarkan T-Statistik
cat("\nKeputusan berdasarkan T-Statistik:\n")
##
## Keputusan berdasarkan T-Statistik:
if (t_hitung <= t_kritis) {
cat(sprintf("Karena T-Hitung (%.4f) <= T-Kritis (%.4f), maka Hipotesis Nol DITOLAK.\n", t_hitung, t_kritis))
} else {
cat(sprintf("Karena T-Hitung (%.4f) > T-Kritis (%.4f), maka Hipotesis Nol GAGAL DITOLAK.\n", t_hitung, t_kritis))
}
## Karena T-Hitung (-1.0310) > T-Kritis (-1.3597), maka Hipotesis Nol GAGAL DITOLAK.
# Kesimpulan Akhir
cat("\nKESIMPULAN:\n")
##
## KESIMPULAN:
if (p_value <= alpha) {
cat("Dengan tingkat kesalahan 10%, terdapat cukup bukti statistik bahwa rata-rata hasil panen tomat dengan pupuk BioSubur lebih rendah dibandingkan dengan NutriPrima.\n")
} else {
cat("Dengan tingkat kesalahan 10%, tidak terdapat cukup bukti statistik bahwa rata-rata hasil panen BioSubur lebih rendah dari NutriPrima.\n")
}
## Dengan tingkat kesalahan 10%, tidak terdapat cukup bukti statistik bahwa rata-rata hasil panen BioSubur lebih rendah dari NutriPrima.