outputlist = inputlist/fac
except TypeError: #not an array
outputlist = [el/fac for el in inputlist]
return outputlist
isi = convert_tdelt(results['tau_valid'], units='minutes')
isi_hr = convert_tdelt(results['tau_valid'], units='hours')
tau_ax = np.arange(0,30*60,.2)
try:
from sklearn.neighbors.kde import KernelDensity
kern_type = 'epanechnikov'
kern_lab = '{0}{1} KDE'.format(kern_type[0].upper(), kern_type[1:])
kernel = KernelDensity(kernel=kern_type, bandwidth=60).fit(isi[:, np.newaxis])
kde_plot = np.exp(kernel.score_samples(tau_ax[:, np.newaxis]))
except ImportError:
from scipy import stats
kern_lab = 'Gaussian KDE'
kernel = stats.gaussian_kde(isi, bw_method='scott')
kde_plot = kernel.evaluate(tau_ax)
fig, ax = splot.set_target(None)
ax.hist(isi_hr, bins=np.arange(0,25,0.5), histtype='step', normed=True, lw=1.5, label='Binned Data')
ax.plot(tau_ax/60., kde_plot*60., lw=1.5, label=kern_lab)
ax.set_xlim([0,25])
ax.set_ylabel('Probability')
ax.set_xlabel(r'Inter-substorm Interval, $\tau$ [hours]') #raw string req'd (else \t in \tau becomes [tab]au
ax.legend()
fig.suptitle('MSM$_{Python}$: ' + '{0} (1998-2002)'.format(satname))
plt.show()
