OK. There are a couple "attributes" to add to the dset to attach the statistical information for the p-to-t calcs. These can be aded in a couple different ways, this (below) seems the easier one than putting in individual attributes.

Let's say you have a dset of 2 subbricks, where the [0]th brick is an effect estimate, and the [1]th brick is a t-stat with 412 degrees of freedom (DOFs). Then the following would take some DSET and populate the [1]th brick with the desired stat information (and put some reference labels for the strings, for good measure):

set DSET      = some_data_set_name.nii.gz
set stat_lab  = fitt
set dof       = 412

3drefit                                  \
    -sublabel 0 EFFECT_ESTIMATE \
    -sublabel 1 STAT_ESTIMATE \
    -substatpar 1 ${stat_lab}  ${dof}  \
    ${DSET}

# optionally, add FDR curve info;  see "-FDRmask .." and "-STATmask .." options for this, as well
3drefit                  \
    -addFDR        \
    ${DSET}

How does that seem?

-pt



Edited 1 time(s). Last edit at 11/08/2022 09:01AM by ptaylor.