利用数据算了一下方差等,但是感觉数挺大的,可能单位不对。
residuals = y_test.values - y_pred.flatten()
# 计算偏差的均值、标准差和方差
mean_residuals = round(np.mean(residuals), 4)
std_residuals = round(np.std(residuals), 4)
var_residuals = round(np.var(residuals), 4)
# 打印偏差的均值、标准差和方差
print("偏差均值:", mean_residuals)
print("偏差的标准差:", std_residuals)
print("偏差的方差:", var_residuals)
# 计算置信区间
alpha = 0.05 # 置信水平为95%
n = len(residuals) # 样本数
t_value = 2.042 # t分布的临界值
ci_lower = round(mean_residuals - t_value * std_residuals / np.sqrt(n), 4)
ci_upper = round(mean_residuals + t_value * std_residuals / np.sqrt(n), 4)
# 打印置信区间
print("偏差的置信区间 ({}%): ({}, {})".format((1 - alpha) * 100, ci_lower, ci_upper))