Estatistica (TEC Disco Rigido MB vs POP Brasil)

x = [0.05,0.24,1.11,5,60,250,500,2000]
y = [51944397,70992343,94508583,120150573,146917459,169590693,190755799,236259570]

media_x = 0
media_y = 0

a = []
b = []
ab = []
soma_ab = 0
soma_aa = 0
soma_bb = 0

#Calculo Media de X
calculo = 0
for valor in x:
    calculo += valor
    media_x = calculo/(len(x))
print("Média de X:", media_x)

#Calculo Media de Y
calculo = 0
for valor in y:
  calculo += valor
  media_y = calculo/(len(y))
print("Média de Y:", media_y)

#Cálculo de A (X-MédiaX)
for valor in x:
  calculo = valor - media_x
  a.append(calculo)
print("A:", a)

#Cálculo de B (Y-MédiaY)
for valor in y:
  calculo = valor - media_y
  b.append(calculo)
print("B:", b)

#Cálculo de A*B e Somatorio A*B
for i in range(len(a)):
  calculo = a[i]*b[i]
  soma_ab += calculo
  ab.append(calculo)
print("A*B:", ab)
print("Somatório A*B:", soma_ab)

#Cálculo A*A
for i in range(len(a)):
  a[i] *= a[i]
  soma_aa += a[i]
print("A*A:",a)
print("Somatório A*A:", soma_aa)

#Cálculo B*B
for i in range(len(b)):
  b[i] *= b[i]
  soma_bb += b[i]
print("B*B:",b)
print("Somatório B*B:", soma_bb)

from math import sqrt
r = soma_ab/(sqrt(soma_aa)*sqrt(soma_bb))
print ("Pearson:", r)