Вот моя проблема: polyfit
не принимает значения datetime, поэтому я преобразовал datetime с mktime
, создавая работы с полиномиальным подходом
z4 = polyfit(d, y, 3)
p4 = poly1d(z4)
Для сюжета, однако, мне хотелось бы описать дату и время на оси и не # выяснить, как это сделать. Можете ли вы мне помочь?
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()
остальная часть кода:
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import axis
from datetime import datetime
from numpy import *
import pylab
import time
мои первые 4 временных точки данных
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()))
мои первые 4 значения данных
y = [0.00134328779552718,
0.00155187668863844,
0.0039431374327427,
0.00780037563783297]
мои рассчитанные стандартные отклонения для баров ошибок
yerr = [0.0000137547160254577,
0.0000225670232594083,
0.000105623642510075,
0.00011343121508]