如何计算和绘制线性回归的系数区间？

Question

我需要用线性回归绘制系数区间。

    model = LinearRegression()
    model.fit(np.array(X).reshape((-1,1)),np.array(Y))
    print(model.score(np.array(X).reshape((-1,1)),np.array(Y)))
    print('Coefficient: \n', model.coef_)
    print(model.intercept_)
    plt.scatter(np.array(X), range(X_max, X_max + len(X)), marker='o', label=train_key, alpha=0.05, color='red') 

我需要提取一个置信区间，而不仅仅是一个值，并像这张图那样绘制它。