Does the brain's E:I balance really shape long-range temporal correlations? Lessons learned from 3T MRI.
The 9 matches · 2 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
- [1] § Methods › MRS ↔ EI_Hurst_Measure_Collection_Expanded.sh, lines 88–135 · score 0.74 · Asc, Asp, GPC, NAAG, GSH, Lac
- [2] § Methods › MRS ↔ 06_MeasureCollection.sh, lines 188–268 · score 0.73 · Asc, Asp, GPC, NAAG, GSH, Lac
- [3] § Methods › MRS ↔ gannetauto.m, the whole file · a weak match · score 0.67 · water suppressed, MEGA PRESS, Philips, co, metabolite, sequences
- [4] § Methods › Hurst exponent calculation ↔ dfa/notebook/dfa_simulations.ipynb, lines 16–105 · score 0.65 · power law, Hurst exponent, windows, algorithms, log, density
- [5] § Methods › MRS ↔ 06_MeasureCollection.sh, lines 44–114 · score 0.64 · water scaling, LCModel, sLASER, quantified, spectra, Glu
- [6] § Methods › Image processing › fMRI ↔ 01_BashPipeline.sh, lines 180–218 · score 0.57 · preprocessed BOLD, fMRIPrep, Tedana, echo
- [7] § Methods › Acquisition details ↔ gannetauto.m, the whole file · a weak match · score 0.54 · water suppression, MEGA PRESS, weighted, sequence, scans, voxel
- [8] § Methods › Hurst exponent calculation ↔ dfa/dfa.py, lines 39–90 · score 0.52 · power law, windows, algorithms, log, density, exponent
- [9] § Methods › Image processing › Structural images ↔ 01_BashPipeline.sh, lines 540–604 · score 0.51 · ANTs, FSL, v2, GM, spaces, mask
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
Shell · 371 lines · 38 KB · no license · 2 matches
- #!/bin/bash
- set -e # exit when any command fails
- BIDSdir=/mnt/WeberLab/Projects/EI_Hurst
- codedir=${BIDSdir}/code
- subjectfile=${codedir}/SubjectDCMInfo.csv
- derivatives=${BIDSdir}/derivatives
- MEGAPRESSdir=${derivatives}/megapress
- SLASERdir=${derivatives}/slaser
- hurstdir=${derivatives}/hurst
- alffdir=${derivatives}/ALFF
- mriqcdir=${derivatives}/mriqc
- tedanadir=${derivatives}/tedana
- workdir=$HOME/work
- fslicense=/mnt/WeberLab
- sed 1d $subjectfile | grep -v '^#' | while read d # ignore first row, ignore rows starting with '#'
- do
- subid=$(echo $d | awk -F',' '{print $1}')
- echo "#------------------------------------------------------------------#"
- echo "Starting sub-${subid}"
- ## Collecing Osprey sLASER data
- echo "#------------------------------------------------------------------#"
- echo "Step 1: Collecting SLASER measures"
- echo "Step 1a: Collecting SLASER Values"
- rest_vals=$(awk -F '\t' '(NR>1) {print}' ${SLASERdir}/sub-${subid}/rest/MarkBasisSet/QuantifyResults/A_TissCorrWaterScaled_Voxel_1_Basis_1.tsv)
- rest_glu_tcr=$(awk -F '\t' '(NR>1) {print $9}' ${SLASERdir}/sub-${subid}/rest/MarkBasisSet/QuantifyResults/A_tCr_Voxel_1_Basis_1.tsv)
- movie_vals=$(awk -F '\t' '(NR>1) {print}' ${SLASERdir}/sub-${subid}/movie/MarkBasisSet/QuantifyResults/A_TissCorrWaterScaled_Voxel_1_Basis_1.tsv)
- movie_glu_tcr=$(awk -F '\t' '(NR>1) {print $9}' ${SLASERdir}/sub-${subid}/movie/MarkBasisSet/QuantifyResults/A_tCr_Voxel_1_Basis_1.tsv)
- rest_vals_old=$(awk -F '\t' '(NR>1) {print}' ${SLASERdir}/sub-${subid}/rest/QuantifyResults/A_TissCorrWaterScaled_Voxel_1_Basis_1.tsv)
- movie_vals_old=$(awk -F '\t' '(NR>1) {print}' ${SLASERdir}/sub-${subid}/movie/QuantifyResults/A_TissCorrWaterScaled_Voxel_1_Basis_1.tsv)
- rest_glx_tcr=$(awk -F '\t' '(NR>1) {print $28}' ${SLASERdir}/sub-${subid}/rest/MarkBasisSet/QuantifyResults/A_tCr_Voxel_1_Basis_1.tsv)
- movie_glx_tcr=$(awk -F '\t' '(NR>1) {print $28}' ${SLASERdir}/sub-${subid}/movie/MarkBasisSet/QuantifyResults/A_tCr_Voxel_1_Basis_1.tsv)
- echo "Step 1b: Collecting FWHM Values"
- rest_fwhm_slaser=$(awk -F '\t' '(NR>1) {print $3}' ${SLASERdir}/sub-${subid}/rest/MarkBasisSet/QM_processed_spectra.tsv)
- movie_fwhm_slaser=$(awk -F '\t' '(NR>1) {print $3}' ${SLASERdir}/sub-${subid}/movie/MarkBasisSet/QM_processed_spectra.tsv)
- rest_fwhm_slaser_old=$(awk -F '\t' '(NR>1) {print $3}' ${SLASERdir}/sub-${subid}/rest/QM_processed_spectra.tsv)
- movie_fwhm_slaser_old=$(awk -F '\t' '(NR>1) {print $3}' ${SLASERdir}/sub-${subid}/movie/QM_processed_spectra.tsv)
- echo "Step 1: Complete"
- ## Collecting LCModel sLASER
- echo "#------------------------------------------------------------------#"
- echo "Step 2: Collecting LCModel sLASER measures"
- rest_lc_vals=$(awk -F ',' '(NR>1) {print $3,$4,$6,$7,$9,$10,$12,$13,$15,$16,$18,$19,$21,$22,$24,$25,$27,$28,$30,$31,$33,$34,$36,$37,$39,$40,$42,$43,$45,$46,$48,$49,$51,$52,$54,$55,$57, \
- $58,$60,$61,$63,$64,$66,$67,$69,$70,$72,$73,$75,$76,$78,$79,$81,$82,$84,$85,$87,$88}' ${SLASERdir}/sub-${subid}/rest/MarkBasisSet/LCModelFiles/EI_Hurst_sub-${subid}_sLASER_LCModel_Mark.csv)
- movie_lc_vals=$(awk -F ',' '(NR>1) {print $3,$4,$6,$7,$9,$10,$12,$13,$15,$16,$18,$19,$21,$22,$24,$25,$27,$28,$30,$31,$33,$34,$36,$37,$39,$40,$42,$43,$45,$46,$48,$49,$51,$52,$54,$55,$57, \
- $58,$60,$61,$63,$64,$66,$67,$69,$70,$72,$73,$75,$76,$78,$79,$81,$82,$84,$85,$87,$88}' ${SLASERdir}/sub-${subid}/movie/MarkBasisSet/LCModelFiles/EI_Hurst_sub-${subid}_sLASER_LCModel.csv)
- ##Old Basis Set
- rest_lc_vals_old=$(awk -F ',' '(NR>1) {print $3,$4,$6,$7,$9,$10,$12,$13,$15,$16,$18,$19,$21,$22,$24,$25,$27,$28,$30,$31,$33,$34,$36,$37,$39,$40,$42,$43,$45,$46,$48,$49,$51,$52,$54,$55,$57, \
- $58,$60,$61,$63,$64,$66,$67,$69,$70,$72,$73,$75,$76,$78,$79,$81,$82,$84,$85,$87,$88}' ${SLASERdir}/sub-${subid}/rest/LCModelFiles/EI_Hurst_sub-${subid}_sLASER_LCModel.csv)
- movie_lc_vals_old=$(awk -F ',' '(NR>1) {print $3,$4,$6,$7,$9,$10,$12,$13,$15,$16,$18,$19,$21,$22,$24,$25,$27,$28,$30,$31,$33,$34,$36,$37,$39,$40,$42,$43,$45,$46,$48,$49,$51,$52,$54,$55,$57, \
- $58,$60,$61,$63,$64,$66,$67,$69,$70,$72,$73,$75,$76,$78,$79,$81,$82,$84,$85,$87,$88}' ${SLASERdir}/sub-${subid}/movie/LCModelFiles/EI_Hurst_sub-${subid}_sLASER_LCModel.csv)
- echo "Step 2: Complete"
- ## Collecting MEGAPRESS data
- echo "#------------------------------------------------------------------#"
- echo "Step 3: Collecting GABA measures -- Osprey"
- echo "Step 3a: Collecting GABA + Glu Values"
- rest_gaba_gabaplus=$(awk -F '\t' '(NR>1) {print $1, $2}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/MarkBasisSet/QuantifyResults/diff1_AlphaCorrWaterScaled_Voxel_1_Basis_1.tsv)
- movie_gaba_gabaplus=$(awk -F '\t' '(NR>1) {print $1, $2}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/MarkBasisSet/QuantifyResults/diff1_AlphaCorrWaterScaled_Voxel_1_Basis_1.tsv)
- rest_glu_megapress=$(awk -F '\t' '(NR>1) {print $9}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/MarkBasisSet/QuantifyResults/A_TissCorrWaterScaled_Voxel_1_Basis_1.tsv)
- movie_glu_megapress=$(awk -F '\t' '(NR>1) {print $9}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/MarkBasisSet/QuantifyResults/A_TissCorrWaterScaled_Voxel_1_Basis_1.tsv)
- rest_gabaplus_tcr=$(awk -F '\t' '(NR>1) {print $31}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/MarkBasisSet/QuantifyResults/diff1_tCr_Voxel_1_Basis_1.tsv)
- movie_gabaplus_tcr=$(awk -F '\t' '(NR>1) {print $31}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/MarkBasisSet/QuantifyResults/diff1_tCr_Voxel_1_Basis_1.tsv)
- rest_gaba_tcr=$(awk -F '\t' '(NR>1) {print $5}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/MarkBasisSet/QuantifyResults/diff1_tCr_Voxel_1_Basis_1.tsv)
- movie_gaba_tcr=$(awk -F '\t' '(NR>1) {print $5}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/MarkBasisSet/QuantifyResults/diff1_tCr_Voxel_1_Basis_1.tsv)
- echo "Step 3b: Collecting FWHM Values"
- rest_fwhm_osprey=$(awk -F '\t' '(NR>1) {print $3}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/MarkBasisSet/QM_processed_spectra.tsv)
- movie_fwhm_osprey=$(awk -F '\t' '(NR>1) {print $3}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/MarkBasisSet/QM_processed_spectra.tsv)
- echo "Step 3c: Collecting GABA measures -- Gannet"
- gannet_rest=$(find ${MEGAPRESSdir}/sub-${subid}/rest_gannet/MarkBasisSet/GannetQuantify_output/ -type f -name '*.pdf')
- pdftotext ${gannet_rest}
- txt_gannet_rest=$(find ${MEGAPRESSdir}/sub-${subid}/rest_gannet/MarkBasisSet/GannetQuantify_output/ -type f -name '*.txt')
- rest_gaba_gannet_temp=$(awk '(NR==26) {print}' ${txt_gannet_rest} | grep -Eo '[0-9]{1,4}')
- rest_gaba_gannet=$(echo $rest_gaba_gannet_temp | tr " " .)
- rest_glx_gannet_temp=$(awk '(NR==32) {print}' ${txt_gannet_rest} | grep -Eo '[0-9]{1,4}')
- rest_glx_gannet=$(echo $rest_glx_gannet_temp | tr " " .)
- gannet_rest_tcr=$(find ${MEGAPRESSdir}/sub-${subid}/rest_gannet/MarkBasisSet/GannetFit_output/ -type f -name '*.pdf')
- pdftotext ${gannet_rest_tcr}
- txt_gannet_rest_tcr=$(find ${MEGAPRESSdir}/sub-${subid}/rest_gannet/MarkBasisSet/GannetFit_output/ -type f -name '*.txt')
- rest_gaba_tcr_gannet=$(sed -n '106p' < ${txt_gannet_rest_tcr})
- gannet_movie=$(find ${MEGAPRESSdir}/sub-${subid}/movie_gannet/MarkBasisSet/GannetQuantify_output/ -type f -name '*.pdf')
- pdftotext ${gannet_movie}
- txt_gannet_movie=$(find ${MEGAPRESSdir}/sub-${subid}/movie_gannet/MarkBasisSet/GannetQuantify_output/ -type f -name '*.txt')
- movie_gaba_gannet_temp=$(awk '(NR==26) {print}' ${txt_gannet_movie} | grep -Eo '[0-9]{1,4}')
- movie_gaba_gannet=$(echo $movie_gaba_gannet_temp | tr " " .)
- movie_glx_gannet_temp=$(awk '(NR==32) {print}' ${txt_gannet_movie} | grep -Eo '[0-9]{1,4}')
- movie_glx_gannet=$(echo $movie_glx_gannet_temp | tr " " .)
- gannet_movie_tcr=$(find ${MEGAPRESSdir}/sub-${subid}/movie_gannet/MarkBasisSet/GannetFit_output/ -type f -name '*.pdf')
- pdftotext ${gannet_movie_tcr}
- txt_gannet_movie_tcr=$(find ${MEGAPRESSdir}/sub-${subid}/movie_gannet/MarkBasisSet/GannetFit_output/ -type f -name '*.txt')
- movie_gaba_tcr_gannet=$(sed -n '106p' < ${txt_gannet_movie_tcr})
- ## Old Basis Set
- echo "Step 3d: Collecting GABA Values -- Old Basis Set"
- rest_gaba_gabaplus_old=$(awk -F '\t' '(NR>1) {print $1, $2}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/QuantifyResults/diff1_AlphaCorrWaterScaled_Voxel_1_Basis_1.tsv)
- movie_gaba_gabaplus_old=$(awk -F '\t' '(NR>1) {print $1, $2}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/QuantifyResults/diff1_AlphaCorrWaterScaled_Voxel_1_Basis_1.tsv)
- echo "Step 3e: Collecting FWHM Values -- Old Basis Set"
- rest_fwhm_osprey_old=$(awk -F '\t' '(NR>1) {print $3}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/QM_processed_spectra.tsv)
- movie_fwhm_osprey_old=$(awk -F '\t' '(NR>1) {print $3}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/QM_processed_spectra.tsv)
- echo "Step 3f: Collecting GABA measures -- Gannet; Old Basis Set"
- gannet_rest_old=$(find ${MEGAPRESSdir}/sub-${subid}/rest_gannet/GannetQuantify_output/ -type f -name '*.pdf')
- pdftotext ${gannet_rest_old}
- txt_gannet_rest_old=$(find ${MEGAPRESSdir}/sub-${subid}/rest_gannet/GannetQuantify_output/ -type f -name '*.txt')
- rest_gaba_gannet_temp_old=$(awk '(NR==26) {print}' ${txt_gannet_rest_old} | grep -Eo '[0-9]{1,4}')
- rest_gaba_gannet_old=$(echo $rest_gaba_gannet_temp_old | tr " " .)
- gannet_movie_old=$(find ${MEGAPRESSdir}/sub-${subid}/movie_gannet/GannetQuantify_output/ -type f -name '*.pdf')
- pdftotext ${gannet_movie_old}
- txt_gannet_movie_old=$(find ${MEGAPRESSdir}/sub-${subid}/movie_gannet/GannetQuantify_output/ -type f -name '*.txt')
- movie_gaba_gannet_temp_old=$(awk '(NR==26) {print}' ${txt_gannet_movie_old} | grep -Eo '[0-9]{1,4}')
- movie_gaba_gannet_old=$(echo $movie_gaba_gannet_temp_old | tr " " .)
- echo "Step 3: Complete"
- ## Collecting LCModel MEGAPRESS data
- echo "#------------------------------------------------------------------#"
- echo "Step 4: Collecting LCModel GABA measures"
- rest_lc_gaba=$(awk -F ',' '(NR>1) {print $15,$16}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/MarkBasisSet/LCModelFiles/EI_Hurst_sub-${subid}_MEGAPRESS68_LCModel_DIFF1_Mark.csv)
- movie_lc_gaba=$(awk -F ',' '(NR>1) {print $15,$16}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/MarkBasisSet/LCModelFiles/EI_Hurst_sub-${subid}_MEGAPRESS68_LCModel_DIFF1_Mark.csv)
- rest_lc_gaba_old=$(awk -F ',' '(NR>1) {print $15,$16}' ${MEGAPRESSdir}/sub-${subid}/rest_osprey/LCModelFiles/EI_Hurst_sub-${subid}_MEGAPRESS68_LCModel_DIFF1.csv)
- movie_lc_gaba_old=$(awk -F ',' '(NR>1) {print $15,$16}' ${MEGAPRESSdir}/sub-${subid}/movie_osprey/LCModelFiles/EI_Hurst_sub-${subid}_MEGAPRESS68_LCModel_DIFF1.csv)
- echo "Step 4: Complete"
- ## Collecting Hurst data
- echo "#------------------------------------------------------------------#"
- echo "Step 5: Collecting Hurst measures"
- rest_hurst_whole=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.101.nii.gz -n -M`
- rest_hurst_fullfreq_whole=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -M`
- rest_hurst_fullfreq=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_hurst_fullfreq_blur=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq_7mmblur.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_hurst=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.151.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_hurst_old=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.101.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_hurst_vold=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.088.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- movie_hurst_whole=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.101.nii.gz -n -M`
- movie_hurst_fullfreq_whole=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -M`
- movie_hurst_fullfreq=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_hurst_fullfreq_blur=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq_7mmblur.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_hurst=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.151.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_hurst_old=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.101.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_hurst_vold=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.088.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- rest_hurst_vn=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.151.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_hurst_fullfreq_vn=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_hurst_fullfreq_blur_vn=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq_7mmblur.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_hurst_vn_old=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.101.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_hurst_vn_vold=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.088.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- movie_hurst_vn=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.151.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_hurst_fullfreq_vn=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_hurst_fullfreq_blur_vn=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq_7mmblur.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_hurst_vn_old=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.101.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_hurst_vn_vold=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_0.013_to_0.088.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- rest_hurst_dmn=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_DMN_to_native.nii.gz -M`
- movie_hurst_dmn=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_DMN_to_native.nii.gz -M`
- rest_hurst_fpn=`fslstats ${hurstdir}/sub-${subid}/rest/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_FPN_to_native.nii.gz -M`
- movie_hurst_fpn=`fslstats ${hurstdir}/sub-${subid}/movie/desc-optcomDenoised_bold_WelchHurst_fullfreq.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_FPN_to_native.nii.gz -M`
- echo "Step 5: Complete"
- ## Collecting BOLD data
- echo "#------------------------------------------------------------------#"
- echo "Step 6: Collecting BOLD data"
- rest_bold=`fslstats ${tedanadir}/sub-${subid}/rest/desc-optcomDenoised_bold_Tmean.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_bold_vn=`fslstats ${tedanadir}/sub-${subid}/rest/desc-optcomDenoised_bold_Tmean.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- movie_bold=`fslstats ${tedanadir}/sub-${subid}/movie/desc-optcomDenoised_bold_Tmean.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_bold_vn=`fslstats ${tedanadir}/sub-${subid}/movie/desc-optcomDenoised_bold_Tmean.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- echo "Step 6: Complete"
- ## Collecting ALFF data
- echo "#------------------------------------------------------------------#"
- echo "Step 7: Collecting ALFF"
- rest_alff=`fslstats ${alffdir}/sub-${subid}/rest/alff.nii -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- movie_alff=`fslstats ${alffdir}/sub-${subid}/movie/alff.nii -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- rest_alff_vn=`fslstats ${alffdir}/sub-${subid}/rest/alff.nii -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- movie_alff_vn=`fslstats ${alffdir}/sub-${subid}/movie/alff.nii -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- rest_alff_blur=`fslstats ${alffdir}/sub-${subid}/rest/alff_7mmblur.nii -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- movie_alff_blur=`fslstats ${alffdir}/sub-${subid}/movie/alff_7mmblur.nii -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- rest_alff_blur_vn=`fslstats ${alffdir}/sub-${subid}/rest/alff_7mmblur.nii -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- movie_alff_blur_vn=`fslstats ${alffdir}/sub-${subid}/movie/alff_7mmblur.nii -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- echo "Step 7: Complete"
- ## Collecting MeanFD data
- echo "#------------------------------------------------------------------#"
- echo "Step 8: Collecting MeanFD"
- rest_meanfd=`jq .fd_mean ${mriqcdir}/sub-${subid}/func/sub-${subid}_task-rest_acq-rest_echo-1_bold.json`
- movie_meanfd=`jq .fd_mean ${mriqcdir}/sub-${subid}/func/sub-${subid}_task-rest_acq-movie_echo-1_bold.json`
- echo "Step 8: Complete"
- ## Parsing TSV file
- echo "#------------------------------------------------------------------#"
- echo "Step 9a: Parsing TSV"
- participants_measures=${BIDSdir}/participant_measures
- [ -d $participants_measures ] || mkdir -p $participants_measures
- echo "Rest_Glx_tCr Movie_Glx_tCr Rest_GABAplus_tCr Movie_GABAplus_tCr Rest_GABA_tCr Movie_GABA_tCr Rest_Glu_tCR Movie_Glu_tCr Rest_Glu_Megapress Rest_Asc Rest_Asp Rest_Cr Rest_CrCH2 Rest_GABA_slaser Rest_GPC Rest_GSH Rest_Gln Rest_Glu \
- Rest_mI Rest_Lac Rest_NAA Rest_NAAG Rest_PCh Rest_PCr Rest_PE Rest_sI Rest_Tau Rest_MM09 Rest_MM12 Rest_MM14 \
- Rest_MM17 Rest_MM20 Rest_Lip09 Rest_Lip13 Rest_Lip20 Rest_tNAA Rest_Glx Rest_tCho Rest_tCr Rest_FWHM_SLASER \
- Movie_Asc Movie_Glu_Megapress Movie_Asp Movie_Cr Movie_CrCH2 Movie_GABA_slaser Movie_GPC Movie_GSH Movie_Gln \
- Movie_Glu Movie_mI Movie_Lac Movie_NAA Movie_NAAG Movie_PCh Movie_PCr Movie_PE Movie_sI Movie_Tau \
- Movie_MM09 Movie_MM12 Movie_MM14 Movie_MM17 Movie_MM20 Movie_Lip09 Movie_Lip13 Movie_Lip20 Movie_tNAA \
- Movie_Glx Movie_tCho Movie_tCr Movie_FWHM_SLASER Rest_GABA Rest_GABAplus Rest_GABA_Gannet Rest_FWHM_Osprey \
- Movie_GABA Movie_GABAplus Movie_GABA_Gannet Movie_FWHM_Osprey Rest_Hurst Rest_Hurst_Old Rest_Hurst_VOld Rest_Hurst_VN \
- Rest_Hurst_FullFreq_VN Rest_Hurst_FullFreq_Blur_VN Rest_Hurst_VN_Old Rest_Hurst_VN_VOld Movie_Hurst Movie_Hurst_Old Movie_Hurst_VOld Movie_Hurst_VN \
- Movie_Hurst_FullFreq_VN Movie_Hurst_FullFreq_Blur_VN Movie_Hurst_VN_Old Movie_Hurst_VN_VOld Rest_MeanFD \
- Movie_MeanFD Rest_ALFF Rest_ALFF_Blur Rest_ALFF_VN Rest_ALFF_Blur_VN Movie_ALFF Movie_ALFF_Blur Movie_ALFF_VN Movie_ALFF_Blur_VN LC_Rest_Ala LC_Rest_AlaSD LC_Rest_Asc LC_Rest_AscSD \
- LC_Rest_Asp LC_Rest_AspSD LC_Rest_Cr LC_Rest_CrSD LC_Rest_GABA_sl LC_Rest_GABASD_sl LC_Rest_GSH LC_Rest_GSHSD \
- LC_Rest_Gln LC_Rest_GlnSD \
- LC_Rest_Glu LC_Rest_GluSD LC_Rest_Lac LC_Rest_LacSD LC_Rest_NAA LC_Rest_NAASD LC_Rest_NAAG LC_Rest_NAAGSD \
- LC_Rest_PCr LC_Rest_PCrSD LC_Rest_PE LC_Rest_PESD LC_Rest_Tau LC_Rest_TauSD LC_Rest_mI LC_Rest_mISD \
- LC_Rest_sI LC_Rest_sISD LC_Rest_CrCH2 LC_Rest_CrCH2SD LC_Rest_CrPCr LC_Rest_CrPrSD LC_Rest_NAANAAG \
- LC_Rest_NAANAAGSD LC_Rest_GluGln LC_Rest_GluGlnSD LC_Rest_Lip13a LC_Rest_Lip13aSD LC_Rest_Lip13b LC_Rest_Lip13bSD \
- LC_Rest_Lip09 LC_Rest_Lip09SD LC_Rest_MM09 LC_Rest_MM09SD LC_Rest_Lip20 LC_Rest_Lip20SD LC_Rest_MM20 \
- LC_Rest_MM20SD LC_Rest_MM12 LC_Rest_MM12SD LC_Rest_MM14 LC_Rest_MM14SD LC_Rest_MM17 LC_Rest_MM17SD \
- LC_Movie_Ala LC_Movie_AlaSD LC_Movie_Asc LC_Movie_AscSD LC_Movie_Asp LC_Movie_AspSD LC_Movie_Cr LC_Movie_CrSD \
- LC_Movie_GABA_sl LC_Movie_GABASD_sl LC_Movie_GSH LC_Movie_GSHSD LC_Movie_Gln LC_Movie_GlnSD \
- LC_Movie_Glu LC_Movie_GluSD LC_Movie_Lac LC_Movie_LacSD LC_Movie_NAA LC_Movie_NAASD LC_Movie_NAAG \
- LC_Movie_NAAGSD LC_Movie_PCr LC_Movie_PCrSD LC_Movie_PE LC_Movie_PESD LC_Movie_Tau LC_Movie_TauSD \
- LC_Movie_mI LC_Movie_mISD \
- LC_Movie_sI LC_Movie_sISD LC_Movie_CrCH2 LC_Movie_CrCH2SD LC_Movie_CrPCr LC_Movie_CrPrSD LC_Movie_NAANAAG \
- LC_Movie_NAANAAGSD LC_Movie_GluGln LC_Movie_GluGlnSD LC_Movie_Lip13a LC_Movie_Lip13aSD LC_Movie_Lip13b \
- LC_Movie_Lip13bSD \
- LC_Movie_Lip09 LC_Movie_Lip09SD LC_Movie_MM09 LC_Movie_MM09SD LC_Movie_Lip20 LC_Movie_Lip20SD LC_Movie_MM20 \
- LC_Movie_MM20SD LC_Movie_MM12 LC_Movie_MM12SD LC_Movie_MM14 LC_Movie_MM14SD LC_Movie_MM17 LC_Movie_MM17SD \
- LC_Rest_GABA \
- LC_Rest_GABASD LC_Movie_GABA LC_Movie_GABASD Rest_BOLD Rest_BOLD_VN Movie_BOLD Movie_BOLD_VN Rest_Hurst_Whole \
- Rest_Hurst_FullFreq Rest_Hurst_FullFreq_Blur Movie_Hurst_Whole Movie_Hurst_FullFreq Movie_Hurst_FullFreq_Blur \
- Rest_Hurst_FullFreq_Whole Movie_Hurst_FullFreq_Whole Rest_Glx_Gannet Movie_Glx_Gannet Rest_Hurst_DMN Movie_Hurst_DMN \
- Rest_Hurst_FPN Movie_Hurst_FPN Rest_GABA_Gannet_tCr Movie_GABA_Gannet_tCr" > ${participants_measures}/sub-${subid}_EI_Hurst_Measures_Expanded.tsv
- echo $rest_glx_tcr $movie_glx_tcr $rest_gabaplus_tcr $movie_gabaplus_tcr $rest_gaba_tcr $movie_gaba_tcr $rest_glu_tcr $movie_glu_tcr $rest_glu_megapress $rest_vals $rest_fwhm_slaser $movie_glu_megapress $movie_vals $movie_fwhm_slaser $rest_gaba_gabaplus $rest_gaba_gannet \
- $rest_fwhm_osprey $movie_gaba_gabaplus $movie_gaba_gannet $movie_fwhm_osprey $rest_hurst $rest_hurst_old $rest_hurst_vold \
- $rest_hurst_vn $rest_hurst_fullfreq_vn $rest_hurst_fullfreq_blur_vn $rest_hurst_vn_old $rest_hurst_vn_vold $movie_hurst $movie_hurst_old $movie_hurst_vold $movie_hurst_vn \
- $movie_hurst_fullfreq_vn $movie_hurst_fullfreq_blur_vn $movie_hurst_vn_old $movie_hurst_vn_vold $rest_meanfd $movie_meanfd \
- $rest_alff $rest_alff_blur $rest_alff_vn $rest_alff_blur_vn $movie_alff $movie_alff_blur $movie_alff_vn $movie_alff_blur_vn $rest_lc_vals $movie_lc_vals $rest_lc_gaba $movie_lc_gaba \
- $rest_bold $rest_bold_vn $movie_bold $movie_bold_vn $rest_hurst_whole $rest_hurst_fullfreq $rest_hurst_fullfreq_blur $movie_hurst_whole \
- $movie_hurst_fullfreq $movie_hurst_fullfreq_blur $rest_hurst_fullfreq_whole $movie_hurst_fullfreq_whole $rest_glx_gannet $movie_glx_gannet $rest_hurst_dmn $movie_hurst_dmn $rest_hurst_fpn \
- $movie_hurst_fpn $rest_gaba_tcr_gannet $movie_gaba_tcr_gannet >> ${participants_measures}/sub-${subid}_EI_Hurst_Measures_Expanded.tsv
- participants_measures=${BIDSdir}/participant_measures
- echo "Step 9b: Parsing TSV: Old Basis Set"
- echo "Rest_Asc Rest_Asp Rest_Cr Rest_CrCH2 Rest_GABA_slaser Rest_GPC Rest_GSH Rest_Gln Rest_Glu \
- Rest_mI Rest_Lac Rest_NAA Rest_NAAG Rest_PCh Rest_PCr Rest_PE Rest_sI Rest_Tau Rest_MM09 Rest_MM12 Rest_MM14 \
- Rest_MM17 Rest_MM20 Rest_Lip09 Rest_Lip13 Rest_Lip20 Rest_tNAA Rest_Glx Rest_tCho Rest_tCr Rest_FWHM_SLASER \
- Movie_Asc Movie_Asp Movie_Cr Movie_CrCH2 Movie_GABA_slaser Movie_GPC Movie_GSH Movie_Gln \
- Movie_Glu Movie_mI Movie_Lac Movie_NAA Movie_NAAG Movie_PCh Movie_PCr Movie_PE Movie_sI Movie_Tau \
- Movie_MM09 Movie_MM12 Movie_MM14 Movie_MM17 Movie_MM20 Movie_Lip09 Movie_Lip13 Movie_Lip20 Movie_tNAA \
- Movie_Glx Movie_tCho Movie_tCr Movie_FWHM_SLASER Rest_GABA Rest_GABAplus Rest_GABA_Gannet Rest_FWHM_Osprey \
- Movie_GABA Movie_GABAplus Movie_GABA_Gannet Movie_FWHM_Osprey LC_Rest_Ala LC_Rest_AlaSD LC_Rest_Asc LC_Rest_AscSD \
- LC_Rest_Asp LC_Rest_AspSD LC_Rest_Cr LC_Rest_CrSD LC_Rest_GABA_sl LC_Rest_GABASD_sl LC_Rest_GSH LC_Rest_GSHSD \
- LC_Rest_Gln LC_Rest_GlnSD \
- LC_Rest_Glu LC_Rest_GluSD LC_Rest_Lac LC_Rest_LacSD LC_Rest_NAA LC_Rest_NAASD LC_Rest_NAAG LC_Rest_NAAGSD \
- LC_Rest_PCr LC_Rest_PCrSD LC_Rest_PE LC_Rest_PESD LC_Rest_Tau LC_Rest_TauSD LC_Rest_mI LC_Rest_mISD \
- LC_Rest_sI LC_Rest_sISD LC_Rest_CrCH2 LC_Rest_CrCH2SD LC_Rest_CrPCr LC_Rest_CrPrSD LC_Rest_NAANAAG \
- LC_Rest_NAANAAGSD LC_Rest_GluGln LC_Rest_GluGlnSD LC_Rest_Lip13a LC_Rest_Lip13aSD LC_Rest_Lip13b LC_Rest_Lip13bSD \
- LC_Rest_Lip09 LC_Rest_Lip09SD LC_Rest_MM09 LC_Rest_MM09SD LC_Rest_Lip20 LC_Rest_Lip20SD LC_Rest_MM20 \
- LC_Rest_MM20SD LC_Rest_MM12 LC_Rest_MM12SD LC_Rest_MM14 LC_Rest_MM14SD LC_Rest_MM17 LC_Rest_MM17SD \
- LC_Movie_Ala LC_Movie_AlaSD LC_Movie_Asc LC_Movie_AscSD LC_Movie_Asp LC_Movie_AspSD LC_Movie_Cr LC_Movie_CrSD \
- LC_Movie_GABA_sl LC_Movie_GABASD_sl LC_Movie_GSH LC_Movie_GSHSD LC_Movie_Gln LC_Movie_GlnSD \
- LC_Movie_Glu LC_Movie_GluSD LC_Movie_Lac LC_Movie_LacSD LC_Movie_NAA LC_Movie_NAASD LC_Movie_NAAG \
- LC_Movie_NAAGSD LC_Movie_PCr LC_Movie_PCrSD LC_Movie_PE LC_Movie_PESD LC_Movie_Tau LC_Movie_TauSD \
- LC_Movie_mI LC_Movie_mISD \
- LC_Movie_sI LC_Movie_sISD LC_Movie_CrCH2 LC_Movie_CrCH2SD LC_Movie_CrPCr LC_Movie_CrPrSD LC_Movie_NAANAAG \
- LC_Movie_NAANAAGSD LC_Movie_GluGln LC_Movie_GluGlnSD LC_Movie_Lip13a LC_Movie_Lip13aSD LC_Movie_Lip13b \
- LC_Movie_Lip13bSD \
- LC_Movie_Lip09 LC_Movie_Lip09SD LC_Movie_MM09 LC_Movie_MM09SD LC_Movie_Lip20 LC_Movie_Lip20SD LC_Movie_MM20 \
- LC_Movie_MM20SD LC_Movie_MM12 LC_Movie_MM12SD LC_Movie_MM14 LC_Movie_MM14SD LC_Movie_MM17 LC_Movie_MM17SD \
- LC_Rest_GABA \
- LC_Rest_GABASD LC_Movie_GABA LC_Movie_GABASD Rest_Hurst_DMN Movie_Hurst_DMN Rest_Hurst_FPN Movie_Hurst_FPN" > ${participants_measures}/sub-${subid}_EI_Hurst_Old_Basis.tsv
- echo $rest_vals_old $rest_fwhm_slaser_old $movie_vals_old $movie_fwhm_slaser_old $rest_gaba_gabaplus_old $rest_gaba_gannet_old \
- $rest_fwhm_osprey_old $movie_gaba_gabaplus_old $movie_gaba_gannet_old $movie_fwhm_osprey_old \
- $rest_lc_vals_old $movie_lc_vals_old $rest_lc_gaba_old $movie_lc_gaba_old $rest_hurst_dmn $movie_hurst_dmn $rest_hurst_fpn $movie_hurst_fpn >> ${participants_measures}/sub-${subid}_EI_Hurst_Old_Basis.tsv
- echo "Step 9: Complete"
- ## Collecting PSD Data
- echo "#------------------------------------------------------------------#"
- echo "Step 10: Collecting PSD Data"
- psd_data=${participants_measures}/psd_data
- [ -d $psd_data ] || mkdir -p $psd_data
- [ -d ${psd_data}/rest ] || mkdir -p ${psd_data}/rest
- [ -d ${psd_data}/movie ] || mkdir -p ${psd_data}/movie
- cp ${hurstdir}/sub-${subid}/rest/psd_fullfreq.csv $psd_data/rest/psd_fullfreq_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/rest/psd_0.013_to_0.151.csv $psd_data/rest/psd_0.013_to_0.151_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/rest/psd_0.013_to_0.101.csv $psd_data/rest/psd_0.013_to_0.101_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/rest/psd_0.013_to_0.088.csv $psd_data/rest/psd_0.013_to_0.088_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/movie/psd_fullfreq.csv $psd_data/movie/psd_fullfreq_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/movie/psd_0.013_to_0.151.csv $psd_data/movie/psd_0.013_to_0.151_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/movie/psd_0.013_to_0.101.csv $psd_data/movie/psd_0.013_to_0.101_sub-${subid}.csv
- cp ${hurstdir}/sub-${subid}/movie/psd_0.013_to_0.088.csv $psd_data/movie/psd_0.013_to_0.088_sub-${subid}.csv
- echo "Step 10: Complete"
- ## Collecting Additional Hurst Data
- echo "#------------------------------------------------------------------#"
- echo "Step 11: Collecting Additional Hurst Data"
- add_hurst_data=${participants_measures}/add_hurst_data
- [ -d $add_hurst_data ] || mkdir -p $add_hurst_data
- [ -d ${add_hurst_data}/rest ] || mkdir -p ${add_hurst_data}/rest
- [ -d ${add_hurst_data}/movie ] || mkdir -p ${add_hurst_data}/movie
- rest_bp_dfa=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_bp_dfa_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_bp_ksghe=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_bp_ksghe_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_bp_rs=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_bp_rs_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_bp_wavelet=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_bp_wavelet_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_bp_welch_gt1div=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_bp_welch_gt1div_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_BP_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_dfa=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_dfa_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_ksghe=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_ksghe_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_rs=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_rs_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_wavelet=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_wavelet_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- rest_welch_gt1div=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/gm_voxel_mask_space-BOLD_rest.nii.gz -M`
- rest_welch_gt1div_vn=`fslstats ${hurstdir}/sub-${subid}/rest/additional/desc-optcomDenoised_bold_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/rest/Yeo2011_Visual_to_native.nii.gz -M`
- movie_bp_dfa=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_bp_dfa_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_bp_ksghe=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_bp_ksghe_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_bp_rs=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_bp_rs_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_bp_wavelet=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_bp_wavelet_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_bp_welch_gt1div=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_bp_welch_gt1div_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_BP_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_dfa=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_dfa_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_DFAHurst_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_ksghe=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_ksghe_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_KSGHE_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_rs=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_rs_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_RS_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_wavelet=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_wavelet_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_Wavelet_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- movie_welch_gt1div=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/gm_voxel_mask_space-BOLD_movie.nii.gz -M`
- movie_welch_gt1div_vn=`fslstats ${hurstdir}/sub-${subid}/movie/additional/desc-optcomDenoised_bold_Welch-gt1div_Hurst.nii.gz -n -k ${hurstdir}/sub-${subid}/movie/Yeo2011_Visual_to_native.nii.gz -M`
- echo "Rest_BP_DFA Rest_BP_DFA_VN Rest_BP_KSGHE Rest_BP_KSGHE_VN Rest_BP_RS Rest_BP_RS_VN \
- Rest_BP_Wavelet Rest_BP_Wavelet_VN Rest_BP_Welch_GT1Div Rest_BP_Welch_GT1Div_VN \
- Rest_DFA Rest_DFA_VN Rest_KSGHE Rest_KSGHE_VN Rest_RS Rest_RS_VN \
- Rest_Wavelet Rest_Wavelet_VN Rest_Welch_GT1Div Rest_Welch_GT1Div_VN \
- Movie_BP_DFA Movie_BP_DFA_VN Movie_BP_KSGHE Movie_BP_KSGHE_VN Movie_BP_RS Movie_BP_RS_VN \
- Movie_BP_Wavelet Movie_BP_Wavelet_VN Movie_BP_Welch_GT1Div Movie_BP_Welch_GT1Div_VN \
- Movie_DFA Movie_DFA_VN Movie_KSGHE Movie_KSGHE_VN Movie_RS Movie_RS_VN \
- Movie_Wavelet Movie_Wavelet_VN Movie_Welch_GT1Div Movie_Welch_GT1Div_VN" > ${add_hurst_data}/sub-${subid}_Add_Hurst.tsv
- echo $rest_bp_dfa $rest_bp_dfa_vn $rest_bp_ksghe $rest_bp_ksghe_vn $rest_bp_rs $rest_bp_rs_vn \
- $rest_bp_wavelet $rest_bp_wavelet_vn $rest_bp_welch_gt1div $rest_bp_welch_gt1div_vn \
- $rest_dfa $rest_dfa_vn $rest_ksghe $rest_ksghe_vn $rest_rs $rest_rs_vn \
- $rest_wavelet $rest_wavelet_vn $rest_welch_gt1div $rest_welch_gt1div_vn \
- $movie_bp_dfa $movie_bp_dfa_vn $movie_bp_ksghe $movie_bp_ksghe_vn $movie_bp_rs $movie_bp_rs_vn \
- $movie_bp_wavelet $movie_bp_wavelet_vn $movie_bp_welch_gt1div $movie_bp_welch_gt1div_vn \
- $movie_dfa $movie_dfa_vn $movie_ksghe $movie_ksghe_vn $movie_rs $movie_rs_vn \
- $movie_wavelet $movie_wavelet_vn $movie_welch_gt1div $movie_welch_gt1div_vn \ >> ${add_hurst_data}/sub-${subid}_Add_Hurst.tsv
- echo "Step 11: Complete"
- ## Collecting Gannet Spectrum Data
- echo "#------------------------------------------------------------------#"
- echo "Step 12: Collecting Gannet Spectrum Data"
- gannet_spectrum_data=${BIDSdir}/gannet_spectrum_data
- [ -d $gannet_spectrum_data ] || mkdir -p $gannet_spectrum_data
- [ -d ${gannet_spectrum_data}/rest ] || mkdir -p ${gannet_spectrum_data}/rest
- [ -d ${gannet_spectrum_data}/movie ] || mkdir -p ${gannet_spectrum_data}/movie
- cp ${MEGAPRESSdir}/sub-${subid}/rest_gannet/gannetout.mat $gannet_spectrum_data/rest/gannet_spectrum_sub-${subid}.mat
- cp ${MEGAPRESSdir}/sub-${subid}/movie_gannet/gannetout.mat $gannet_spectrum_data/movie/gannet_spectrum_sub-${subid}.mat
- echo "Step 12: Complete"
- echo "Finished sub-${subid}"
- done
06_MeasureCollection.sh at commit 3e7f845, no license · at the source
Overview
- School of Biomedical Engineering, The University of British Columbia, Vancouver, BC, Canada
- Department of Radiology, Weill Cornell Medicine, New York, NY, United States
- BC Children’s Hospital Research Institute, The University of British Columbia, Vancouver, BC, Canada
- Pediatrics, The University of British Columbia, Vancouver, BC, Canada
Abstract
A 3T multimodal MR study of healthy adults (n = 18; 10 female; 21.3–53.4 years) was performed to investigate the relationship between functional magnetic resonance imaging (fMRI) long-range temporal correlations and excitatory/
Reproduced under the paper's license (CC BY), from the paper cited above.
Repository
Its files are read in the Code ↔ Paper reader above, with 9 matches between paragraphs and lines of code.
WeberLab/EI_Hurst_Analysis
3e7f84551f18f19b9361a0586b1ff5a0f0231687, 17 February 2025Availability: 1 check, the latest on 28 September 2026: the link answers
- 28 September 2026: the link answers
33 files
- .ipynb_checkpoints/
OliviaScriptTest-checkpo , Jupyter, 1 lineint.ipynb - .ipynb_checkpoints/
PSD_Analysis-checkpoint. , Jupyter, 134 linesipynb - 01_BashPipeline.sh, Shell, 747 lines, 2 matches
- 02_VisualNetwork.sh, Shell, 59 lines
- 03_7mmSmooth.sh, Shell, 54 lines
- 04_fALFF_7mmSmooth.sh, Shell, 48 lines
- 05_MarkBasisSet.sh, Shell, 375 lines
- 06_MeasureCollection.sh, Shell, 371 lines, 2 matches
- 07_TestNewOsprey.sh, Shell, 67 lines
- AverageVoxelLocation.sh, Shell, 34 lines
- BIDS-to-Methods/
bids-to-methods.ipynb , Jupyter, 131 lines - BIDS-to-Methods/
bids-to-methods.py , Python, 164 lines - BIDS_Pilot01.sh, Shell, 68 lines
- EI_Hurst_Analysis_Prelim
.R , R, 179 lines - EI_Hurst_Measure_Collect
ion.sh , Shell, 58 lines - EI_Hurst_Measure_Collect
ion_Expanded.sh , Shell, 135 lines, 1 match - EI_Hurst_Measure_Collect
ion_Expanded_AW.sh , Shell, 38 lines - GannetSpectrum.m, MATLAB, 55 lines
- OliviaScriptTest.ipynb, Jupyter, 114 lines
- PSD_Analysis.ipynb, Jupyter, 231 lines
- dfa/
dfa.py , Python, 101 lines, 1 match - dfa/
generate.py , Python, 74 lines - dfa/
notebook/ , Jupyter, 296 lines, 1 matchdfa_simulations.ipynb - gannetauto.m, MATLAB, 53 lines, 2 matches
- grab_results.sh, Shell, 24 lines
- renameseries.sh, Shell, 23 lines
- test_CT.sh, Shell, 11 lines
- varioush.py, Python, 173 lines
- welch.py, Python, 181 lines
- welch_freq_sel.py, Python, 224 lines
- welch_freq_sel_v2.py, Python, 226 lines
- welch_fullfreq.py, Python, 201 lines
- README.md, Text, 35 lines
The paper's code and data availability statement is in the Data section.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 32 scripts, each with its path and the digest of its content;
- 9 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
No dataset and no data link were found in the paper.
Data and Code Availability
All code used in this paper is available at github.com/
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 2, 28 September 2026
- Authors: added Jessica Archibald (0000-0001-6651-183X); removed Jessica Archibald
Version 1, 28 September 2026: the first record
Recorded: type, language, journal, volume, pages, dates, 3 authors, 9 keywords, 139 references, 3 RRIDs.
Cite
This paper
Sochan, L., Archibald, J., & Weber, A. M. (2026). Does the brain's E:I balance really shape long-range temporal correlations? Lessons learned from 3T MRI. Imaging neuroscience (Cambridge, Mass.), 4, IMAG.a.1252. https://
BibTeX
@article{sochan2026does,
author = {Sochan, Lydia and Archibald, Jessica and Weber, Alexander Mark},
title = {{Does the brain's E:I balance really shape long-range temporal correlations? Lessons learned from 3T MRI}},
journal = {Imaging neuroscience (Cambridge, Mass.)},
year = {2026},
month = may,
volume = {4},
pages = {IMAG.a.1252},
publisher = {MIT Press},
issn = {2837-6056},
doi = {10.1162/
url = {https://
pmid = {42232071},
pmcid = {PMC13224315}
}
RIS
TY - JOUR
AU - Sochan, Lydia
AU - Archibald, Jessica
AU - Weber, Alexander Mark
TI - Does the brain's E:I balance really shape long-range temporal correlations? Lessons learned from 3T MRI
T2 - Imaging neuroscience (Cambridge, Mass.)
J2 - Imaging Neurosci (Camb)
PY - 2026
DA - 2026/
VL - 4
SP - IMAG.a.1252
SN - 2837-6056
PB - MIT Press
DO - 10.1162/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1162/
"type": "article-journal",
"title": "Does the brain's E:I balance really shape long-range temporal correlations? Lessons learned from 3T MRI",
"container-title": "Imaging neuroscience (Cambridge, Mass.)",
"author": [
{
"family": "Sochan",
"given": "Lydia"
},
{
"family": "Archibald",
"given": "Jessica"
},
{
"family": "Weber",
"given": "Alexander Mark"
}
],
"container-title-short":
"volume": "4",
"page": "IMAG.a.1252",
"DOI": "10.1162/
"PMID": "42232071",
"PMCID": "PMC13224315",
"ISSN": "2837-6056",
"publisher": "MIT Press",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
5,
29
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1162/netn.a.547 [code]
- An evaluation of the efficacy of single-echo and multi-echo fMRI denoising strategies.Journal: Network neuroscience (Cambridge, Mass.)In common: fMRIPrep, tedana, ANTs, 9 other tools, fMRI, 9 references
- [2] doi:10.1162/imag.a.1347 [code]
- Neural and behavioural correlates of theory of mind reasoning in five-year-old children born preterm.Journal: Imaging neuroscience (Cambridge, Mass.)In common: MRIQC, Dcm2Bids, fMRIPrep, 12 other tools, fMRI, 1 reference
- [3] doi:10.1016/j.dcn.2026.101765 [code]
- Fusiform face area development correlates with development in higher-order social brain regions.Journal: Developmental cognitive neuroscienceIn common: MRIQC, Dcm2Bids, fMRIPrep, 12 other tools, fMRI, 1 reference
- [4] doi:10.1038/s41467-026-71151-2 [code]
- Common and distinct neural correlates of social interaction processing and theory of mind in narratives.Journal: Nature communicationsIn common: fMRIPrep, AFNI, ANTs, 9 other tools, 7 references
- [5] doi:10.1162/imag.a.1198 [code]
- MEPrep: A robust pipeline for multi-echo fMRI denoising and preprocessing.Journal: Imaging neuroscience (Cambridge, Mass.)In common: fMRIPrep, tedana, AFNI, 8 other tools, fMRI, 6 references
- [6] doi:10.1002/mrm.70380 [code]
- The Impact and Reliability of Tissue Segmentation on In Vivo Magnetic Resonance Spectroscopy Metabolite Quantification.Journal: Magnetic resonance in medicineIn common: structural MRI / diffusion, 8 references, author Jessica Archibald
- [7] doi:10.1162/imag.a.1275 [code]
- Glucose metabolism echoes long-range temporal correlations in the human brain.Journal: Imaging neuroscience (Cambridge, Mass.)In common: NiBabel, scikit-learn, pandas, 3 other tools, fMRI, 8 references
- [8] doi:10.1016/j.neuron.2026.04.011 [code]
- Precision fMRI reveals densely interdigitated network patches with conserved motifs in the lateral prefrontal cortex.Journal: NeuronIn common: tedana, AFNI, ANTs, 7 other tools, fMRI, 4 references
- [9] doi:10.1002/hbm.70621 [code]
- Optimising 7T-fMRI for Imaging Regions of Magnetic Susceptibility.Journal: Human brain mappingIn common: fMRIPrep, tedana, AFNI, 2 other tools, fMRI, 6 references
- [10] doi:10.1038/s41467-026-73668-y [code]
- Convergent and divergent brain-cognition development in early adolescence.Journal: Nature communicationsIn common: tedana, AFNI, ANTs, 10 other tools, fMRI, 1 reference
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 1 repository of the authors' code, each at its verified commit and with its license, 32 scripts, and 9 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:ab5484ed6ab83aac…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
