Dear AFNI experts,
I was doing pretest for SEM analysis with time series (250 time points) from 3 ROIs at each subject (2 subjects) So I run SEM_script using a command like below.
tcsh -x AFNI_SEM_script_treegrowth.csh c03.1D c05.1D thetas.1D
Unfortunately, I face an error message at very begining of process of SEM.
thetas.1D is below
# RIFG RIFG2 RIns
RIFG 0 2 2
RIFG2 2 0 2
RIns 2 2 0
#One of timeseries files,c03.1D, for 2 subjects is like below.
7465.09 7273.43 11149.5
7556.84 7364.05 11151.7
7434.16 7281.24 11201.4
7377.94 7332.38 11123.7
7451.46 7273.34 11150.6
7530.23 7252.7 11251.2
7490.49 7288.49 11198.1
.....
# Here is the error message
[root@gwd TimeSeries]# tcsh -x AFNI_SEM_script_treegrowth.csh c03.1D c05.1D thetas.1D
@ last = 3 - 1
set minModel = thetas.1D
set subj_dsets = ( c03.1D c05.1D )
set NSUBJS = 2
glob c03.1D c05.1D
c03.1Dc05.1Decho
set NVALS = `1dcat ${subj_dsets[1]}'[0]' | wc -l`
1dcat c03.1D[0]
wc -l
set NROIS = `1dcat ${subj_dsets[1]}'{0}' | wc -w`
1dcat c03.1D{0}
wc -w
@ NP = 250 + 2
@ NVM1 = 250 - 2
rm: No match.
rm: No match.
rm: No match.
rm: No match.
rm: No match.
rm: No match.
rm -i -f SEM_resvar.1D
rm -i -f SEM_resvarcorr.1D
rm: No match.
rm: No match.
rm: No match.
rm -i -f SEM_autocorr.1D
rm -i -f SEM_cov.1D
rm -i -f SEM_corr.1D
rm -i -f SEM_temp.1D
rm -i -f SEM_DF.1D
foreach roi ( `count 1 $NROIS` )
count 1 3
@ roi_index = 0001 - 1
foreach subj ( `count 1 $NSUBJS` )
count 1 2
1dcat c03.1D[0]
end
1dcat c05.1D[0]
end
1dcat SEM_roi0001_subj0001_extract.1D SEM_roi0001_subj0002_extract.1D
1dsvd SEM_roi0001_allsubjs.1D
grep -A 252 -m 1 ++ Left Vectors SEM_svd_roi0001.txt
tail -250
tr -d :
1dcat SEM_temp.1D[1]
grep -A 1 -m 1 ++ Left Vectors SEM_svd_roi0001.txt
tail -1
1dtranspose SEM_temp.1D
1deval -a SEM_sv_roi0001.1D -expr a*a
3dTstat -sum -prefix SEM_resvar_roi0001 SEM_sv2_roi0001.1D{1..$}'
++ 3dTstat: AFNI version=AFNI_2007_05_29_1644 (Aug 10 2007) [64-bit]
++ Authored by: KR Hammett & RW Cox
** FATAL ERROR: Can't use dataset with < 2 values per voxel!
1dcat SEM_resvar_roi0001.1D
** Can't read input file SEM_resvar_roi0001.1D
Thank you,
Sean
