History of AFNI updates
Suggested reading for beginners
Add yourself to the AFNI map!
|
Author: bob cox (---.nimh.nih.gov)
Date: 11-05-09 10:33
If you want to see the FDR q-value in the AFNI GUI, run the command
3drefit -addFDR ANOVAdataset+tlrc
If you need to use a mask, add the option
-FDRmask maskdataset+tlrc
before the input dataset name.
See the output of '3drefit -help' for a few more details. For lots
more information on FDR in AFNI, see slides 44-53 of the presentation
below:
http://afni.nimh.nih.gov/pub/dist/edu/latest/afni_handouts/afni07_advanced.pdf
http://afni.nimh.nih.gov/pub/dist/edu/latest/afni_handouts/afni07_advanced.pdf
The '3drefit -addFDR' command will compute a curve giving q as a
function of F for each F statistic sub-brick, and then as you change
the F slider threshold, the corresponding q-value is displayed at the
bottom of the slider. In addition, if you hover the mouse cursor over
AFNI's p=/q= label you will see a 'popup hint' that shows the MDF = an
estimate of the 'missed detection fraction' = fraction of true
positives that lie BELOW the current threshold. (MDF=60% is NOT saying
that 60% of the voxels below threshold are true positives; it is
estimating that only 40% of the true positives are above threshold and
the rest are below threshold.)
MDF is sort of the converse to FDR = fraction of positives ABOVE the
current threshold that are false. Estimating MDF reliably is much
harder than FDR, so I wouldn't put a huge amount of reliance on it --
but it can be instructive, since at usual thresholds, MDF is roughly
50-70%, indicating that we might be missing a lot of stuff with our
current FMRI analyses .
3dFDR is less useful -- it is an older program. It replaces the F-
statistic volume with the equivalent z(q)-values. You CAN threshold
with this, but I much prefer thresholding with F (or t) and just
seeing the q-value below the slider.
|
|
|
Helen Weng |
11-04-09 17:54 |
|
Gang Chen |
11-05-09 02:03 |
   |
Re: 3dFDR question on Group statistics from groupana |
|
bob cox |
11-05-09 10:33 |
|
Helen Weng |
11-05-09 14:46 |
|
ziad |
11-05-09 22:35 |
