#!/bin/tcsh

# ----------------------------------------------------------------------
# get an idea of where this has an effect on the variance by computing
# the ratio of (variance without) to (variance with)

# perform the computation in the results directory from script s1
cd s620.1.RICOR.perrun.results

3dTstat -sos -prefix sos.1.perrun errts.s620.1.RICOR.perrun+orig

# sos already computed
# 3dTstat -sos -prefix ./sos.2.base               \
#         ../s620.2.base.results/errts.s620.2.base+orig

# take a simple ratio of variances (equal to ratio of sum of squares)
3dcalc -a sos.1.perrun+orig -b ../sos.results/sos.2.base+orig  \
       -expr 'step(a)*b/a' -prefix ./sos.ratio