Hi - indeed - actually it's slightly simpler than that - see the FDR
manual page:
http://www.fmrib.ox.ac.uk/fsl/randomise/fdr.html
Cheers, Steve.
On 8 Jun 2006, at 19:39, Greg Scott wrote:
> Greetings fmrib,
>
> I want to run false discovery rate thresholding on TBSS/Randomise
> output and
> came up with the analysis steps mentioned below. When you have the
> time,
> could someone double-check the following. Thank you yet again. Greg
>
>
> After TBSS and Randomise using default parameters.
>
> 1.) Turn unthresholded tstat into logp image:
>
> avwmaths tbss_tstat1 -mul 0 -add 1 dummyvar
> ttologp dummyvar tbss_tstat1 1000
> avwmaths logps -exp ptstat
>
> 2.) FDR:
>
> fdr -i ptstat -m mean_FA_skeleton_mask -q 0.05
>
> 3.) Threshold p values according to FDR output
>
> avwmaths ptstat -mul -1 -add 1 -thr {(1 - FDR OUTPUT)} resultname
------------------------------------------------------------------------
---
Stephen M. Smith, Professor of Biomedical Engineering
Associate Director, Oxford University FMRIB Centre
FMRIB, JR Hospital, Headington, Oxford OX3 9DU, UK
+44 (0) 1865 222726 (fax 222717)
[log in to unmask] http://www.fmrib.ox.ac.uk/~steve
------------------------------------------------------------------------
---
|