L1_Surrogates #85

SOI = pyleo.utils.load_dataset('SOI')
psd = SOI.standardize().spectral(method='mtm')
psd.plot()

fig, axs = plt.subplots(1,2,figsize=(10,4),sharey=True)
axs = axs.flatten()
psd_s = {}
for i, method in enumerate(['ar1sim','CN']):
    psd_s[method] = psd.signif_test(method=method,number=1000)
    psd_s[method].plot(ax=axs[i])
    axs[i].set_title(method, fontweight='bold')