Uji T Satu Sampel: Evaluasi Kadar COD PT. Tirta Jernih

Memasukkan Data

# Data hasil pengukuran
n <- 9                  # jumlah sampel
rata_sampel <- 84       # rata-rata kadar COD (mg/L)
sd_sampel <- 9          # simpangan baku (mg/L)
mu0 <- 80               # rata-rata populasi yang diujikan (batas maksimum)
alpha <- 0.05           # tingkat signifikansi

# Hitung t-statistik
t_hitung <- (rata_sampel - mu0) / (sd_sampel / sqrt(n))

# Derajat bebas
df <- n - 1

# Hitung p-value (satu sisi, karena ingin menguji apakah melebihi 80 mg/L)
p_value <- pt(t_hitung, df = df, lower.tail = FALSE)

Menghitung Nilai T-Kritis

# t-kritis untuk uji satu sisi (kanan)
t_kritis <- qt(1 - alpha, df = df)

Interpretasi Hasil

cat("--- Hasil Uji & Interpretasi Lengkap ---\n\n")

## --- Hasil Uji & Interpretasi Lengkap ---

cat(sprintf("Rata-rata Sampel (x̄): %.2f mg/L\n", rata_sampel))

## Rata-rata Sampel (x̄): 84.00 mg/L

cat(sprintf("Simpangan Baku Sampel (s): %.2f mg/L\n", sd_sampel))

## Simpangan Baku Sampel (s): 9.00 mg/L

cat(sprintf("Jumlah Sampel (n): %d\n", n))

## Jumlah Sampel (n): 9

cat("--------------------------------------------\n")

## --------------------------------------------

cat("Hipotesis Nol (H0): μ ≤ 80 mg/L\n")

## Hipotesis Nol (H0): μ ≤ 80 mg/L

cat("Hipotesis Alternatif (H1): μ > 80 mg/L\n")

## Hipotesis Alternatif (H1): μ > 80 mg/L

cat("--------------------------------------------\n")

## --------------------------------------------

cat(sprintf("Tingkat Signifikansi (α): %.3f\n", alpha))

## Tingkat Signifikansi (α): 0.050

cat(sprintf("T-Hitung (t-statistic): %.4f\n", t_hitung))

## T-Hitung (t-statistic): 1.3333

cat(sprintf("T-Kritis (t-critical): %.4f\n", t_kritis))

## T-Kritis (t-critical): 1.8595

cat(sprintf("P-value: %.4f\n", p_value))

## P-value: 0.1096

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.1096) > alpha (0.050), 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.3333) <= T-Kritis (1.8595), maka Hipotesis Nol GAGAL DITOLAK.

# Kesimpulan akhir
cat("\nKESIMPULAN:\n")

## 
## KESIMPULAN:

if (p_value <= alpha) {
  cat("Dengan tingkat kesalahan 5%, terdapat cukup bukti statistik bahwa kadar COD melebihi ambang batas 80 mg/L. BLHD dapat mempertimbangkan tindakan hukum.\n")
} else {
  cat("Dengan tingkat kesalahan 5%, tidak terdapat cukup bukti statistik bahwa kadar COD melebihi ambang batas 80 mg/L. BLHD belum memiliki dasar kuat untuk tindakan hukum.\n")
}

## Dengan tingkat kesalahan 5%, tidak terdapat cukup bukti statistik bahwa kadar COD melebihi ambang batas 80 mg/L. BLHD belum memiliki dasar kuat untuk tindakan hukum.