add SEM error in plot and threshold to find modulation
new script