Menghitung korelasi secara akurat dengan metode numerik

Berikut adalah algoritma (dalam pseudocode) yang akan mengestimasi korelasi dengan menggunakan metode mumerik

sum_sq_x = 0

sum_sq_y = 0

sum_coproduct = 0

mean_x = x[1]

mean_y = y[1]

last_x = x[1]

last_y = y[1]

for i in 2 to N:

    sweep = (i - 1.0) / i

    delta_x = x[i] - mean_x

    delta_y = y[i] - mean_y

    sum_sq_x += delta_x * delta_x * sweep

    sum_sq_y += delta_y * delta_y * sweep

    sum_coproduct += delta_x * delta_y * sweep

    mean_x += delta_x / i

    mean_y += delta_y / i 

pop_sd_x = sqrt( sum_sq_x / N )

pop_sd_y = sqrt( sum_sq_y / N )

cov_x_y = sum_coproduct / N

correlation = cov_x_y / (pop_sd_x * pop_sd_y)