Como eu calcularia os intervalos de confiança para uma regressão LOWESS em Python? Gostaria de adicioná-los como uma região sombreada ao gráfico LOESS criado com o código a seguir (outros pacotes além do statsmodels também são bons).

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

Adicionei um exemplo de plotagem com intervalo de confiança abaixo no blogEstatísticas sérias (é criado usando ggplot em R).