def seed_connectivity_overview(
template='/usr/share/mouse-brain-atlases/dsurqec_40micron_masked.nii',
cut_coords=[None, [0, -4.9, -3.3]],
plot=False,
):
import numpy as np
from labbookdb.report.tracking import treatment_group, append_external_identifiers
from samri.plotting.overview import multiplot_matrix, multipage_plot
db_path = '~/syncdata/meta.db'
groups = treatment_group(db_path, ['cFluDW', 'cFluDW_'], 'cage')
groups = append_external_identifiers(db_path, groups, ['Genotype_code'])
all_subjects = groups['ETH/AIC'].unique()
treatment = groups[(groups['Genotype_code'] == "eptg")
& (groups['Cage_TreatmentProtocol_code'] == "cFluDW"
)]['ETH/AIC'].tolist()
no_treatment = groups[(groups['Genotype_code'] == "eptg")
& (groups['Cage_TreatmentProtocol_code'] == "cFluDW_"
)]['ETH/AIC'].tolist()
negative_controls = groups[groups['Genotype_code'] ==
"epwt"]['ETH/AIC'].tolist()
print(treatment, no_treatment, negative_controls)
substitutions = bids_substitution_iterator(
["ofM", "ofM_aF", "ofM_cF1", "ofM_cF2", "ofM_pF"],
all_subjects,
[
"EPI_CBV_chr_longSOA",
],
"~/ni_data/ofM.dr/",
"composite",
)
fc_results = fc.seed_based(
substitutions,
"~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz",
'/usr/share/mouse-brain-atlases/dsurqec_200micron_mask.nii',
ts_file_template=
"~/ni_data/ofM.dr/preprocessing/{preprocessing_dir}/sub-{subject}/ses-{session}/func/sub-{subject}_ses-{session}_task-{task}.nii.gz",
)
print([i['subject'] for i in fc_results])
return
if plot:
multipage_plot(
fc_results,
treatment,
figure_title="Chronic Fluoxetine (drinking water) Treatment Group",
template=template,
threshold=0.1,
base_cut_coords=cut_coords,
save_as="fc_treatment.pdf",
overlays=['~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz'],
scale=0.4,
)
multipage_plot(
fc_results,
no_treatment,
figure_title="Chronic Fluoxetine (drinking water) Treatment Group",
template=template,
threshold=0.1,
base_cut_coords=cut_coords,
save_as="fc_no_treatment.pdf",
overlays=['~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz'],
scale=0.4,
)
multipage_plot(
fc_results,
negative_controls,
figure_title="Chronic Fluoxetine (drinking water) Treatment Group",
template=template,
threshold=0.1,
base_cut_coords=cut_coords,
save_as="fc_negative_control.pdf",
overlays=['~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz'],
scale=0.4,
)
def seed_connectivity_overview(
template="/usr/share/mouse-brain-atlases/dsurqec_40micron_masked.nii",
cut_coords=[None,[0,-4.9,-3.3]],
plot=False,
):
import numpy as np
from labbookdb.report.tracking import treatment_group, append_external_identifiers
from samri.plotting.overview import multiplot_matrix, multipage_plot
db_path = '~/syncdata/meta.db'
groups = treatment_group(db_path, ['cFluDW','cFluDW_'], 'cage')
groups = append_external_identifiers(db_path, groups, ['Genotype_code'])
all_subjects = groups['ETH/AIC'].unique()
treatment = groups[
(groups['Genotype_code']=="eptg")&
(groups['Cage_TreatmentProtocol_code']=="cFluDW")
]['ETH/AIC'].tolist()
no_treatment = groups[
(groups['Genotype_code']=="eptg")&
(groups['Cage_TreatmentProtocol_code']=="cFluDW_")
]['ETH/AIC'].tolist()
negative_controls = groups[groups['Genotype_code']=="epwt"]['ETH/AIC'].tolist()
print(treatment, no_treatment, negative_controls)
substitutions = bids_substitution_iterator(
["ofM","ofM_aF","ofM_cF1","ofM_cF2","ofM_pF"],
all_subjects,
["EPI_CBV_chr_longSOA",],
"~/ni_data/ofM.dr/",
"composite",
)
fc_results = fc.seed_based(substitutions, "~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz", "/usr/share/mouse-brain-atlases/dsurqec_200micron_mask.nii",
ts_file_template="~/ni_data/ofM.dr/preprocessing/{preprocessing_dir}/sub-{subject}/ses-{session}/func/sub-{subject}_ses-{session}_task-{task}.nii.gz",
)
print([i['subject'] for i in fc_results])
return
if plot:
multipage_plot(fc_results, treatment,
figure_title="Chronic Fluoxetine (drinking water) Treatment Group",
template=template,
threshold=0.1,
base_cut_coords=cut_coords,
save_as="fc_treatment.pdf",
overlays=['~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz'],
scale=0.4,
)
multipage_plot(fc_results, no_treatment,
figure_title="Chronic Fluoxetine (drinking water) Treatment Group",
template=template,
threshold=0.1,
base_cut_coords=cut_coords,
save_as="fc_no_treatment.pdf",
overlays=['~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz'],
scale=0.4,
)
multipage_plot(fc_results, negative_controls,
figure_title="Chronic Fluoxetine (drinking water) Treatment Group",
template=template,
threshold=0.1,
base_cut_coords=cut_coords,
save_as="fc_negative_control.pdf",
overlays=['~/ni_data/templates/roi/DSURQEc_dr_xs.nii.gz'],
scale=0.4,
)