Просмотр исходного текста страницы Генералов Кирилл N2

import numpy as np
import math
f = open("order.txt", "r")
X = []
Y = []
for line in f:
    q, x1, x2, x3, x4, x5, x6, y = map(float, line.split())
    x = []
    x.append(x1)
    x.append(x2)
    x.append(x3)
    x.append(x4)
    x.append(x5)
    x.append(x6)
    x.append(x1 * x1)
    x.append(x1 * x2)
    x.append(x1 * x3)
    x.append(x1 * x4)
    x.append(x1 * x5)
    x.append(x1 * x6)
    x.append(x2 * x2)
    x.append(x2 * x3)
    x.append(x2 * x4)
    x.append(x2 * x5)
    x.append(x2 * x6)
    x.append(x3 * x3)
    x.append(x3 * x4)
    x.append(x3 * x5)
    x.append(x3 * x6)
    x.append(x4 * x4)
    x.append(x4 * x5)
    x.append(x4 * x6)
    x.append(x5 * x5)
    x.append(x5 * x6)
    x.append(x6 * x6)
    X.append(x)
    Y.append(y)
f.close()
Mx = []
for j in range(27):
    Mx1 = np.sum(X[j]) / len(X)
    Mx.append(Mx1)
My = np.sum(Y)/len(Y)
print(Mx)
print("Мат.ожидание(Mx) =", Mx)
Dx = []
Dx11 = 0
Dy = 0


for i in range(len(X)):
    Dy += (Y[i] - My) ** 2
for i in range(27):
    Dxs = 0.0
    for j in range(len(X)):
        Dxs += (X[j][i] - Mx[i]) ** 2
    print(Dxs)

    Dx.append(Dxs/len(X))

print("Дисперсия(Dx) = ", Dx)
kx = []
for i in range(27):
    kxs = 0.0
    for j in range(len(X)):
        kxs += ((X[j][i] - Mx[i]) * (Y[j] - My))
    kx.append(kxs)
print("Коэффициент корреляции (kx)=", kx)
r = []
rABS = []
for i in range(len(kx)):
    r1 = (kx[i] / (len(X)-1)) / (math.sqrt(Dx[i]) * math.sqrt(Dy))
    r.append(r1)
    rABS.append(math.fabs(r1))
print("r =", r)
print(np.max(rABS))
print("Первый индекс = ", rABS.index(np.max(rABS)) + 1)