Python 统计学：计算 R 平方值并保存到文件

import numpy as np

c_pred = [1.2] c_true = [1.1] d_pred = [1.5] d_true = [1.6] e_pred = [2.3] e_true = [2.1] f_pred = [3.1] f_true = [3.0] g_pred = [4.7] g_true = [4.8] h_pred = [5.2] h_true = [5.1] i_pred = [6.3] i_true = [6.4] j_pred = [7.5] j_true = [7.6] k_pred = [8.9] k_true = [9.0]

R_squared = 1 - ((c_pred[0] - c_true[0]) ** 2 + (d_pred[0] - d_true[0]) ** 2 + (e_pred[0] - e_true[0]) ** 2 + (f_pred[0] - f_true[0]) ** 2 + (g_pred[0] - g_true[0]) ** 2 + (h_pred[0] - h_true[0]) ** 2 + (i_pred[0] - i_true[0]) ** 2 + (j_pred[0] - j_true[0]) ** 2 + (k_pred[0] - k_true[0]) ** 2) ** 2 / ((c_pred[0] - 1 / c_true[0]) ** 2 + (d_pred[0] - 1 / d_true[0]) ** 2 + (e_pred[0] - 1 / e_true[0]) ** 2 + (f_pred[0] - 1 / f_true[0]) ** 2 + (i_pred[0] - 1 / i_true[0]) ** 2 + (j_pred[0] - 1 / j_true[0]) ** 2 + (k_pred[0] - 1 / k_true[0]) ** 2) ** 2

print('R^2的值为:', R_squared)

B1 = np.array([R_squared]) np.savetxt("R_squared.txt", B1, fmt="%f")