Hi Gang,
The SEMscript.csh script has the first few lines extracting the eigen values, commented out in the section of the code that computes residual variance. Instead the calculation seems to be done using the first eigen vector. wonder whys that so?
# compute residual variance
# grep -A 1 -m 1 "++ Left Vectors" SEM_svd_roi${roi}.txt | \
# tail -1 > SEM_temp.1D
# 1dtranspose SEM_temp.1D > SEM_sv_roi${roi}.1D
# 1deval -a SEM_sv_roi${roi}.1D -expr 'a*a' > SEM_sv2_roi${roi}.1D
1dcat SEM_temp.1D'[1]' > SEM_sv_roi${roi}.1D
1deval -a SEM_sv_roi${roi}.1D -expr 'a*a' > SEM_sv2_roi${roi}.1D
Thanks,
Pradeep