如果你只是要求数据V>0.5的概率,我会找到它的独立概率。正如我所评论的,我不认为数据之间存在相关性,因为data-V有一个峰值,而data-R根本没有反应。
import numpy as np
types = {"data_v": "data_v.csv"} # only considering data-V
for protname, fname in types.items(): # your code to load the data
col_time, col_window = np.loadtxt(fname, delimiter=',').T
trailing_window = col_window[:-1] # "past" values at a given index
leading_window = col_window[1:] # "current values at a given index
decreasing_inds = np.where(leading_window < trailing_window)[0]
quotient = leading_window[decreasing_inds] / trailing_window[decreasing_inds]
quotient_times = col_time[decreasing_inds]
# Now we'll go through the data for quotients and find how many meet the > 0.5 criteria
count = 0
for quot in quotient:
if quot > 0.5:
count +