end = start + WINDOW_LENGTH difference_data = np.array(difference_data) meanvar_data = np.array(meanvar_data) print( "[%s] Wavelet analysis on data done. Starting analysis on surrogates..." % (str(datetime.now()))) #if PHASE_ANALYSIS_YEAR == last_mid_year: # # (bins, cond_means, cond_means_surr, phase, dates, subtit = '', fname = None): # fn = ('debug/detail/%s_phase_bins_%d_time_point.png' % ('32to16' if WINDOW_LENGTH > 16000 else '16to14', last_mid_year)) # render_phase_and_bins(phase_bins, cond_means, cond_means_surrs, phase, # [g_working.get_date_from_ndx(0), g_working.get_date_from_ndx(-1)], fname = fn) if SURR_TYPE == 'AR': sg.prepare_AR_surrogates() if AMPLITUDE: sg_amp.prepare_AR_surrogates() def _cond_difference_surrogates(sg, sg_amp, a, a2, jobq, resq): mean, var, trend = a mean2, var2, trend2 = a2 last_mid_year = first_mid_year cond_means_out = np.zeros((8, )) while jobq.get() is not None: if SURR_TYPE == 'MF': sg.construct_multifractal_surrogates() sg.add_seasonality(mean, var, trend) if AMPLITUDE: sg_amp.construct_multifractal_surrogates()
plot_vars.append(cnt) start = g.find_date_ndx(date(y1 + cnt*WINDOW_SHIFT, 7, 28)) end = start + WINDOW_LENGTH difference_data = np.array(difference_data) meanvar_data = np.array(meanvar_data) print("[%s] Wavelet analysis on data done. Starting analysis on surrogates..." % (str(datetime.now()))) #if PHASE_ANALYSIS_YEAR == last_mid_year: # # (bins, cond_means, cond_means_surr, phase, dates, subtit = '', fname = None): # fn = ('debug/detail/%s_phase_bins_%d_time_point.png' % ('32to16' if WINDOW_LENGTH > 16000 else '16to14', last_mid_year)) # render_phase_and_bins(phase_bins, cond_means, cond_means_surrs, phase, # [g_working.get_date_from_ndx(0), g_working.get_date_from_ndx(-1)], fname = fn) if SURR_TYPE == 'AR': sg.prepare_AR_surrogates() if AMPLITUDE: sg_amp.prepare_AR_surrogates() def _cond_difference_surrogates(sg, sg_amp, a, a2, jobq, resq): mean, var, trend = a mean2, var2, trend2 = a2 last_mid_year = first_mid_year cond_means_out = np.zeros((8,)) while jobq.get() is not None: if SURR_TYPE == 'MF': sg.construct_multifractal_surrogates() sg.add_seasonality(mean, var, trend) if AMPLITUDE: sg_amp.construct_multifractal_surrogates()