可以采用两种方法,一种是直接调用像seaborn这样的库,另一种是在matplotlib的基础上根据置信带的原理自己完善。
一、调用seaborn
可以使用sns.regplot()这个函数,其中参数ci为置信水平,默认为95%,我们可以设置为99%或者其他值。
调用方式为:
sns.regplot(x=x, y=y, ci=95)
具体示例如下所示:
代码
# Import standard packages
import seaborn as sns
import numpy as np
import matplotlib.pyplot as plt
# Set time seed
np.random.seed(8)
# Generate data
mean, cov = [4, 6], [(1.5, .7), (.7, 1)]
x, y = np.random.multivariate_normal(mean, cov, 80).T
# Plot figure
ax = sns.regplot(x=x, y=y, ci=95)
plt.show()
lr.png
二、自己实现
在理解置信区间(confidence interval)的基础上,尝试实现
例如下面的代码中