In [71]:
import pandas
import numpy as np

In [72]:
%cd C:\Users\da.angulo39\Documents\MATLAB\kmc400_mars


C:\Users\da.angulo39\Documents\MATLAB\kmc400_mars

In [73]:
thr = 0.05

In [74]:
p_file="at_corrected_p.csv"
effect_file="at_conts.csv"
excel_file="atencion_mb.xlsx" 
excel_file_p="atencion_mb_ps.xlsx"
filter_name_char="_"

#p_file="co_corrected_p.csv"
#effect_file="co_conts.csv"
#excel_file="coordinacion_mb.xlsx" 
#excel_file_p="coordinacion_mb_ps.xlsx"
#filter_name_char="-"

In [75]:
ef = pandas.read_csv(effect_file,index_col=0)
ef.head()


Out[75]:
FEF-left_roi.mat FEF-right_roi.mat MT-left_roi.mat MT-right_roi.mat SP-Left_roi.mat SP-Right_roi.mat Unnamed: 7
subject
9 2.202807 1.870303 -0.130362 1.233423 2.787738 2.194587 NaN
15 0.794069 0.582817 -0.375399 0.493381 1.192017 1.568221 NaN
19 1.063118 0.280254 0.563899 0.060583 1.721303 1.666498 NaN
25 0.741467 1.688934 1.246920 2.502616 1.809077 2.404814 NaN
29 0.749212 0.277349 0.904308 1.247580 0.457264 0.723541 NaN

In [76]:
ps = pandas.read_csv(p_file,index_col=0)
ps.head()


Out[76]:
FEF-left_roi.mat FEF-right_roi.mat MT-left_roi.mat MT-right_roi.mat SP-Left_roi.mat SP-Right_roi.mat Unnamed: 7
subject
9 0.000047 0.000351 0.996224 0.031766 0.000001 0.007581 NaN
15 0.110243 0.179627 0.999797 0.787808 0.061839 0.094904 NaN
19 0.013756 0.589113 0.419284 0.974934 0.000009 0.005824 NaN
25 0.337045 0.012857 0.124466 0.005812 0.008698 0.028932 NaN
29 0.528113 0.704117 0.201008 0.101751 0.691210 0.717842 NaN

In [77]:
ps2=ps.copy()
ps3=ps2>=thr
ps3.sum()/float(len(ps3.index))


Out[77]:
FEF-left_roi.mat     0.669492
FEF-right_roi.mat    0.567797
MT-left_roi.mat      0.838983
MT-right_roi.mat     0.627119
SP-Left_roi.mat      0.457627
SP-Right_roi.mat     0.584746
Unnamed: 7           0.000000
dtype: float64

In [78]:
ps2[ps3]=float("nan")
ps2.head()


Out[78]:
FEF-left_roi.mat FEF-right_roi.mat MT-left_roi.mat MT-right_roi.mat SP-Left_roi.mat SP-Right_roi.mat Unnamed: 7
subject
9 0.000047 0.000351 NaN 0.031766 0.000001 0.007581 NaN
15 NaN NaN NaN NaN NaN NaN NaN
19 0.013756 NaN NaN NaN 0.000009 0.005824 NaN
25 NaN 0.012857 NaN 0.005812 0.008698 0.028932 NaN
29 NaN NaN NaN NaN NaN NaN NaN

In [79]:
ef2=ef.copy()
ef2[ps3]=float("nan")
ef2.head()


Out[79]:
FEF-left_roi.mat FEF-right_roi.mat MT-left_roi.mat MT-right_roi.mat SP-Left_roi.mat SP-Right_roi.mat Unnamed: 7
subject
9 2.202807 1.870303 NaN 1.233423 2.787738 2.194587 NaN
15 NaN NaN NaN NaN NaN NaN NaN
19 1.063118 NaN NaN NaN 1.721303 1.666498 NaN
25 NaN 1.688934 NaN 2.502616 1.809077 2.404814 NaN
29 NaN NaN NaN NaN NaN NaN NaN

In [80]:
#remove last column
ef3=ef2[ef2.columns[:-1]]
col_names = [c.split(filter_name_char)[0] for c in ef3.columns]
ef3.columns=col_names
ef3.head()


Out[80]:
FEF-left FEF-right MT-left MT-right SP-Left SP-Right
subject
9 2.202807 1.870303 NaN 1.233423 2.787738 2.194587
15 NaN NaN NaN NaN NaN NaN
19 1.063118 NaN NaN NaN 1.721303 1.666498
25 NaN 1.688934 NaN 2.502616 1.809077 2.404814
29 NaN NaN NaN NaN NaN NaN

In [81]:
ef3.to_excel(excel_file,merge_cells=False)

In [82]:
ps4=ps2[ps2.columns[:-1]]
ps4.columns=col_names
ps4.to_excel(excel_file_p,merge_cells=False)

In [82]: