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)) df_info_lr['phi_lr'] = phi_lr print(df_info_lr)