Intervalle de confiance pour LOWESS en Python

comment calculer les intervalles de confiance pour une régression LOWESS en Python? Je voudrais les ajouter comme une région ombrée au graphique de LOESS créé avec le code suivant (d'autres paquets que des modèles de statistiques sont aussi bien).

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()

j'ai ajouté un exemple de graphique avec intervalle de confiance ci-dessous à partir du webblog Statistiques Sérieuses (il est créé en utilisant ggplot dans R).