calc_PhiID return array

[Lasa937]

Attempting to calculate information atoms on two timeseries derived from BOLD signals, the calc_PhiID returns an array of atoms which has N - tau dimensions, N being the number of timepoints in the timeseries. My understanding would be that the function would return a single value per information atom instead, is that correct?

Hereby attached the steps required to elicit the same behavior:

path_file = '/XXX/preprocessing/ROI_Subject001_Condition001.mat'
mat_file = loadmat(path_file)

data = mat_file['data']
data_src = data[0,3]
data_trg = data[0,4]

aaa, bbb = calc_PhiID(data_src, data_trg, tau=1, kind='gaussian', redundancy='MMI')

where both src and trg are (152,1) and the return information atom is (151,1)

[liuzhenqi77]

Hey, thanks for the issue. Currently, the function returns the "local" PhiID atoms. You can take their mean to get the single value. We may make changes in future versions.

[Lasa937]

Thank you!!