Hello,
I am trying to run 2 correlations using 3dRegAna. I am trying to correlate individual percent signal change maps (averaged across 2 runs) with 2 separate behavioral scores. Here is the script I am using:
######## Correlation of extaversion ########################
3dRegAna \
-rows 16 \
-cols 1 \
-xydata 22 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj2_Hap_Neu+tlrc \
-xydata 25 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj3_Hap_Neu+tlrc \
-xydata 30 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj4_Hap_Neu+tlrc \
-xydata 28 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj6_Hap_Neu+tlrc \
-xydata 35 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj11_Hap_Neu+tlrc \
-xydata 12 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj12_Hap_Neu+tlrc \
-xydata 26 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj13_Hap_Neu+tlrc \
-xydata 28 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj14_Hap_Neu+tlrc \
-xydata 39 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj15_Hap_Neu+tlrc \
-xydata 31 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj16_Hap_Neu+tlrc \
-xydata 31 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj17_Hap_Neu+tlrc \
-xydata 38 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj18_Hap_Neu+tlrc \
-xydata 29 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj20_Hap_Neu+tlrc \
-xydata 30 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj21_Hap_Neu+tlrc \
-xydata 30 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj22_Hap_Neu+tlrc \
-xydata 25 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj24_Hap_Neu+tlrc \
-rmsmin 1.0 \
-fdisp 2 \
-model 1 : 0 \
-bucket 0 percent_STATE_STAI_Hap_Neu
####Hap-Fix
3dRegAna \
-rows 16 \
-cols 1 \
-xydata 28 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj2_Hap_Fix+tlrc \
-xydata 28 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj3_Hap_Fix+tlrc \
-xydata 28 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj4_Hap_Fix+tlrc \
-xydata 31 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj6_Hap_Fix+tlrc \
-xydata 31 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj11_Hap_Fix+tlrc \
-xydata 34 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj12_Hap_Fix+tlrc \
-xydata 34 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj13_Hap_Fix+tlrc \
-xydata 34 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj14_Hap_Fix+tlrc \
-xydata 30 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj15_Hap_Fix+tlrc \
-xydata 24 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj16_Hap_Fix+tlrc \
-xydata 34 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj17_Hap_Fix+tlrc \
-xydata 39 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj18_Hap_Fix+tlrc \
-xydata 33 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj20_Hap_Fix+tlrc \
-xydata 37 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj21_Hap_Fix+tlrc \
-xydata 37 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj22_Hap_Fix+tlrc \
-xydata 19 /study/hnmri/processed/group_analysis_hnmri/percent_change/avg_runs_percent_subj24_Hap_Fix+tlrc \
-rmsmin 1.0 \
-fdisp 2 \
-model 1 : 0 \
-bucket 0 percent_extra_Hap_Fix
I have tried to run this script multiple times and it appears that 3dRegAna is executed twice, but there is no output or error message. This all happens within about 1 second. Here are the command lines I get:
[somerville@cayenne hnmri_analysis_scripts]$ sh percent_corr_extra_STAIState
Program: 3dRegAna
Author: B. Douglas Ward
Initial Release: 10 Oct 1997
Latest Revision: 15 Aug 2001
Program: 3dRegAna
Author: B. Douglas Ward
Initial Release: 10 Oct 1997
Latest Revision: 15 Aug 2001
[somerville@cayenne hnmri_analysis_scripts]
Basically, this script gave me no output but also no error message, so I don't know what is wrong. I'd appreciate any input you have!
Thanks,
Leah Somerville
Lab of Paul Whalen, Ph.D.
University of Wisconsin