def p_clusters(mask):
from samri.plotting import summary, timeseries
substitutions = bids_substitution_iterator(
["ofM", "ofM_aF", "ofM_cF1", "ofM_cF2", "ofM_pF"],
["4011", "4012", "5689", "5690", "5691"],
# ["4007","4008","4011","4012","5689","5690","5691"],
["EPI_CBV_jb_long", "EPI_CBV_chr_longSOA"],
"~/ni_data/ofM.dr/",
"composite",
l1_dir="dr",
)
timecourses, designs, stat_maps, events_dfs, subplot_titles = summary.p_filtered_ts(
substitutions,
ts_file_template=
"{data_dir}/preprocessing/{preprocessing_dir}/sub-{subject}/ses-{session}/func/sub-{subject}_ses-{session}_task-{scan}.nii.gz",
beta_file_template=
"{data_dir}/l1/{l1_dir}/sub-{subject}/ses-{session}/sub-{subject}_ses-{session}_task-{scan}_cope.nii.gz",
# p_file_template="{data_dir}/l1/{l1_dir}/sub-{subject}/ses-{session}/sub-{subject}_ses-{session}_task-{scan}_pstat.nii.gz",
p_file_template=
"{data_dir}/l1/{l1_dir}/sub-{subject}/ses-{session}/sub-{subject}_ses-{session}_task-{scan}_pfstat.nii.gz",
design_file_template=
"{data_dir}/l1/{l1_workdir}/_subject_session_scan_{subject}.{session}.{scan}/modelgen/run0.mat",
event_file_template=
"{data_dir}/preprocessing/{preprocessing_dir}/sub-{subject}/ses-{session}/func/sub-{subject}_ses-{session}_task-{scan}_events.tsv",
brain_mask=mask,
p_level=0.05,
)
timeseries.multi(timecourses,
designs,
stat_maps,
events_dfs,
subplot_titles,
figure="timecourses")
plt.show()
def p_clusters(mask):
from samri.plotting import summary, timeseries
substitutions = bids_substitution_iterator(
["ofM","ofM_aF","ofM_cF1","ofM_cF2","ofM_pF"],
["4011","4012","5689","5690","5691"],
# ["4007","4008","4011","4012","5689","5690","5691"],
["EPI_CBV_jb_long","EPI_CBV_chr_longSOA"],
"~/ni_data/ofM.dr/",
"composite",
l1_dir="dr",
)
timecourses, designs, stat_maps, events_dfs, subplot_titles = summary.p_filtered_ts(substitutions,
ts_file_template="{data_dir}/preprocessing/{preprocessing_dir}/sub-{subject}/ses-{session}/func/sub-{subject}_ses-{session}_task-{scan}.nii.gz",
beta_file_template="{data_dir}/l1/{l1_dir}/sub-{subject}/ses-{session}/sub-{subject}_ses-{session}_task-{scan}_cope.nii.gz",
# p_file_template="{data_dir}/l1/{l1_dir}/sub-{subject}/ses-{session}/sub-{subject}_ses-{session}_task-{scan}_pstat.nii.gz",
p_file_template="{data_dir}/l1/{l1_dir}/sub-{subject}/ses-{session}/sub-{subject}_ses-{session}_task-{scan}_pfstat.nii.gz",
design_file_template="{data_dir}/l1/{l1_workdir}/_subject_session_scan_{subject}.{session}.{scan}/modelgen/run0.mat",
event_file_template="{data_dir}/preprocessing/{preprocessing_dir}/sub-{subject}/ses-{session}/func/sub-{subject}_ses-{session}_task-{scan}_events.tsv",
brain_mask=mask,
p_level=0.05,
)
timeseries.multi(timecourses, designs, stat_maps, events_dfs, subplot_titles, figure="timecourses")
plt.show()