dates = np.load('dates.npy', allow_pickle=True)
dates = dates.flatten()

print(type(dates))
print(np.shape(dates))

l = []
for i in range(len(dates)):
    l.append(datetime.datetime.strptime(dates[i], '%Y-%m-%d %H:%M'))

print(type(l))
print(type(l[0]))
print(np.shape(l))

l = mdates.date2num(l)

ax = sns.boxenplot(l)
ax.set(xlabel='MY X AXIS', ylabel='MY Y AXIS')
ax.yaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M'))
ax.yaxis.set_major_locator(mdates.YearLocator())
plt.savefig('BOXENPLOT.pdf', bbox_inches='tight')

plt.hist(l, bins=10)
axes = plt.gca()
axes.set_xlabel('MY X AXIS')
axes.set_ylabel('MY Y AXIS')
axes.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m'))
plt.gcf().autofmt_xdate()
plt.show()

ax = sns.violinplot(l)
ax.set(xlabel='MY X AXIS', ylabel='MY Y AXIS')
ax.yaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d %H:%M'))
ax.yaxis.set_major_locator(mdates.YearLocator())
plt.savefig('VIOLINPLOT.pdf', bbox_inches='tight')