# 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) for i, flag in enumerate(mark): if flag == -12: phi = FUNC.LR_negative(para[i]) phi_lr = np.append(phi_lr, phi) elif flag == -5: phi = FUNC.LR_postive(para[i]) phi_lr = np.append(phi_lr, phi) print(phi_lr) print(len(phi_lr)) """ 保存LR计算出的phi_lr 计算 flux_lr (log) """ df_info_lr = df_info.loc['1977-06-01':'2020-01-01', :] print(len(df_info_lr))