Aquí está mi problema:polyfit no toma valores de fecha y hora, por lo que convertí fecha y hora conmktime produciendo los trabajos de ajuste polinomial

z4 = polyfit(d, y, 3) 
p4 = poly1d(z4)

Sin embargo, para la trama, me gustaría la descripción de fecha y hora en el eje y no # descubrí cómo hacerlo. ¿Me puedes ayudar?

fig = plt.figure(1)
cx= fig.add_subplot(111) 

xx = linspace(0,  d[3], 100)
pylab.plot(d, y, '+', xx, p4(xx),'-g')
cx.plot(d, y,'+', color= 'b', label='blub')
plt.errorbar(d, y,
           yerr,
           marker='.',
           color='k',
           ecolor='b',
           markerfacecolor='b',
           label="series 1",
           capsize=0,
           linestyle='')

cx.grid()
cx.set_ylim(0,0.03)
plt.show()

resto del código:

import numpy as np
import matplotlib.pyplot as plt
from matplotlib import axis
from datetime import datetime
from numpy import *
import pylab
import time

mis primeros 4 puntos de datos de tiempo

x = [datetime(1978, 7, 7), 
     datetime(1980, 9, 26), 
     datetime(1983, 8, 1), 
     datetime(1985,8,8)]

d=[]
for i in x:
    d.append(time.mktime(i.timetuple()))

mis primeros 4 valores de datos

y = [0.00134328779552718,
     0.00155187668863844,
     0.0039431374327427,
     0.00780037563783297]

mis desviaciones estándar calculadas para barras de error

yerr = [0.0000137547160254577,
        0.0000225670232594083,
        0.000105623642510075,
        0.00011343121508]