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()
