AFNI Message Board

History of AFNI updates  

|
November 25, 2002 12:32PM
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
Subject Author Posted

3dRegAna

Leah Somerville November 25, 2002 12:32PM

Re: 3dRegAna

bob cox November 25, 2002 05:16PM

Re: 3dRegAna

Leah Somerville November 25, 2002 05:34PM

Re: 3dRegAna

B. Douglas Ward November 26, 2002 10:05AM

Re: 3dRegAna

Leah Somerville November 26, 2002 12:23PM



Sorry, only registered users may post in this forum.

Click here to login