Laporan Komstat

nilai_korelasi <- function(x, y) {
  # Menghitung komponen-komponen utama
  jumlah_data <- length(x)
  sum_x <- sum(x)
  sum_y <- sum(y)
  sum_produk_deviasi <- sum(x * y)
  sum_kuadrat_deviasi_x <- sum(x^2)
  sum_kuadrat_deviasi_y <- sum(y^2)
  
  # Menghitung korelasi Pearson
  pembilang <- (jumlah_data * sum_produk_deviasi) - (sum_x * sum_y)
  penyebut <- sqrt((jumlah_data * sum_kuadrat_deviasi_x - sum_x^2) * (jumlah_data * sum_kuadrat_deviasi_y - sum_y^2))
  korelasi <- pembilang / penyebut
  
  # Menampilkan hasil korelasi
  cat("Korelasi Pearson:", korelasi, "\n")
  
  return(korelasi)
}