Tablas de correlacion de estadistica en python

#tablas de correlacion
print("----------------TABLAS DE CORRELACIÓN-------------------------")
#obtenemos el tamaño del vector
tama = int(input("digite la cantidad de datos que va a ingresar"))
datos = [0]*tama #creamos el vector de la cantidad que tiene la variable tama
datos2 = [0]*tama
i= 0
SumatoriaX2=0
SumatoriaY2=0
SumatoriaY=0
SumatoriaX=0

sx2=0#sumatoria de x - xpromedio todo al cuadrado
sy2=0
sxy=0
#terminamos de definir algunas variables
############################################################

print ("ingrese datos de x para hallar correlacion")
while i < tama:
    datos[i]=float (input ("dato: "))
    SumatoriaX += datos[i] #sumamos todos los x
    SumatoriaX2+=(datos[i]**2)
    i+=1
print ("LOS DATOS INGRESADOS SON ",datos)

##########################################################

print ("ingrese datos de y ")
print ("ingrese sus datos")
Xpromedio= SumatoriaX/tama
i =0
while i < tama:
    datos2[i]=float (input ("dato: "))
    SumatoriaY+=datos2[i]
    SumatoriaY2+=(datos[i]**2)
    i+=1
Ypromedio= SumatoriaY/tama
i=0
print("-----------------------------------------------------------------------------------------------")


while i < tama:
    print("--------intervalo ",i,"----------------------------------------------------------------")
    print("|x| ",datos[i])
    print("|y| ",datos2[i])
    print("|X*Y| ",datos[i]*datos2[i])
    print("|Xª2| ",datos[i]**2)
    print("|Yª2| ",datos2[i]**2)
    print("|x -X|",datos[i]-Xpromedio)
    print("|y-Y|",datos2[i]-Ypromedio)
    print("|x -X|ª2",(datos[i]-Xpromedio)**2)
    print("|y-Y|ª2",(datos2[i]-Ypromedio)**2)
    sx2+= (datos[i]-Xpromedio)**2
    sy2+= (datos2[i]-Ypromedio)**2
    sxy += datos[i]*datos2[i]  
    i+=1
print("-----------------------------------------------------------------------------------------------")
print("Fin Tabla")
print("-----------------------------------------------------------------------------------------------")
sx=(sx2/(tama-1))**(1/2)
print("sx=",sx)

#la raiz cuadrada de la sumatoria de x-xpromedio elevado al cuadrado dividido el numero de elementos
sy=(sy2/(tama-1))**(1/2)
b =((tama*sxy) - (SumatoriaX*SumatoriaY))/((tama*SumatoriaX2) - (SumatoriaX**2))
a = (SumatoriaY/tama)- (b*SumatoriaX/tama)
sxyy=(sxy/(tama))-(Xpromedio*Ypromedio)
print("sy=",sy)
print("sumatoria xy=",sxy)
print("promedio x",Xpromedio)
print("promedio y",Ypromedio)

r= sxyy/(sx*sy)
print ("la ecuacion de la linea recta es ")
print ("Y = ",a ," + ",b,"X")
print ("la correlacion es de",r,"sxy= ", sxyy)
if r>(-1) and r<=(-0.9):
    print("!!!la correlacion es Excelente!!!!")
if r>=(0.8) and r<=(0.9):
    print("!!!la correlacion es buena!!!")
if r>(0.6) and r<(0.8):
    print("!!!la correlacion es regular!!!")
if r>(0.3) and r<(0.6):
    print("!!!mala!!!")