Estoy tratando de obtener un intervalo de confianza en un ajuste exponencial para algunosx,y datos disponiblesaquí) Aquí está el MWE que tengo para encontrar el mejor ajuste exponencial a los datos:

from pylab import *
from scipy.optimize import curve_fit

# Read data.
x, y = np.loadtxt('exponential_data.dat', unpack=True)

def func(x, a, b, c):
    '''Exponential 3-param function.'''
    return a * np.exp(b * x) + c

# Find best fit.
popt, pcov = curve_fit(func, x, y)
print popt

# Plot data and best fit curve.
scatter(x, y)
x = linspace(11, 23, 100)
plot(x, func(x, *popt), c='r')
show()

que produce:

¿Cómo puedo obtener el intervalo de confianza del 95% (o algún otro valor) en este ajuste, preferiblemente usando puropython, numpy oscipy (¿Cuáles son los paquetes que ya he instalado)?