For opening all scripts we used in virtual subway navigation study. Here, we are grateful for all the supports given by Peter Zeidman in DCM analysis.
We used Psychopy v3.0 to present the whole task.
You can find the python scripts in Task_Program_Psychopy/.
Directory named "formal_experiment" contains the program used in formal experiment when fMRI scanning.
Directory named "train_stage" contains the program used in the training stage a day before scanning.
The linear mixed model which modeling the behavioral data of participants was built using R package lmerTest.
You can find the R scritps in Behavior_analysis/.
Step-1_Pre-analysis_data_resahpe.R is the script to shape the data recorded from the tasks which is located in behav_data_clean/. And subway_behaviordata_assumble_DA.csv is an example of the data after shaping by the pipeline.
Step-2_LMMs_Model_Construction.R is the script for two aims: (1) Estimating the effect of multiple planning complexities to the participants' RT (please see Fig.3a on our paper), and (2) testing the effect synchronicity of the planning complexities for influencing the RT and brain activity by Pearson's correlation (please see Fig.3c on our paper).
GLM1_roi_behavior/ contains the z-map of brain activity estimated by the GLM1, and the ROIs were identified by the DL modulation.
MRIqc_group/ contains the html report for the quality functional and structural images.
To detect the response of brain areas to computational costs, we built a general linear model 1 which contains four costs for parameteric modulation.
The GLM1 was conducted using FSL's feat. All *.fsf files in GLM1_related/ are design files using in run-level, subject-level and group-level.
The bash script BashRun_first_level_GLM1_Feat.sh is used for batching process of run-level GLM1.
The time series could be found in the directory named TimePointExtraction_GLM1.
To investgate the neural architecture underlying hierarchical planning and test the connectivity in dorsomedial frontal cortext contributed to the individual difference in planning efficiency, DCM and PEB were performed. The scripts we ued are placed in the directory DCM_related/.
We first construct a subject-level GLM using SPM. Then we extracted the time series in each VOI, including the dorsomedial prefrontal cortex, the bilateral premotor cortex and the superior parietal cortex. The script used to build GLM and extract VOIs are Step1_GLM_all_pipeline_concate.m and Step2_GLM_pipeline_VOI_extraction.m.
The fist level DCM was specified using Step3_first_level_DCM.m.
The seond level DCM, including PEB analysis, auto-BMR, cross-validation and family-wise model comparison, was estimated using Step4_second_level_DCM.m.
DCM result visualization was conducted using R. The script is Step5_DCM_PEB_visualization.R.
All DCM pipeline was done using SPM12, following the guidance of part 1: First level analysis with DCM for fMRI and part 2: Second level analysis with PEB
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。