Try this:
1) Find the degrees of freedom for the t-stat sub-brick using a command like this (assuming the sub-brick number is 9):
3dAttribute BRICK_STATAUX myFile+tlrc'[9]'
The last number is the degrees of freedom.
2) Compute the critical value corresponding to the p-value of, for example, 0.05 (DF is the degrees of freedom from 1) above):
cdf -p2t fitt 0.05 DF
3) Find out the upper bound of 95% uncertainty interval (assuming the critical value from 2) is CR and the beta and t-stat sub-bricks are 8 and 9):
3dcalc -a myFile+tlrc'[8]' -b myFile+tlrc'[9]' -expr 'a+CR*a/b' -prefix upper
and lower bound
3dcalc -a myFile+tlrc'[8]' -b myFile+tlrc'[9]' -expr 'a-CR*a/b' -prefix lower
Gang
Edited 1 time(s). Last edit at 02/03/2022 04:42PM by Gang.