예제 #1
0
# 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))