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()