pam_phi = pam.iloc[:, 0].values df_ams_ssn = pd.read_csv(r'../ssn_flux_2/output/info_test.csv', header=0, index_col=0) ams_ssn = df_ams_ssn.loc['2011-06-01':'2017-03-01', 'ssn_delay'].values + fix_ssn ams_phi = ams.iloc[:, 0].values pam_log_ssn = pam_ssn**0.5 + fix_log_ssn ams_log_ssn = ams_ssn**0.5 + fix_log_ssn # print(pam_ssn) pam_end = df_alpha.loc[:'2014-01', alpha_type] pam_dealy = np.ones(len(pam_ssn), dtype='int') * -12 pam_alpha = FUNC.delay_ssn(pam_end, pam_dealy) + fix_alpha ams_end = df_alpha.loc[:'2017-03', alpha_type] ams_dealy = np.ones(len(ams_ssn), dtype='int') * -12 ams_alpha = FUNC.delay_ssn(ams_end, ams_dealy) + fix_alpha print(ams_alpha) # pam_alpha = np.log(pam_alpha) # ams_alpha = np.log(ams_alpha) df_pam = pd.DataFrame(np.vstack((pam_phi, pam_ssn, pam_log_ssn, pam_alpha)).T, columns=['phi_0', 'ssn', 'log_ssn_fix', alpha_type], index=df_pam_ssn.loc['2006-08-01':'2014-01-01', :].index) df_ams = pd.DataFrame(np.vstack((ams_phi, ams_ssn, ams_log_ssn, ams_alpha)).T, columns=['phi_0', 'ssn', 'log_ssn_fix', alpha_type], index=df_ams_ssn.loc['2011-06-01':'2017-03-01', :].index)
df_info = pd.read_csv(r'../ssn_flux_2/output/info_test.csv', header=0, index_col=0) df_alpha = pd.read_csv(r'../mon_alpha_1/output/alpha_mon.csv', header=0, index_col=0) print(df_alpha) phi_c = df_info.loc['1977-06-01':'2020-01-01', 'phi_c'].values Phi_c = df_info.loc['1977-06-01':'2020-01-01', 'Phi_c'].values ssn_delay = (df_info.loc['1977-06-01':'2020-01-01', 'ssn_delay'].values + fix_ssn) ** 0.5 + fix_log_ssn print(ssn_delay) alpha = df_alpha.loc[:, alpha_type].values dealy = np.ones(len(alpha) - 12, dtype='int') * -12 # dealy[:24] = np.ones(24, dtype='int') * -8 alpha_dealy = FUNC.delay_ssn(alpha, dealy) + fix_alpha A_values = df_alpha.loc['1977-05-01':, 'A'].values # alpha_dealy = np.log(alpha_dealy) print(len(A_values)) mark = df_info.loc['1977-06-01':'2020-01-01', 'mark'].values para = np.vstack((ssn_delay, alpha_dealy)).T phi_lr = np.array([]) # for i in range(len(mark)): # phi = FUNC.LR_all(para[i]) # phi_lr = np.append(phi_lr, phi)