# 1. Input data manual
data_bayi <- data.frame(
  Y  = c(57.5, 52.8, 61.3, 67.0, 53.5, 62.7, 56.2, 68.5, 69.2),
  X1 = c(78, 69, 77, 88, 67, 80, 74, 94, 102),
  X2 = c(2.75, 2.15, 4.41, 5.52, 3.21, 4.32, 2.31, 4.30, 3.71),
  X3 = c(29.5, 26.3, 32.2, 36.5, 27.2, 27.7, 28.3, 30.3, 28.7)
)

# 2. Hitung matriks korelasi
matriks_korelasi <- cor(data_bayi)

# 3. Install package jika belum ada
if (!require("pheatmap")) install.packages("pheatmap")
## Loading required package: pheatmap
## Warning: package 'pheatmap' was built under R version 4.5.3
# 4. Panggil library & buat heatmap
library(pheatmap)
pheatmap(matriks_kriks_korelasi <- matriks_korelasi, 
         display_numbers = TRUE, 
         cluster_rows = FALSE, 
         cluster_cols = FALSE)

# =========================================================
# 5. TAMBAHAN: SINTAKS INTERPRETASI DI CONSOLE
# =========================================================

# Keterangan nama variabel lengkap
ket_var <- c(
  Y  = "Panjang bayi saat lahir (cm)",
  X1 = "Umur (hari)",
  X2 = "Berat badan (kg)",
  X3 = "Lingkar dada (cm)"
)

# Fungsi kategorisasi nilai korelasi
kategori_korelasi <- function(r) {
  abs_r <- abs(r)
  if (abs_r == 1) return("Sempurna")
  else if (abs_r >= 0.80) return("Sangat Kuat")
  else if (abs_r >= 0.60) return("Kuat")
  else if (abs_r >= 0.40) return("Sedang")
  else if (abs_r >= 0.20) return("Lemah")
  else return("Sangat Lemah")
}

# Cetak interpretasi di Console
cat("\n--- INTERPRETASI KORELASI ---\n\n")
## 
## --- INTERPRETASI KORELASI ---
vars <- colnames(matriks_korelasi)

for (i in 1:(length(vars)-1)) {
  for (j in (i+1):length(vars)) {
    v1 <- vars[i]
    v2 <- vars[j]
    val <- round(matriks_korelasi[i, j], 2)
    kat <- kategori_korelasi(val)
    
    cat(sprintf("- Hubungan %s (%s) dan %s (%s) = %.2f: Korelasi Positif %s\n", 
                ket_var[v1], v1, ket_var[v2], v2, val, kat))
  }
}
## - Hubungan Panjang bayi saat lahir (cm) (Y) dan Umur (hari) (X1) = 0.95: Korelasi Positif Sangat Kuat
## - Hubungan Panjang bayi saat lahir (cm) (Y) dan Berat badan (kg) (X2) = 0.76: Korelasi Positif Kuat
## - Hubungan Panjang bayi saat lahir (cm) (Y) dan Lingkar dada (cm) (X3) = 0.56: Korelasi Positif Sedang
## - Hubungan Umur (hari) (X1) dan Berat badan (kg) (X2) = 0.53: Korelasi Positif Sedang
## - Hubungan Umur (hari) (X1) dan Lingkar dada (cm) (X3) = 0.39: Korelasi Positif Lemah
## - Hubungan Berat badan (kg) (X2) dan Lingkar dada (cm) (X3) = 0.78: Korelasi Positif Kuat