Как бы я рассчитал доверительные интервалы для регрессии LOWESS в Python? Я хотел бы добавить их как заштрихованную область в график LOESS, созданный с помощью следующего кода (другие пакеты, кроме statsmodels, тоже подойдут).

import numpy as np
import pylab as plt
import statsmodels.api as sm

x = np.linspace(0,2*np.pi,100)
y = np.sin(x) + np.random.random(100) * 0.2
lowess = sm.nonparametric.lowess(y, x, frac=0.1)

plt.plot(x, y, '+')
plt.plot(lowess[:, 0], lowess[:, 1])
plt.show()

Я добавил пример графика с доверительным интервалом ниже из веб-блогаСерьезная статистика (он создается с использованием ggplot в R).